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ABSTRACT 


This  report  it  a  review  of  currently  avoi table  techniques  for  studying  the  properties  of  objects 
of  interest  in  radar  astronomy,  which  are  usually  spread  radar  targets.  The  report  is  intended 
as  a  sequel  to  Lincoln  Laboratory  Technical  Report  234,  in  which  R.  Price  has  presented  a 
rather  complete  treatment  of  the  theory  and  methodology  of  the  detection  of  such  spread  tar¬ 
gets.'  Here  we  are  interested  in  measurement  (or  estimation)  rather  than  detection. 

Spread  targets  are  defined  to  be  those  that  produce  an  observable  smearing  of  the  echo  in 
range,  or  an  observable  rote  of  echo  fluctuation,  or  both.  (Spread  targets  may  or  may  not  be 
at  the  same  time  extended  targets,  i.e.,  ones  that  fill  the  antenna  beamwidth.) 

The  approach  taken  here  in  presenting  the  available  material  on  spread  target  measurements 
is  to  follow  what  happens  to  the  various  attributes  (amplitude,  delay,  phase,  frequency  shift, 
polarization)  of  an  incident  signal  upon  reflection  from  such  a  target.  Then  the  question  is 
turned  around  by  detailing  the  ways  in  which  the  study  of  rhese  signal  attributes  by  appro¬ 
priate  receiver  processing  can  be  used  to  infer  target  properties.  In  the  case  of  such  radar 
astronomy  targets  as  rotating  planets,  these  target  properties  might  include  range,  velocity, 
shape,  size,  rotation  vector,  and  surface  characteristics.  (The  last  of  these  may  be  studied 
either  as  average  behavior  over  the  targat  surface  or  as  a  function  of  location  on  the  target 
surface.)  The  possible  transmitted  signals  include  simple  (unity  time-bandwidth  product)  sig¬ 
nals,  or  complicated  (large-TW) signals,  including  frequency-spaced  setsof  sinusoids.  The  re¬ 
ceiver  operations  considered  include  processing  of  the  echo  signal  received  at  spaced  receivers 
(interferometry),  and  processing  of  separate  campon  -’.rs  received  at  the  same  point  (polar- 
imetry),  as  well  as  the  much  more  completely  understood  case  of  processing  of  the  output  of  a 
single  receiving  antenna  by  various  methods. 

The  concluding  section  is  a  summary  of  presently  available  information  on  the  variance  of  the 
error  in  making  such  measurements  (including  some  very  rer.«.,ii  results  obtained  by  M.  J.  Levin 
and  R.  "rice). 

Much  c-f  this  report  consists  of  previous  results  that  have  not  been  placed  in  a  connected  story 
before.  One  topic  that  is  new  here  is  the  extension  of  Manosse's  study  of  radar  interferom¬ 
etry  to  include  the  effected  target  rotation,  and  the  resulting  motion  of  the  diffraction  pattern 
observed  at  the  spaced  rodar  receivers  on  earth. 
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RADAR  ASTRONOMY  MEASUREMENT  TECHNIQUES 


I.  INTRODUCTION 

This  report  deals  with  the  use  of  radar  systems  to  measure  quantities  of  astronomical 
interest.  Radar  astronomy  targets,  which  may  be  rotating  solid  bodies,  moving  ionized  regions, 
or  clouds  of  particles,  are  examples  of  spread  radar  targets,  i.e.,  objects  that  impose  on  the 
probing  radar  signal  either  or  both  of  the  following  two  effects:  (a)  an  observable  range  smear¬ 
ing,  and  (b)  an  observable  fluctuation  rate.  It  is  easy  to  imagine  other  examples  of  spread  tar¬ 
gets,  for  example  aircraft,  tumbling  nose  cones,  missile  wakes,  and  so  forth,  to  which  the  ideas 
presented  here  may  be  more  or  less  relevant,  but  we  shall  ignore  these  possibilities  here  and 
talk  in  terms  of  radar  astronomy  targets. 

The  theory  and  practice  of  optimizing  the  detectability  of  spread  targets  has  been  treated 
rather  completely  by  Price1  in  Lincoln  Laboratory  Technical  Report  234,  to  which  this  report  is 
a  sequel.  The  procedures  developed  there  are  applicable  interchangeably  to  the  detection  of 
spread  objects  as  radar  targets  or  the  use  of  them  as  communication  channels.  In  this  report 
we  ignore  the  radar  detection  as  well  as  the  communication  aspects  of  the  spread  target  problem 
and  deal  solely  with  the  use  of  the  radar  system  for  making  measurements  of  the  properties  of 
the  target. 

The  outgoing  and  returning  electromagnetic  signals  are  characterized  by  certain  quantities, 
and  we  would  like  to  investigate  how  the  parameters  of  the  outgoing  signal  are  altered  in  the  re¬ 
flection  process,  and  then  how  this  alteration  may  be  studied.  If  we  regard  the  electrical  signal 
as  a  scalar  quantity,  i.e.,  a  narrow-band  time  function  appearing  at  the  two  terminals  of  the  an¬ 
tenna,  then  a  description  of  the  amplitude,  frequency,  and  phase  as  a  function  of  time  completely 
specifies  the  signal.  If  in  addition  we  are  willing  to  take  into  account  the  vector  direction  of  the 
electrical  signals  (that  is,  the  instantaneous  plane  of  polarization  as  a  function  of  time),  and  the 
way  in  which  all  the  quantities  just  listed  vary  from  point  to  point  in  space,  still  more  sensitive 
measurements  on  the  target  can  be  made.  As  we  shall  see,  certain  interesting  things  may  be 
learned  from  such  interferometric  measurements  (correlation  between  the  waves  ooserved  at 
two  different  points)  and  polarlmetric  measurements  (correlation  between  components  observed 
at  the  same  point). 

A  complete  understanding  of  all  these  factors  embraces  a  very  large  area  indeed,  and  only 
parts  of  the  subject  are  well  understood.  What  we  shall  be  doing  in  this  report  to  systematize 
these  scattered  notions  can  be  stated  this  way:  The  radar  signal  that  is  transmitted  can  be 
written 

x(t)  =  2(t)  cos[2wfjt  +  v>(t)l  ,  (1) 
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which  is  seen  to  be  characterized  by  (a)  scalar  amplitude  E,  (b)  origin  of  the  time  scale  t, 

(c)  phase  <p,  (d)  average  rate  of  change  of  phase,  i.e.,  frequency  2*^,  and  (e)  instantaneous 

direction  of  polarization  x/x,  (The  arrow  indicates  a  vector  quantity.)  Corresponding  to  the 
— * 

transmitted  signal  x(t)  there  is  received  at  a  given  point  an  echo  signal 

w(t)  =  W(t)  cos  [27rfQt  +  ©<t)) 

likewise  characterized  by  five  quantities.  It  is  the  fact  that  these  five  parameters  can  be  con¬ 
trolled  at  will  in  the  transmission  (as  well  as  measured  in  various  ways  on  reception)  that  em¬ 
bodies  any  advantage  that  radar  astronomy  might  have  relative  to  radio  astronomy. 

It  will  be  found  most  useful  to  employ  the  common  representation  of  the  real  signal  x(t),  not 
in  the  trigonometric  form  (Eq.  (1)],  but  rather  as  the  real  part  of  a  complex  signal  &(t) 

x(t)  =  Re  {$?(!)}  =  Re  {X(t)  exp  [j2irfot]}  (2) 

(in  other  words  Re  {&(t)}  =  E(t)  cos  The  quantity  X  (t)  is  called  the  complex  envelope. 

Manipulation  of  TCi t)  is  much  less  cumbersome  than  manipulation  of  x(t)  since  the  usually  irrele¬ 
vant  carrier  frequency  has  been  removed.  B>  exercising  reasonable  care,  the  "real-part-of" 
operation  can  be  reinserted  at  the  conclusion  of  the  manipulation  to  give  the  correct  answer.  The 
hypothetical  imaginary  part  of  2(t)  is  usually  chosen  to  be  the  Hilbert  transform  of  our  real  ob¬ 
servable  waveform  x(t);  then  for  x(t)  a  narrow-band  function  (center  frequency  much  greater 
than  bandwidth)  Im  {  2(t)}  always  represents  the  quadrature  or  out-of -phase  component  of  the  real 
signal  x(t).  Roman  letters  will  be  used  to  indicate  the  real-time  function  and  script  letters  will 
denote  the  corresponding  complex  envelope. 

In  what  follows  we  shall  divide  radar  targets  into  two  classes,  hard  and  soft  targets,  and  in¬ 
vestigate  one  by  one  the  influence  of  the  target  properties  on  each  of  our  attributes: 

amplitude  (or  actually  average  power)  (SecB.  II- B  and  III-A), 

delay  (II-C  and  HI-B), 

phase  and  frequency  (II-D  and  IH-C),  and 

polarization  direction  (II-E  and  III-D). 

After  a  general  discussion  of  the  instrumentation  of  radar  measurements  (Sec.  IV),  we  will  pro¬ 
ceed  in  Sec.  V  to  turn  things  around  and  ask  how  various  target  properties  of  astronomical  in¬ 
terest  may  be  studied  by  studying  the  above  attributes  in  the  received  signal  w(t).  Section  VI 
discusses  errors  in  such  measurements. 

It  would  be  desirable  to  have  a  concise,  and  perfectly  general  notion  of  how  optimally  to 
measure  the  interesting  target  properties  in  the  presence  of  the  inevitable  noise,  finite  observa¬ 
tion  bandwidth,  finite  observing  time,  and  other  factors  causing  errors.  Using  the  statistical 
discipline  known  as  hypothesis  testing,  the  presence  or  absence  of  the  radar  astronomy  target 
has  been  put  on  a  sound  enough  theoretical  basis*  Unfortunately  the  theory  of  measurement  (or, 
in  the  language  of  the  statistician,  estimation  of  target  parameters)  has  been  placed  on  a  similar 
strong  foundation  only  for  certain  very  simple  types  of  measurements.  What  we  shall  attempt 
to  do  here  in  the  case  of  most  of  the  measurement  procedures  is  to  present  them  systematically 
as  a  collection  of  hoc  ideas  about  reasonable  ways  to  measure  things  about  radar  astronomy 
targets,  and,  whenever  possible,  to  give  an  idea  of  the  measurement  accuracy. 
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H.  HARD  TARGETS 

By  hard  targets  we  mean  bodies  that,  although  they  may  be  moving  or  changing  their  orien¬ 
tation,  are  nevertheless  rigid  and  of  fixed  shape.  The  term  soft  targets  will  refer  here  to  those 
that  can  change  shape  with  time.  The  most  Interesting  hard  radar  astronomy  targets  are  ob¬ 
viously  the  moon,  planets,  and  minor  planets.  These  are  all  bodies  whose  size  is  much  larger 
than  any  reasonable  radar  wavelength,  and  which  have  rough  surfaces,  in  the  sense  that  they 
possess  surface  irregularities  having  a  scale  size  of  the  order  of  magnitude  of  a  wavelength. 

With  the  exception  of  some  smaller  minor  planets,  these  hard  targets  may  be  safely  assumed 
to  have  approximately  spherical  shapes  with  rough  surfaces,  and  in  our  treatment  of  hard  targets 
we  will  be  most  interested  in  the  spherical  case.  Where  the  radar  properties  depend  on  a  spher¬ 
ical  target  shape  the  fact  will  be  pointed  out. 

A.  The  Angular  Power  Spectrum  aQ  of  an  Incremental  Portion  of  the  Target 

When  an  incident  plane  wave  strikes  a  rough  surface,  energy  is  reflected  in  various  direc¬ 
tions  {see  Fig.  1).  In  adding  up  the  individual  contributions2  one  finds  that  across  the  wavefront 
of  the  reflection  there  are  variations  in  delay,  phase,  amplitude,  and  direction  of  the  polariza¬ 
tion  vector.  The  distant  observer  receives  a  plane  wave  which  likewise  contains  irregularities 
as  a  function  of  position  laterally  on  the  advancing  plane  wavefront. 


Fig.  1 .  The  manner  In  which  an  irregular  surface  produces 
irregularities  in  the  reflected  wavefront. 


As  we  shall  see,  it  is  possible  to  work  backwards  from  the  statistics  of  these  irregularities 
at  the  distant  observing  point  to  get  the  statistics  of  the  irregularities  at  same  smooth  quasi- 
planar  boundary  near  but  not  on  the  target  surfaces  (dotted  line  of  Fig.  1),  but  the  question  of  how 
these  wavefront  irregularities  relate  to  the  surface  irregularities  themselves  is  a  much  more 
complicated  one.  The  factors  that  enter  include  whether  there  is  an  undulating  or  jagged  (dis¬ 
continuous)  surface,  whether  irregularities  are  smaller  or  greater  than  the  wavelength  in  the 
lateral  direction  and  in  the  up-down  direction,  whether  the  surface  material  is  uniform  in  di¬ 
electric  properties  but  of  a  rough  shape  or  whether  the  body  is  smooth  but  of  irregular  dielectric 
composition,  and  so  forth.  It  is  safe  to  say  that  a  large  number  of  physical  surface  types  could 
give  rise  to  the  same  type  of  irregularities  of  delay,  phase,  amplitude,  and  direction  across 
both  the  dotted  boundary  and  the  received  wavefront.  The  reader  is  referred  to  Ref.  3  for  a  dis¬ 
cussion  of  this  question  of  the  relation  of  terrain  properties  to  properties  of  the  reflected  wavefront. 
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In  this  report  we  shall  be  content  to  start  with  the  irregularities  in  the  reflected  field  produced 
at  the  dotted  boundary  of  Fig.  1  by  terrain  irregularities,  and  not  with  these  latter  irregularities 
themselves.  In  the  rest  of  this  report  when  we  speak  of  the  target  surface  we  shall  mean,  strictly 
speaking,  the  wavefront  of  the  reflected  signal  as  examined  over  a  quasi-planar  region  near  to 
but  not  actually  on  the  surface  of  the  body. 

The  radar  receiver  at  a  great  distance  from  the  target  is  exposed  to  an  electromagnetic 
signal  formed  by  the  propagation  outward  of  the  wavefront  with  all  its  irregularities.  As  we  shall 
see  in  the  next  several  sections,  radar  measurements  on  hard  targets,  whether  performed  by  a 
single  receiver  or  a  set  of  spaced  receivers  (interferometer)  often  amount  simply  to  the  meas¬ 
urement  of  a  quantity  vQ(.<p)  called  the  radar  angular  power  spectrum,  which  expresses  the  sta¬ 
tistics  of  the  wavefront  irregularities  in  a  particular  way. 

In  the  next  several  pages,  we  shall  take  the  familiar  Fourier  transform  relationship  between 
field  intensity  across  an  aperture  (at  the  target)  and  the  remotely  received  signal,  and  derive  an 
expression  for  <^0(<p)  that  shows  just  what  statistical  properties  of  the  wavefront  irregularities  it 
measures.  While  we  are  at  it,  we  shall  define  several  other  terms  and  draw  several  conclusions 
that  will  be  useful  later  on. 


Fig.  2.  Coordinate  syjfem  for  derivation  of  <r„($),  where  $  =  the  angle 
that  an  incident  ray  makes  with  the  normal  to  a  plane  target  region 
(aperture),  $'  =  similar  angle  for  ray  to  receiver,  0  =  equatorial  angle 
measured  in  the  aperture  plane  between  the  projections  of  the  two  rays. 

Figure  2  shows  at  z  =  0  a  plane  region  of  the  target  surface  (in  the  sense  just  mentioned) 
with  the  irregular  pattern  of  the  electric  field  intensity  of  the  echo  being  £  (x,  y,  o).  Denote  the 
component  of  £  that  lines  up  with  the  receiving  antenna  on  the  ground  as  €.  The  script  letter 
indicates  a  complex  envelope  as  before,  and  the  distances  x,  y,  and  z  are  expressed  in  multiples 
of  wavelengths.  For  the  moment  let  us  assume  that  £  is  non-time-varying.  The  well-known 
relationship  giving  the  electric  field  ~Yf  at  the  distant  receiver  point  (xQ,  y  ,  zQ)  on  the  ground 
plane  as  the  Fourier  transform  of  £  is4 

m)  =  ^  £(x,y,  o)  exp[-j2ir(fx  +  my)]  dx  dy 

Aperture 
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(3) 


or,  more  compactly, 

'W'(l.m)  <=  f(x,y)  , 
where  /,  m,  and  n  are  the  direction  cosines 

*  =  <Xo  -  *>  [ZQ2  +  (x0  -  X)2  +  (yo  -  y)2]-1/2 

«» =  (y0  - y)  fzo  +  (x0  “  x)Z  +  ^o  -  y>2f 1/2 

n  =  sin  <p'  =  zQ[z2  +  (x0  —  x)2  +  (y0-y)2]-1^2 

The  equation  represents  the  sum  of  contributions  falling  on  an  incremental  portion  of  the  ground 
plane  having  direction  cosines  in  the  range  I  to  t  +  df  and  m  to  m  +  dm.  The  Fourier  inversion 
of  Eq.  (3)  is,  of  course, 

OO 

f  (x,  y,  o)  =  ^ m)  exp  [j2?r(fx  +  my)]  d t  dm 

-oo 

ort 

f(x,  y)  m)  (4) 

Imagine  now  that  at  the  target  the  aperture  field  distribution  £  is  actually  a  reflection 
caused  by  a  plane  wave  incident  along  the  ray  path  from  the  direction  shown  in  the  figure,  and 
that  this  wave  arrives  on  the  target  at  a  power  level  of  one  watt  per  square  meter  of  the  incident 
wavefront.  is  then  seen  to  be  the  radar  echo  arriving  at  the  receiving  location.  We  want  to 
determine  the  number  of  watts  per  steradian  of  time-averaged  echo  power  per  square  meter  of 
target  when  the  power  level  of  the  incident  signal  is  one  watt  per  square  meter  of  wavefront. 

This  can  be  obtained  by  taking  one-half  the  squared  magnitude  of  the  complex  envelope  "W]  suit¬ 
ably  corrected  for  the  number  of  steradians  subtended  by  each  incremental  area  df  x  dm.  This 
ratio  is 

22-  l>^f,m)|2  (1  -,’2-m2)l/2  . 

o 

If  the  field  £  on  the  aperture  is  randomly  time-varying,  then  will  likewise  vary  ran¬ 
domly  with  time  and  so  will  the  power  ratios  just  defined.  We  define  the  quantity  to  be  4* 
times  this  average  ratio: 

<70  =  ^-<|T^(f,m)|2>  (1  _i2_m2)l/Z  ,  (5) 

where  <>  denotes  the  infinite  time  average  and  ZQ  is  the  characteristic  impedance  of  free  space; 
<ro  is  called  the  radar  angular  power  spectrum  of  the  target  surface,  and  its  dimensions  are 
square  meters  of  radar  cross  section  per  square  meter  of  target  surface.  The  factor  of  4rr  en¬ 
ters  simply  by  an  arbitrary  convention  whereby  all  radar  cross  sections  are  defined  as  a  ratio 
of  powers  such  as  Eq.  (5)  multiplied  by  4 rr.  As  it  stands,  <r  is  a  function  of  <p,  <p\  6,  /,  and  m. 


t  Note  that  Eqt.  (3)  and  (4)  constitute  a  Fourier  transform  pair  all  right,  but  in  terms  of  direction  cosines  t  and 
m,  and  not  in  either  angles  such  as  $'  or  in  distance  on  the  ground  plane  (x^  and  yQ).  Only  near  perpendicularity 
(x^  as  x  or  yQ  m  y)  does  one  have  ground  distances  as  an  argument. 
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If  one  rules  out  an  anisotropic  target  material,  then  the  dependence  is  only  on  <p,  <p',  and  6.  If, 
further,  a  monostatic  rather  than  bistatic  radar  situation  is  assumed  (that  is,  the  transmitter 
and  receiver  are  collocated)  then  only  one  variable,  <p  =  <p',  is  left.  In  most  of  this  report  this 
assumption  will  be  made  and  we  shall  be  dealing  with  the  radar  angular  spectrum  <j Q{<p). 


(d)  Powsr  vs  range.  (e)  Power  vs  frequency. 

Fig.  3.  The  manner  In  which  the  radar  angular  spectrum  a0(tf)  influences  the  radar  echo  returned 
from  a  uniform  spherical  hard  ta.  »t. 

An  example  of  a  radar  angular  spectrum  is  sketched  in  Fig.  3(a).  Perhaps  more  familiar 
examples  are  the  Lambert  Law  ( aQ  ~  cos2  <p)  and  Lommel-Seeliger  or  Euler  Law  (cr^  ~  cos  <p). 

That  aQ(<p)  expresses  certain  statistics  of  the  echo  wavefront  irregularities  is  now  clear 
from  the  following  reasoning.  From  Eq.  (5),  we  see  that  <?0W)  is  proportional  to  the  time- 
averaged  squared  magnitude  of  the  received  electric  field,  but  this  electric  field  in  turn  is  the 
Fourier  transform  of  the  aperture  distribution  [Eq.  (3)  J.  Thus,  since  the  squared  magnitude  of 
the  function  in  one  domain  is  the  Fourier  transform  of  the  correlation  function  of  the  transformed 
function  in  the  other  domain,  it  follows  that  the  power  observed  as  a  function  of  direction  cosines 
is  proportional  to  the  Fourier  transform  of  the  spatial  correlation  function  of  the  electric  field 
across  the  aperture.  Equations  (3)  and  (4)  are  for  non-time-varying  fields.  As  shown  in  Appen¬ 
dix  A  [in  particular,  Eq.  (A-9)],  if  we  introduce  the  assumption  that  £  is  time-varying,  then  the 
time-averaged  received  power  as  a  function  of  direction  cosines  is  proportional  to  the  inverse 
transform  of  the  spatial  correlation  function  where  the  average  is  over  time,  that  is,t 

cr  =  <  W*,  m>  I  2>  =*  Const  x  $  C  f  <£(u,v)f*(u  +  x,  v  +  y)>  du  dv  .  (6) 

t  A s  discussed  in  Appendix  A,  this  relation  depends  on  the  reasonable  assumption  of  space  stotionoritv  of  the 
aperture  distribution  £  .  That  is,  the  power  level  may  vary  across  the  aperture,  but  apart  from  this,  the  cross 
correlation  of  the  random  time  variations  occurring  at  two  separate  points  on  the  aperture  depends  only  on  the 
separation  of  the  points  and  not  on  their  location. 
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.So  in  measuring  a  ^  q>),  we  are  really  measuring  the  spatial  correlation  function  of  the  wavefront 
irregularities  at  the  target;  this  is  as  good  a  statistical  description  of  the  irregularities  as  has 
been  considered  to  date.t 

By  similar  reasoning  we  can  show  an  equally  important  fact  -  that  the  spatial  correlation 
function  of  the  signal  received  on  the  ground  (time  average  of  the  product  of  two  receiver  outputs) 
is  the  Fourier  transform  of  the  average  power  as  a  function  of  the  coordinates  x  and  y  at  the 
target.  This  fact  is  of  obvious  relevance  to  radar  interferometry,  and  we  shall  Investigate  its 
consequences  in  Sec.  U-D.  The  specific  relationship  [Eq.  (A-10)  of  Appendix  A]  goes  as  follows: 

<|  f(x,  y)  |2>  <=<V(p,  q)>/*(p  +  f,  q  +  m)  dp  dq>  a  Z  (Z,  0)  ,  (7) 

where  we  define  the  (complex)  number  Z  (7,  At)  as  the  cross  correlation  obtained  by  time¬ 
averaging  the  product  (taken  with  a  relative  delay  At)  of  the  two  voltages  received  at  points 
spaced  by  a  vector  separation  7  on  the  ground  plane.  The  station  separation  7  is  characterized 
by  the  direction  cosine  differences  l  and  m. 

B.  Power  Reflected  From  the  Target  as  a  Whole 

Radar  astronomy  targets  are  of  course  very  large  bodies,  usually  spherical.  The  received 
echo  is  made  up  of  the  superposition  of  individual  echoes  from  a  very  large  number  of  individual 
regions,  each  of  which  can  be  regarded  as  approximately  planar,  so  that  the  development  of  the 
preceding  section  is  valid  for  each  region.  Moreover,  if  the  echoes  from  different  regions 
fluctuate  independently,  superposition  of  power  (rather  than  voltages)  is  permissible. 

The  received  power  (per  steradian)  per  square  meter  of  surface  at  a  given  point  on  the  sur¬ 
face  will  be  given,  of  course,  by  the  radar  angular  power  spectrum  <r  (<p)/ 4ir,  where  <p  is  now 
the  tilt  of  the  radar  line  of  sight  from  perpendicularity  with  the  local  horizontal.  (Transmitter 
and  receiver  are  assumed  collocated  and  one  watt  per  square  meter  is  assumed  incident  on  the 
target.)  Figure  3(a)  shows  a  hypothetical  radar  angular  spectrum  plotted  in  polar  coordinates. 

An  ocean  area  on  the  target  could  be  expected  to  give  a  large  narrow  polar  plot  of  &0(<p)  concen¬ 
trated  near  small  values  of  ip,  whereas  the  <r  (?)  function  of  rough  terrain  would  be  lower  in 
magnitude  and  extend  out  to  greater  values  of  <p. 

There  have  been  two  basic  approaches  to  radar  study  of  lunar  and  planetary  surfaces.  The 
first  is  to  assume  that  the  body  is  uniform;  that  is,  on  the  target  surface  <JQW)  is  everywhere 
the  same.  Attention  then  centers  on  finding  out  what  this  vQ(<P)  function  is,  and  then  inferring 
average  terrain  characteristics  by  using  the  methods  discussed  in  Ref.  3,  based  on  Eq.  (6).  The 
second  approach  is  somehow  to  isolate  different  regions  of  the  target  surface,  presumably  hav¬ 
ing  different  aQ(<p)  functions,  and  to  study  the  power  reflected  from  each  region.  The  methods 
of  Ref.  3  then  re-enter  the  picture,  but  now  as  a  tool  in  determining  average  features  of  just  the 
local  region.  (We  shall  see  shortly  that  it  is  not  i.ecessary  to  perform  the  Isolation  by  using 
narrow  antenna  beams;  other  methods  are  available.) 


t  if  the  spatial  wavefront  undulations  at  the  target  are  a  gaussian  random  process,  then  statistical  information  of 
a  higher  order  can  be  inferred  from  first-order  and  second-order  data,  and  these  data  thus  tell  the  whole  story. 

If  non-gaussian  surface  variations  are  suspected,  it  would  be  important  to  consider  some  form  of  experiment  to 
develop  higher-order  statistical  information,  since  in  non-gaussion  cases,  first-order  statistics  ond  second-order 
statistics  (such  as  oq($))  do  not  specify  the  higher-order  statistics. 
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We  now  introduce  the  assumption  that  the  target  is  spherical  and  see  what  can  be  said  about 
the  total  received  echo  power  from  such  a  target.  If  the  spherical  surface  is  uniform  In  the  sense 
that  all  parts  of  the  surface  have  the  same  angular  power  spectrum  <rQ(<p),  say  that  shown  in  polar 
form  in  Fig.  3(a),  then  a  head-on  view  of  the  body  would  appear  as  shown  in  Fig.  3(b),  where  the 
darkness  of  the  shading  is  proportional  to  reflected  power  per  square  meter  of  incident  wavefront. 
This  function  P(x,  y),  representing  echo  power  density  across  a  hypothetical  transverse  plane  in 
front  of  any  target,  is  called  its  radar  brightness  distribution.  The  circular  region  on  this  plane 
projected  by  the  spherical  target  is  called  the  target's  apparent  disk. 

At  a  point  on  the  target  where  the  normal  to  the  spherical  surface  makes  an  angle  <p  with  the 
line  of  sight  from  transmitter  and  receiver  to  the  target,  one  square  meter  of  the  apparent  disk 
will  return  toward  the  receiver  l<rQ(<p)  sec  <s]/4»  watts  per  steradian  for  every  watt  per  square 
meter  of  wavefront  incident.  The  sec  <p  factor  is  necessary  because  we  have  defined  <ro(?)  as 
cross  section  per  unit  area  of  surface  rather  than  per  unit  area  of  wavefront;  there  will  be  sec  <p 
square  meters  of  surface  for  every  square  meter  of  incident  wavefront. 

Figure  3(c)  shows  the  way  in  which  the  polar  plot  of  <7Q(<p)  in  Fig.  3(a)  expresses  the  intensity 
of  the  reflection  at  various  points  on  the  surface  like  point  P.  Note  carefully  that  this  <XQ(q»)  tells 
the  story  only  for  reflection  back  along  the  direction  of  the  incident  radiation;  vQ(^)  does  not  say 
anything  about  the  power  reflected  off  in  other  directions.  (For  this  we  would  need  the  complete 
angular  power  spectrum  oQ(<p,  <p',  9).] 

The  brightness  distribution  P(x,  y)  for  a  uniform  sphere  is  obtained  trivially  from  the  fact 
that  there  are  a0(<p)  sec  ip  square  meters  of  cross  section  per  square  meter  of  wavefront.  Thus, 
P(x,  y)  dfl  is  defined  as  the  number  of  watts  per  square  meter  received  on  the  ground  in  a  cone  of 
unit  solid  angle  da  with  its  vertex  on  the  ground;  all  this  assuming  still  that  the  target  is  illumi¬ 
nated  with  a  plane  wave  of  power  level  one  watt  per  square  meter.  Since  a  square  meter  of  tar¬ 
get  subtends  the  same  solid  angle  at  the  ground  as  vice  versa,  we  have 


■  .  8  ,  x2  +  y2  <  a2  ,  (8) 

/  2  2  2 
va  - x  - y 

where  a  is  the  target  radius. 

The  total  effective  cross  section  <7  of  the  target  can  be  obtained  by  adding  up  all  the  <7Q(<p) 
contributions  over  the  entire  surface  illuminated  by  the  incident  signal.  We  have  for  <7  in 
square  meters 

r a  2 

<7  =  \  2*  RdR  (7  (^>)  sec  ip  =  2ira  \  cr  (^>)  sin  <p  d <p  (9) 

J0  0  J0  ° 

since  R  =  a  sin  < p.  (R  is  the  radius  of  an  annular  strip  of  integration  on  the  apparent  disk  of  the 
body.)  That  is,  if  there  is  an  incident  transmitted  wave  at  a  power  level  of  one  watt  per  square 
meter  of  wavefront,  then  <7/4*  watts  per  steradian  will  be  reflected  back  toward  the  transmitter. 

It  is  interesting  to  rewrite  this  expression  for  total  cross  section  a  in  terms  of  two  param¬ 
eters  of  the  surface  itself,  called  the  reflectivity  p,  and  the  gain  or  directivity  g.  If  the  incident 
transmitted  signal  reaches  the  target  at  a  power  level  of  one  watt  per  square  meter,  then  ob¬ 
viously  a  total  of  ira2  watts  is  incident  on  the  target.  The  total  power  reflected  from  the  target  in 
all  directions  is  [by  steps  similar  to  those  leading  to  Eq.  (9)] 


P(x,  y)  =  a Q(<p)  sec  <p  =  <7Q 
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Cu>/ > el  ect  w.;iy  p  is  defined  as  the  ratio  01  toti  l  reflected  power  to  total  incident  power,  and  is 
fjr  j  from  Eq.  (10), 


i  r,/2  r'/2  .',2' 

=  5^  \  d(p  \  dp1  \  7 „(t»,  v',  8)  siny>  sinyi1 
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if  nttwjty  j  U  defined  in  just  the  same  way  as  is  the  more  familiar  antenna  gain  or  di- 

■  :  i  •')  v.  nf  sif  ly  r  /'it  •  itio: 

T  otal  reflected  power  if  rejection  were  as  strong 

_  _  _ in  all  dir  ections  as  it  iajoward  the  transmitter) 

8  ~  | Actual  total  reflected  power) 


.  Eg.  M>  .  _  *  J9r/2  g0<»> 

Eq.  (Iff)  r-r/2  .  r*/2  .  p2»  ,  ,  .  ,  ...  *2 

Jo  d'p  Jo  d*  J o  %(»>.  V  •  6)  sinyi  sinv>'  de 
We  see  from  Eqa.  (9),  (11)..  and  (12)  that  the  over-all  cross  section  a  can  be  expressed  as 
a  =  xa2  pg  , 

i.e.,  as  the  product  of  three  factors,  ta2  (the  arei  of  the  apparent  disk  at  the  body),  the  quantity 
p,  which  depends  only  on  the  dielectric  properties  of  the  body,  and  g,  which  depends  on  the 
terrain  geometry.  Measurement  of  a  and  its  subsequent  separation  into  these  three  factors  is 
an  important  radar  astronomical  problem  since  knowledge  of  g  gives  considerable  knowledge 
about  terrain  shape,  and  knowledge  of  p  reveals  much  about  the  material  of  which  the  terrain  is 
composed.  However,  as  Eqs.  (11)  and  (12)  show,  joth  p  and  g  depend  on  eQ(y>,  <p',  6),  and  the 
simpler  does  not  give  enough  information  to  compute  p  and  g  separately.  As  Pettengill 
has  pointed  out,^  in  order  to  separa  te  p  and  g,  it  is  necessary  to  make  bistatic  measurements 
over  a  variety  of  angles  ip,  <p\  and  9,  w'.iich  is  not  possible  on  radar  astronomy  targets  from 
the  earth  alone,  but  instead  requires  same  sort  of  radar  space  probe. 


C.  Time  Delay  (Range) 

We  have  just  looked  at  one  way  of  spitting  up  the  echo  power,  saying  that  so  much  of  the 
total  radar  cross  section  may  be  attributed  to  the  projected  area,  so  much  to  reflectivity  and  so 
much  to  directivity.  There  are  other  ways  of  divid  ng  up  the  effect  of  the  target  on  the  echo  pro¬ 
duced,  for  example  according  to  the  time  delay  or  frequency  shift,  and  we  shall  find  these  to  be 
much  more  important.  We  are  particularly  interested  in  any  method  that  will  allow  resolution 
of  portions  of  the  echo  coming  from  different  places  on  the  target. 

Different  portions  of  a  deep  target  will  return  energy  at  different  delays.  Loci  of  all  points 
having  constant  delay  t  will  be  concuntri;  spheres  of  radius  ct/2  about  the  transmitter-receiver 
location,  or  if  transmit  tin-  and  r -.reiver  are  not  collocated,  these  loci  will  be  ellipsoids  of  revo¬ 
lution  with  transmitter  and  receiver  locations  as  foci,  and  a  semimajor  axis  equal  to  cr/2.  If 
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the  target  is  remote  enough  these  concentric  spheres  or  ellipsoids  can  be  approximated  by  par¬ 
allel  planes,  and  we  shall  assume  that  this  is  the  case.  Thus,  at  great  distances  for  a  body  of 
any  shape,  an  observer  standing  at  a  point  perpendicular  to  the  line  of  sight  from  radar  to  target 
center  will  see  a  scale  of  range  reading  linearly  along  the  line  of  sight  [the  solid  lines  in  Fig.  4(b)). 
If  the  target  is  a  sphere,  these  loci  of  constant  range  are  circles  as  one  views  the  target  head-on 
from  the  radar  [the  solid  lines  of  Fig.  4(a)). 


if 


(a) 


Fig.  4.  Contours  of  equal  range  (solid  lines)  and  equal 
doppler  offset  (dotted  lines)  fora  spherical  target.  O 
is  the  rotation  vector  defining  the  polar  axis  of  the 
body.  A  pair  of  regions  (heavily  shaded)  can  be  asso¬ 
ciated  with  the  echo  returning  wilha  given  delay  t  and 
doppler  offset  f. 


Clearly,  one  can  isolate  returns  from  regions  lying  at  or  near  one  of  these  loci  of  constant 
range  by  transmitting  narrow  pulses  and  observing  the  received  return  after  an  appropriately 
adjusted  delay. 

Going  back  to  Fig.  3,  the  sketches  (c)  and  (d)  indicate  how  the  average  power  response  to 
narrow  pulses  would  look  if  the  target  were  uniform  with  a  radar  angular  spectrum  aQ(<p)  of  the 
form  shown  in  (a).  In  Appendix  B  it  is  shown  that  the  distribution  of  power  (actually  cross  sec¬ 
tion)  as  a  function  of  time  delay  t  for  a  uniform  spherical  target  (expressed  in  square  meters) 
is  actually 

<t(t)  =  xacc^cos'1  |^)  ,  (13) 

where  c  is  propagation  velocity  and  a  is  the  target  radius. 

D.  Frequency  and  Phase 

It  has  been  mentioned  several  times  that,  because  of  the  inhomogeneous  nature  of  the  target 
terrain,  there  will  be  irregularly  different  amplitude,  phase,  and  polarization  direction  of  the 
electromagnetic  signal  at  different  places  on  the  wavefront  of  the  reflected  signal,  even  if  these 
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quantities  are  constant  across  the  transmitted  wavefront.  As  long  as  the  transmitter,  receiver, 
and  target  are  rigidly  fixed  relative  to  one  another,  there  being  no  relative  rotation  of  the  target, 
then  each  transmitted  Fourier  component  sinusoid  arrives  at  the  receiver  as  a  sinusoid  with  a 
different  but  constant  amplitude,  phase  and  vector  direction.  However,  at  different  receiver 
points  these  quantities  will  be  randomly  different  constant  values,  if  the  target  is  rough. 

If  now  the  target  rotates,  these  spatial  irregularities  in  the  wavefront  will  pass  sideways 
over  the  receiving  point  at  a  certain  speed  in  a  certain  direction  and  will  give  rise  to  fluctuations 
with  time  of  the  amplitude,  phase,  and  polarization  direction  of  the  sinusoidal  electromagnetic 
vector.?  This  means  that  what  started  out  from  the  transmitter  as  a  sinusoid  will  appear  at  the 
receiving  antenna  output  to  be  a  narrow-band  noise  of  some  type.  The  faster  the  rotation  of  the 
body,  the  wider  the  bandwidth  of  this  narrow-band  noise. 

There  is  a  second  way  of  putting  this:  since  some  portions  of  the  rotating  and  presumably 
rough  target  are  moving  toward  the  observer  and  some  away  from  him  there  will  be  a  spread  or 
smearing  of  the  frequency  spectrum  of  any  sinusoid  incident  on  the  target. 

There  are,  then,  these  two  ways  of  looking  at  the  behavior  of  the  phase  and  frequency  [the 
trigonometric  argument  of  Eq.  (1)]  in  the  received  signal.  Interferometer  techniques,  those  in¬ 
volving  more  than  one  receiving  point,  are  more  conveniently  analyzed  by  considering  the  motion 
sideways  across  the  line  of  sight  of  the  irregularities  in  amplitude,  phase,  and  direction  of 
polarization.  What  happens  at  a  single  radar  receiver  can  be  treated  more  readily  by  looking  at 
the  velocity  of  motion  along  the  line  of  sight  of  various  parts  of  the  target  producing  echoes. 

These  are  just  two  different  ways  of  looking  at  the  same  set  of  phenomena. 

Let  us  look  first  at  the  velocity  of  radial  motion  of  various  points  on  the  target.  It  can  be 
shown  fairly  simply?  that  on  a  distant  rotating  body  the  locus  of  all  points  that  have  equal  com¬ 
ponents  of  velocity  along  the  line  of  sight  to  the  target  center  is  the  line  formed  by  the  inter¬ 
section  of  the  surface  of  the  body  with  a  plane  parallel  to  the  plane  of  line  of  sight  and  polar  axis. 
Moreover,  the  separation  between  this  plane  and  the  plane  of  line  of  sight  and  axis  is  proportional 
to  the  velocity.  This  means  that  a  head-on  view  [such  as  that  in  Fig.  4(a)]  of  a  body  of  any  shape 
would  show  loci  of  target  points  producing  equal  doppler  shift  in  an  incident  sinusoid  as  straight 
lines  running  parallel  to  the  line  formed  by  the  polar  axis.  These  lines  are  spaced  proportion¬ 
ately  with  doppler  shift. 

Clearly,  one  can  isolate  returns  from  regions  lying  on  or  near  one  of  these  loci  of  constant 
doppler  offset  by  transmitting  a  narrow-spectrum  sinewave  and  observing  the  received  return  at 
an  appropriately  adjusted  frequency  offset.  These  loci  are  shown  dotted  in  Figs.  4(a)  and  (b). 

Looking  at  the  target  from  a  distant  point  on  a  line  perpendicular  to  the  radar  line  of  sight 
and  in  the  plane  containing  line  of  sight  and  polar  axis  [Fig.  4(b)],  these  same  loci  of  equal  ve¬ 
locity  toward  the  radar  are  still  straight  lines  spaced  equally  with  doppler  shift.  It  is  now  clear 
why  we  have  investigated  what  the  loci  of  constant  range  and  constant  frequency  look  like  from 
this  particular  distant  point  at  right  angles  to  the  radar  line  of  sight;  for  this  situation  only, 
range  and  doppler  frequency  both  read  linearly  along  two  perpendicular  axes,  irrespective  of  the 


t  The  sinusoid  is  generated  by  motion  ot  the  propagation  velocity  of  the  electremognetic  field  of  wavelength  X 
along  the  line  of  sight.  This  motion  should  not  be  confused  with  the  (much  slower)  sideways  motion  of  the  irregu¬ 
larities  we  hove  been  discussing. 

t  For  example  from  Eq.  (B-6)  of  Appendix  B. 
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Fig.  5.  Method  of  deducing  direction  of  motion  on  the  ground  of  corrugations  in  received  wavefronts 
(a)  hypothetical  cate  of  transmitter  motion  locked  to  target  rotation;  (b)  actual  case  of  transmitter  and 
receiver  stationary  and  approximately  collocated;  (c)  view  of  ground  plane  showing  diffraction  pattern. 
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shape  of  the  body.  This  distortionless  mapping  property6  (of  points  on  the  body  onto  points  in 
the  range -vs -doppler  plane)  will  prove  to  be  a  convenience  later  on. 

If  one  looks  head-on  at  the  target,  as  in  Fig.  4(a),  contours  of  equal  doppler  shift  intersect 
contours  of  equal  range  (which  for  a  sphere  are  circles,  as  described  in  the  preceding  section) 
in  at  most  two  points  as  shown  by  the  shaded  regions.  Thus  we  can  isolate  and  study  such  a  pair 
of  points  on  the  target  surface  simply  by  isolating  that  part  of  the  echo  having  a  given  time  delay 
and  simultaneously  a  given  frequency  shift. 

For  a  uniform  spherical  target,  the  distribution  of  power  as  a  function  of  f,  the  frequency 
offset  from  a  "zero  doppler"  represented  by  the  radial  velocity  of  the  target  center,  is  shown 
in  Appendix  B  to  be 


<r(f)  = 


2f  fi  cos 
o 


(Real  values  only) 


<To(cos'1  fa* 


1  -  (%)z  - 

(  fc  \  2  i/2 

\2f  aO  cos  a / 

0 

dr 


(14) 


(expressed  in  square  meters),  where  f  is  the  carrier  frequency,  Q  is  the  rotation  speed  ii 
radians  per  second  and  a  is  the  tilt  of  the  polar  axis  from  perpendicularity  with  the  line  of  sight. 
Figure  3(e)  illustrates  this  quantity  <r(f)  for  the  particular  <rQ(<p)  sketched  in  Fig.  3(a). 

Now  let  us  look  into  the  other  method  of  visualizing  the  behavior  of  the  phase  and  rate  of 
change  of  phase  of  the  received  signal.  We  wish  to  examine  the  motion  perpendicular  to  the  line 
of  sight  of  the  corrugations  in  the  wavefront  reaching  the  radar  receiver  after  reflection  from  a 
distant  target,  and  we  especially  want  to  know  the  degree  of  correlation  of  the  signals  .  eceived 
at  two  points  spaced  by  a  (vector)  separation  J  as  a  function  of  relative  delay  At  of  /ne  station 
behind  the  other  one  <f  distant.  This  is  expressed  by  the  quantity  Z  introduced  er>lier  in  con¬ 
nection  with  Eq.  (7).  / 

/ 

The  method  we  shall  use  to  deduce  the  motion  of  the  wavefront  irregularities  at  the  receiver 
is  depicted  in  Figs.  5(a)  and  (b).  We  stall  consider  two  situations,  ar.  ■'•'artificial"  but  easily 
analyzed  caae  [Fig.  5(a)],  In  which  th .•  transmitter  rotates  with  the  target.,  /in d  the  "actual"  real- 
life  situation  [Fig.  5(b)],  in  which  oily  the  target  rotates  and  the  transmit  er  and  receiver  are 
cloee  together  and  are  fixed  relative  to  the  target  center.  / 

Three  arbitrary  points  forming  a  triangle  are  chosen  on  the  target  Cody  itself.  Consider 
points  1  and  2  of  Fig.  5(a)  or  (b).  Line  AB  is  the  locus  of  all  points  */ a  plane  at  the  receiver 
(perpendicular  to  the  line  of  sight)  for  which  the  phase  angle  of  the  composite  sinusoidal  signal 
reflected  by  points  1  and  2  is  the  same  value  as  at  the  receiver.  That  is,  it  is  the  locus  of  all 
points  having  the  same  difference  in  time-of-flight  from  trar  am liter  to  receiver  via  the  two 
pointe  1  and  2.  Similar  ly,  the  line  CD  is  the  locus  of  all  rvoints  satisfying  a  similar  condition 
regarding  propagation  via  points  1  and  3.  Considering  »ow  that  points  1,  2,  and  3  are  actually 
in  motion  in  different  directions  due  to  target  rotat'  ^h,  it  is  possible  to  find  the  direction  and 
speed  of  motion  of  the  wavefront  irregularities  ^ifthe  receiver  by  finding  the  direction  and  speed 
of  motion  of  point  Q»  the  intersection  of  AF  and  CD. 

Before  getting  to  the  desired  cose  transmitter  and  receiver  collocated  [Fig.  5(b)],  let  us 
analyze  the  artificial  situation  depicte  i  in  Fig.  5(a).  Here  the  receiver  is  fixed,  but  the  trans¬ 
mitter  is  imagined  to  swing  around  as  the  target  rotates,  always  being  at  the  same  distance  from 


the  target  and  with  the  same  orientation  relative  to  it *  Thus,  as  the  target  rotates,  the  reflected 
wavefront  at  the  target  has  its  phase,  amplitude  and  polarization  corrugations  "frozen";  that  is, 
referring  to  Fig,  5(a),  the  times  of  flight  from  transmitter  to  points  such  as  1,  2,  and  3  do  not 
change.  A  moment's  reflection  convinces  one  that  the  direction  of  motion  of  the  corrugations  as 
observed  on  the  ground  must  be  in  the  direction  in  whicli  the  line  OQ  would  swing  if  it,  too,  were 
frozen  to  the  target  at  point  O. 

The  important  thing  about  this  artificial  example  is  that  if  we  imagine  two  points  Q  and  Q1, 
equidistant  from  the  target  and  mutually  separated  by  a  distance  QQ'  lying  parallel  to  the  direc¬ 
tion  of  pattern  motion,  then  the  following  effect  takes  place.  The  pattern  of  irregularities  that 
sweeps  over  the  receiver  point  Q  (and  thus  generates  effectively  a  noise-like  fluctuation  in  the 
receiver  output  at  Q)  is  identical  to  the  pattern  at  the  point  Q1  at  a  later  time  equal  to  the  time 
taken  for  the  pattern  to  move  the  distance  QQ1.  The  velocity  of  pattern  motion  in  the  neighbor¬ 
hood  of  Q  is  clearly 

v  =  OR  cos  a 

where  R  is  the  range  and  a  is  the  polar  angle,  so  that  the  time  delay  is 

Ar  =  . 

o  SIR  cos  o/ 

For  radar  astronomy  targets,  all  the  target  points  are  essentially  at  the  same  range  R.  Fig¬ 
ure  5(c)  shows  a  plan  view  of  the  ground  plane  with  the  vector  v  indicated  by  the  arrow,  and  a 
snapshot  of  the  variations  of  amplitude  and  phase  across  the  wavefront  schematized  by  the  curved 
contours.  Under  certain  conditions  it  may  be  assumed  that, on  the  average,  these  contours  are 
circles;  that  is,  the  pattern  changes  equally  rapidly  with  distance  in  all  directions  on  the  ground 
plane.  Among  these  conditions  is  the  case  of  a  spherical  target  with  uniform  surface  properties, 
i.e.,  uniform  cr Q(<p)- 

Imagine  now  that  we  have  two  receivers  R^  and  R^  spaced  by  the  vector  distance  7,  making 
an  angle  6  with  the  pattern  velocity  v,  and  that  we  wish  to  compute 

Z(  d,Ar)  =  <V(t)\^*(t  +  At)>  , 

the  cross  correlation  of  the  electric  field and  observed  at  the  two  receivers,  as  a  function 
of  the  separation  7  and  a  lag  Ar  of  the  first  signal  behind  the  second.  The  quantity  Z  is  complex 
in  general;  that  is,  it  can  be  viewed  as  embodying  the  phase  and  amplitude  of  a  sinusoid.  (Just 
how  this  sinusoid  can  actually  be  generated  in  an  interferometer  experiment  will  be  discussed  in 
Sec.  IV-F.)  We  will  now  describe  heuristically  how  a  Z{ 7,  At),  such  as  is  depicted  in  Fig.  6(a) 
expresses  certain  target  properties.  (This  is  an  extension  of  Eq.  (7)  in  which  Z (7,  0)  was  re¬ 
lated  to  the  power  across  the  target  aperture.) 

For  7  =  0,  this  function  Z  is  the  autocorrelation  function  of  a  single  receiver  output,  and 
we  already  know  the  Fourier  transform  of  this;  it  is  just  the  power  spectrum  of  the  fluctuation 
<r(f)  as  expressed  in  Eq.  (14)  and  depicted  in  Fig.  3(e),  except  that  in  this  case  this  power  spec¬ 
trum  is  half  as  wide,  because  only  the  target -to-receiver  distance  is  changing,  the  target-to- 
transmitter  distance  being  constant.  Denoting  the  transform  of  cr(f)  by  «(At),  we  have! 


t  Equivalently,  the  transmitter  could  be  fixed  and  the  receiver  move  synchronously  with  target  rotation, 
t  A(At)  is  a  quantity  to  be  discussed  further  in  Sec.  IV-E. 
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Z  (0,  At)  =  «(At/2)  . 

It  la  almost  equally  simple  to  determine  the  zero-delay  correlation  Z(S,  0),  still  sticking  to 
our  "artificial"  example  of  the  synchronously  rotating  transmitter.  The  function  Z(S,  0)  is  de¬ 
rivable  from  Eq.  (7),  which  states  a  Fourier  transform  relationship  between  ground  space  cor¬ 
relation  function  (with  zero  relative  delay)  and  power  across  a  plane  aperture.  If  we  imagine 
the  aperture  to  be  the  apparent  disk  of  the  target,  then  the  distribution  of  power  across  the  ap¬ 
erture  <|  £  (x,  y)  |  2>  in  Eq.  (7)  is  actually  the  brightness  distribution  P(x,  y),  defined  earlier. 
Thus 


P(x, y)  exp  [j2w(px  +  qy)]  dx  dy  =  (W^(f,  m)  W*(t  +  p,  m  +  q)>  d/dm  *  Z(S,  0) 

Aperture 

That  is,  the  inverse  transform  of  the  brightness  distribution  is  the  ground  correlation  Z  taken 
with  time  shift  At  =  0,  and  with  station  separation  d  given  by  the  two  direction  cosine  arguments 
p  and  q.  With  no  loss  of  generality,  we  can  take  d  parallel  to  one  of  the  axes  of  the  aperture  in 
Fig.  2,  say  the  x  axis  (so  that  q  =  0);  then 

C  exp  [j2irpx]  dx  C  P(x,  y)  dy  =  Z (3,  0) 

•'Aperture  •'Aperture 

In  other  words,  the  right  side  is  given  by  the  single  inverse  transform  (in  x)  of  the  integral  along 
y  of  the  brightness  distribution  P(x,  y).  Now  notice  that  for  a  uniform  spherical  target  the  shape 
of  Z(S,  0)  is  the  same  (except  for  a  scale  factor)  as  Z( 0,  At).  This  has  to  be  so  because  the  sin¬ 
gle  Fourier  transform  with  respect  to  At  of  the  latter  is  just  the  echo  power  spectrum  o(f),  which 
in  turn  is  simply  the  integral  of  the  brightness  distribution  P(x,  y)  along  one  of  the  constant  fre¬ 
quency  strips  in  Fig.  4(a),  i.e.,  the  integral  /Aperture  P(x,y)  dy. 

About  the  only  other  thing  we  shall  need  to  know  about  Z(S,  At)  in  order  to  discuss  inter¬ 
ferometry  experiments  later  is  the  value  of  At  that  will  maximize  Z  for  a  given  station  spac¬ 
ing  <T.  Referring  to  Fig.  5(c),  in  which  the  spacing  7  lies  between  points  Rj  and  R2>  we  see  that 
if  the  ground  pattern  has  the  generally  isotropic  character  assumed  earlier,  maximum  correla¬ 
tion  will  be  achieved  if,  in  performing  the  correlation,  the  first  signal  is  delayed  by  a  At  equal 
to  the  time  taken  for  the  pattern  to  move  the  projected  distance  Rjft^;  that  is, 

*  _  d  cos  5  _  d  cos  6 

Tmax  _  v  "OR  cos  a 

We  shall  make  use  of  such  a  relation  later  in  deducing  the  planetary  rotation  vector  S  from  v, 
the  ground  pattern  velocity.  Note  that  the  more  nearly  d  lies  along  v,  the  larger  the  value  of 
Z  obtained  by  inserting  the  best  delay  ATmax  in  processing  the  two  receiver  outputs;  indeed, 
when  6  =  0,  the  full  correlation  Z(S,  ^Tmax)  =  Z{ 0,  0)  is  obtained  because,  since  the  diffraction 
pattern  is  frozen,the  two  received  voltages  are  Identical  apart  from  a  time  shift. 

With  this  "artificial"  example  now  disposed  of,  most  of  the  groundwork  has  been  laid  for 
discussing  the  "actual"  case.  Figure  5(b)  shows  the  transmitter,  target,  and  ground  plane  for 
this  situation,  and  now  the  transmitter  is  assumed  fixed.  To  find  the  direction  and  speed  of 
motion  of  the  ground  pattern  we  once  again  examine  the  motion  of  point  Q,  the  intersection  of 
equiphase  lines  AB  (due  to  reflection  from  target  points  i  and  2),  and  CD  (due  to  reflection  from 
1  and  3). 
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The  actual  rotation  vector  of  the  target  ia  5.  We  can  express  this  as  the  vector  sum  of  two 
components  ftj  and  fl2  along  the  line  of  sight  and  perpendicular  to  it,  respectively.  We  have 
f)1  =  0  sina,  and  f)2  =  0  cos  a.  We  can  get  the  over-all  behavior  of  lines  AB  and  CD  by  con¬ 
sidering  these  two  rotation  velocities  separately,^  knowing  that  the  true  motion  of  points  1,  2, 
and  3  is  the  vector  sum  of  motions  due  to  two  rotations  (1^  and  n2.  Now  if  we  seek  the  motion 
of  point  Q,  the  intersection  of  AB  and  CD,  it  is  clear  that  component  0^  has  no  effect  on  lines 
AB  and  CD  except  that  they  rotate  at  an  angular  rate  0^,  the  point  of  intersection  Q  remaining 
fixed.  This  is  true  because  the  propagation  path  lengths  are  not  changing.  Therefore  we  need 
only  consider  the  motion  of  point  Q  due  to  fi2^  Since  02  is  an  angular  velocity  that  the  target 
would  have  if  its  polar  axis  were  perpendicular  to  the  line  of  sight  and  it  rotated  with  a  speed 
a2  =  0  cos  a,  we  can  state  that  the  motion  of  point  Q  is  exactly  the  same  in  our  "actual"  ex¬ 
ample  as  it  was  in  the  "artificial”  example  of  Pig.  5(a),  except  that  the  speed  is  twice  what  it 
was  before,  since  both  the  propagation  paths  (up  and  back)  to  each  target  point  are  changing  their 
length  with  time,  so  that  point  Q  must  move  twice  as  fast  toward  Q'  to  equalise  path  length  dif¬ 
ference  via  points  1  and  2  and  via  1  and  3. 

Note  that  the  contribution  to  the  ground  pattern  produced  by  reflection  from  any  point  on  the 
target  has  the  same  velocity  of  motion  as  that  produced  by  a  reflection  from  any  other  point, 
namely  in  direction  QQ',  and  with  twice  the  speed  that  a  radius  vector  OQ  from  the  target  center 
would  swing, 

v  =  2QR  cos  a 

This  means,  among  other  things,  that  there  is  no  hope  of  distinguishing  the  echoes  from  differ¬ 
ent  portions  of  the  target  on  the  basis  of  different  velocities  of  their  contributions  to  the  ground 
pattern.  We  must  go  back  to  resolution  according  to  range  (Sec.  II-B)  or  doppler  (Sec.  II-C)  or 
the  brute  force  method  of  using  sharper  antenna  beams. 

There  is  a  second  difference  between  the  "artificial"  and  "actual"  cases  that  must  be  pointed 
out  besides  the  double  velocity  of  the  ground  diffraction  pattern  motion.  Since  the  propagation 
path  lengths  from  transmitter  to  target  are  constantly  changing  as  the  target  rotates  in  the 
"actual"  case,  the  wavefront  irregularities  of  amplitude,  phase,  and  polarization  direction  at  the 
target  are  not  "frozen"  as  was  true  in  the  "artificial"  case,  but  are  constantly  changing  with  time. 
Thus  the  contours  of  Fig.  5(c)  are  not  only  in  motion,  as  before,  but  are  changing  with  time. 

This  means  that  there  will  be  an  additional  diminution  of  the  correlation  Z  (7,  At)  from  the  maxi¬ 
mum,  over  and  above  the  factors  discussed  for  the  artificial  case.  (These  factors  were  the  fail¬ 
ure  of  the  direction  of  7  to  coincide  with  that  of  v,  and  then,  of  course,  the  nonzero  values  of 
variables  d  and  At.)  In  particular,  even  if  one  inserts  the  At  that  maximizes  Z  for  a  given 
7,  namely  ATmax,  which  is  now  given  by 


a_  _  d  cos  i 
Tmax  “  2(1  It  cos  a 


(15) 


he  does  not  get  the  full  correlation  Z{0,  0 )  when  7  lies  along  v  because  the  two  received  volt¬ 
ages  are  no  longer  exactly  identical  apart  from  a  time  shift. 


tThe  author  It  indebted  to  J.  V.  Event  for  pointing  out  thii  timple  method  of  analyzing  the  motion  of  point  Q  in 
the  "actual"  cate. 

t  True  except  for  the  lingular  cate  of  widely  ipaced  trammitter  and  receiver  pointt  and  the  target  axit  appearing 
almoit  exactly  end-on  from  the  ground  (flj  » Clj). 
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If  the  time  rate  of  this  pattern  change  were  a  sufficiently  large  fraction  of  the  fading  rate 

produced  in  a  single  antenna  due  to  motion  of  the  pattern,  then  we  would  know  that  Eq.  (IS)  would 

not  be  valid.  The  best  value  of  At,  At  would  be  shifted  toward  zero  because,  since  the  pat- 

max  r 

tern  is  decomposing  as  it  moves,  a  higher  degree  of  correlation  would  be  obtained  by  advancing 
At  slightly  to  head  off  the  decomposition.  Thus,  before  we  can  be  sure  Eq.  (15)  is  valid,  we 
must  be  sure  that  the  pattern  decomposition  can  be  neglected  compared  to  pattern  motion. 

We  can  determine  how  big  a  part  pattern  fluctuation  plays  relative  to  pattern  motion  by  in¬ 
quiring  how  much  of  the  signal  fluctuation  at  a  single  receiver  is  due  to  motion  and  how  much  is 
due  to  pattern  changes  with  time.  It  will  now  be  shown  (in  the  next  paragraph)  that  al).  the  fluc¬ 
tuation  rate  observed  at  a  single  receiver  is  due  to  pattern  motion;  it  then  being  a  safe  assump¬ 
tion  that  for  spacings  d  of  the  order  of  several  multiples  of  the  correlation  distance,  pattern 
decomposition  is  small.  Equation  (IS)  is  then  valid  except  for  the  case  where  one  simultaneously 
has  very  small  values  of  the  angle  A  and  large  values  of  d. 

Comparison  of  Figs.  6(a)  and  (b)  shows  that  the  profile  of  Z (0,  At)  for  the  "artificial"  case 
of  a  frozen  pattern  of  irregularities  is  twice  as  wide  as  the  same  profile  for  the  "actual"  case 
about  which  we  are  concerned.  This  is  known  from  considering  the  radial  motion  relative  to  the 
observer  of  the  different  points  on  the  target;  a  given  such  velocity  [Fig.  4(a))  will  produce  one- 
half  the  doppler  shift  when  the  transmitter-to-target  distance  is  not  allowed  to  change  [Fig.  5(a)  — 
the  "artificial"  case]  as  when  the  transmitter  and  receiver  are  collocated  [Fig.  5(b)).  Yet  we  have 
already  seen  from  different  considerations  that  the  sideways  pattern  velocity  across  the  ground 
in  the  "actual"  case  is  exactly  twice  that  of  the  "artificial"  case.  But  the  physical  size  of  the 
pattern  corrugations  is  the  same  in  both  cases,  as  ensured  by  the  fact  that  Z(A,  0)  is  the  same 
for  both  cases.  Thus,  if  the  physical  size  of  the  irregularities  is  the  same,  and  the  velocity  of 
motion  along  the  ground  is  twice  as  large  in  the  "act  1"  case,  then  this  completely  accounts  for 
the  fact  that  the  echo  power  spectrum  is  twice  as  wide  in  the  actual  case.  The  pattern  decompo¬ 
sition  makes  a  negligible  contribution  to  the  fading  rate,  relative  to  the  pattern  motion. 

It  is  instructive  to  compare  these  notions  on  rotating  planets  with  equivalent  ideas  that  have 
been  applied  for  some  years  to  the  ionosphere.  If  one  imagines  a  plane  wave  incident  on  a  plane 
ionosphere  overhead,  then  if  the  ionospheric  irregularities  are  frozen,  the  horizontal  motion  of 
the  ionosphere  (an  ionospheric  wind)  causes  a  motion  along  the  ground  of  the  frozen  Irregularities 
in  the  wavefront  pattern.  In  actuality,  however,  since  the  ionosphere  is  a  soft  target  in  which 
there  are  turbulent  effects,  constant  random  alteration  of  the  wavefront  irregularities  is  present. 
With  the  ionosphere,  the  inability  to  recover  complete  correlation  by  readjusting  r  even  when 
7  coincides  with  v  is  due  to  the  changing  shape  of  the  target;  with  a  planet  it  is  due  to  the  fact 
that  the  transmitter  is  stationary  instead  of  moving  synchronously  with  target  rotation.  Briggs, 

•7 

et  al.,  have  made  a  lengthy  analysis  of  the  use  of  such  interferometric  methods  to  study  ground 
pattern  velocity  and  the  way  in  which  the  moving  ground  pattern  itself  fluctuates  with  time.  Their 
analysis  goes  a  good  deal  further  for  the  ionospheric  case  than  the  notions  we  have  presented 
here  for  the  rotating  sphere  case.  The  aim  here  has  been  simply  to  determine  i-‘Tmax  [Eq.  (15)) 
without  deriving  the  complete  function Z (?,  At). 

E.  Polarisation 

Most  of  the  radio  and  visible  energy  from  cosmic  sources  is  relatively  unpolarized;  that  is, 
most  of  the  energy  in  the  arriving  signal  w(t)  not  only  has  a  randomly  time-varying  amplitude  and 
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phase  \W(t)  a  random  complex  variable],  but  also  a  random  time-varying  direction  of  polariza¬ 
tion.  By  contrast,  in  making  radar  observations  the  character  of  the  polarization  is  controlled 
on  transmission  for  the  dual  purposes  of  conserving  system  efficiency  (by  not  wasting  power  into 
unnecessary  degrees  of  freedom  of  the  signal)  and  of  studying  those  target  properties  that  affect 
the  polarization.  Since  the  target  does  not  in  general  produce  a  complete  randomization  of  the 
incident  wave,  it  is  often  found  that  a  large  fraction  of  the  received  energy  arrives  systematically 
polarized  (see,  for  example,  Refs.  5,  8  and  9  for  certain  lunar  results). 

Such  partly  polarized  radiation  can  be  resolved  uniquely  into  two  parts,  the  "unpolarized" 
(i.e.,  randomly  polarized)  component,  and  the  (systematically)  "polarized"  component.  The  un¬ 
polarized  component  will  have  a  narrow-band  noise  spectrum  (due,  for  example,  to  the  doppler 
spreading  discussed  in  Sec.  Il-D)  and  if  resolved  (for  example  into  two  linearly  polarized  com¬ 
ponents  at  right  angles  to  one  another)  will  exhibit  zero  correlation  between  the  two.  The  polari¬ 
zed  component,  on  the  other  hand,  may  or  may  not  have  a  nonzero  bandwidth.  In  either  case  the 
components  at  right  angles  are  mutually  coherent.  That  is,  the  electric  vector  traces  out  an 
ellipse  (Fig.  7)  which  may  be  continually  fluctuating  in  size  (at  the  bandwidth  rate),  but  maintains 
a  constant  orientation,  ratio  of  axes  and  sense  of  rotation. 


Fig.  7.  Path  traced  by  head  of  a  vector  representing 
received  electric  field  strength  as  a  function  of  time 
(complete  polarization). 


H-H-Wt  I 


The  polarized  component  is  specified  by  its  power  P  ,  orientation  angle  x  and  axial  ratio  r. 
The  sense  of  rotation  is  contained  in  the  sign  of  r;  that  is,  r  >  0  for  left-hand  circular  polariza¬ 
tion  and  r  <  0  for  right-hand  circular.  In  Fig.  7,  the  wave  is  propagating  along  the  z  axis  into 
the  paper;  thus  x  is  clockwise  from  the  x-axis,  and,  as  the  arrow  indicates,  r  is  negative. 

The  unpolarized  part  is  specified  by  a  single  parameter  P  ,  its  power.  The  degree  of  polar¬ 
ization  m  is  simply  Pe/P,  where  P  equals  the  total  power  Pe  +  Pu>  Since  in  radar  experiments 
the  wave  is  transmitted  completely  polarized  (m  =  i),  the  term  depolarization  (1  -  m)  is  perhaps 
more  frequently  encountered  in  radar  usage  than  is  the  degree  of  polarization  m.  A  monochro¬ 
matic  (i.e.,  single-frequency)  signed  is  completely  polarized,  by  definition.  It  should  be  noted 
that  one  sometimes  encounters  other  definitions  of  depolarization  in  which  a  monochromatic 
wave  can  be  said  to  be  "depolarized."  For  example,  if  vertical  polarization  is  transmitted,  the 
horizontally  polarized  received  component  is  sometimes  called  the  "depolarized"  component. 

Such  a  definition  obviously  is  not  equivalent  to  that  employed  here. 
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There  ere  eeverel  different  weye  of  rearranging  the  quantities  P,  m,  x  and  r  into  a  more 
usable  form,  and  the  preferred  rearrangement  depends  on  the  particular  properties  of  the  radia¬ 
tion  that  are  most  interesting.  For  radio  engineering  purposes  the  most  interesting  choice  is 
usually  the  four  Stokes  parameters10 

I  =  P 

Q  =  ml  cos  2(tan_1  r)  cos  2x 

U  =  ml  cos  2(tan_1  r)  sin  2x 

V  =  ml  sin  2(tan_1  r)  (16) 

because  of  the  ease  with  which  the  effect  of  various  forms  of  targets  and  receiving  antennas  may 
be  calculated.  The  parameters  have  the  following  physical  interpretation:  I  is  of  course  the 
total  power,  Q  is  the  excess  of  x  *  0  linearly  polarized  power  over  that  at  right  angles  to  it,  U 
is  the  excess  of  x  =45*  linearly  polarized  power  over  that  at  right  angles  to  it,  and  V  is  the  ex¬ 
cess  of  left-hand  circular  over  right-hand  circular  power.  The  four  quantities  are  often  written 
as  a  column  matrix  or  vector,  and  are  then  referred  to  as  the  Stokes  vector  S .  (Note  that  this 
is  a  vector  in  the  mathematical  sense  and  is  not  a  vector  property  of  the  wave  in  physical  space.) 

The  usefulness  of  the  Stokes  representation  in  deducing  the  effects  of  signals  on  specified 
antennas  is  clear  from  the  above  description  of  the  parameters.  The  effect  of  a  target  on  the  in¬ 
cident  signal  turns  out  to  be  simply  defined  too.  Specifically,  the  Stokes  vector  Sr  of  the  re¬ 
flected  wave  is  related  to  the  Stokes  vector  S.  of  the  incident  wave  by  the  set  of  simple  linear 
equations  given  by 

Sr  =  MSj  ,  (17) 

where  M  is  a  four-by-four  matrix  called  the  Mueller  matrix.* 

The  Mueller  matrix  representation  has  proved  quite  useful  in  studies  of  simple  radar  targets 
such  as  cones,  comer  reflectors,  and  so  forth,11  especially  when  the  target  is  maintaining  a  fixed 
aspect,  thus  returning  a  monochromatic  (completely  polarized)  signal. 

Of  central  importance  in  the  study  of  a  hard  rotating  radar  astronomy  target  is  the  degree  of 
polarization 


Only  if  the  body  is  extremely  simple  in  shape  or  not  rotating  at  all  relative  to  the  line  of  sight 
does  one  expect  m  to  be  unity.  The  more  realistic  situation  involves  complicated  diffractive 
and  multiple  reflection  effects  giving  rise  to  the  wavefront  irregularities  we  have  spoken  of 
earlier.  The  way  in  which  m  drops  from  the  unity  value  it  would  have  for  an  ideally  smooth 
shape  (for  example  a  sphere)  as  the  surface  roughness  increases  appears  not  to  have  been  studied 


t  A  motter  of  terminology  it  worth  mentioning  here.  The  elementt  of  the  Mueller  matrix  (which,  like  thote  of  the 

I  Stokes  vector,  deal  with  powers)  are  tumt  of  products  of  the  elements  of  a  matrix  of  voltooes  that  characterizes 

the  target  alto.  This  two-by-two  matrix  of  voltages  is  colled  the  scattering  matrix  of  the  target  and  It  something 
entirely  different  from  the  taroet  scattering  function  to  be  defined  in  Sec.  IV.  The  reader  may  already  be  familiar 
with  the  scattering  matrix;  it  should  not  be  confused  with  the  scattering  function. 

f 
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in  the  literature.  Indeed  it  is  not  entirely  clear  that  m  should  drop  completely  to  zero  as  sur¬ 
face  roughness  increases.^ 

The  theory  of  depolarization  of  electromagnetic  waves  by  reflection  from  rough  surfaces 
seems  to  have  been  given  very  little  attention, t  although  a  great  deal  of  experimental  data  exists. 

P.  Beckmann  and  his  coworkers  have  made  an  encouraging  start  on  the  theory.  Beckman12  has 
treated  the  case  of  random  surfaces  with  radii  of  curvature  much  larger  than  the  wavelength,  and 
Chytil13  has  examined  a  variety  of  target  shapes  (such  as  disks,  strips,  cylinders,  spheres  and 
prisms)  in  the  region  X  ~  a  to  X  =  a/10  (where  a  is  the  target  size).  Chytil  finds  that  with  linearly 
polarized  radiation  incident,  as  X  -»  a  the  echo  power  becomes  about  equally  divided  between  com¬ 
ponents  parallel  with  and  perpendicular  to  the  incident  polarization  direction  for  all  shapes. 

m.  SOFT  TARGETS 

According  to  the  definition  adopted  earlier  in  this  chapter,  soft  targets  include  those  that  are 
changing  shape,  such  as  turbulent  clouds  of  one  sort  or  another.  The  soft  targets  of  greatest  ra¬ 
dar  interest  are  the  ionized  regions  (sun's  corona,  planetary  ionospheres,  and  the  interplanetary 
medium),  and  this  part  of  the  present  chapter  will  focus  on  the  influence  of  such  regions  on  the 
transmitted  signal. 

On  the  other  hand,  clouds  of  material  particles,  such  as  the  zodiacal  particles,  terrestrial 
raindrops,  etc.,  are  of  only  minor  radar  astronomical  interest  and  it  is  not  necessary  to  say 
very  much  about  them  here.  They  are  usually  so  rarified  that  multiple  reflections  can  be  neg¬ 
lected,  and  a  computation  of  radar  cross  section  of  a  cloud  or  portions  of  it  consists  simply  of 
multiplying  the  cross  section  of  a  single  typical  particle  by  the  number  of  particleB  in  the  region. 

If  the  pSrticles  are  perfectly  conducting  spheres,14  for  example,  the  cross  section. per  particle 
is  that  given  by  Fig.  8,  in  which  it  is  seen  that  when  the  radius  a  is  much  smaller  than  the  wave¬ 
length  X,  the  cross  section  falls  off  as  the  fourth  power  of  a/x  (the  "Rayleigh  scattering  law")  § 
This  a/X  «  1  condition  holds  with  such  soft  targets  as  the  zodiacal  regions  and  terrestrial  rain 
clouds,  and  it  also  holds  in  a  rough  quantitative  way  with  ionospheric  free-electron  scattering.15 
(The  cross  section  per  free  electron  is  obtained  by  substituting  the  classical  electron  diameter 
in  the  Rayleigh  law.)  It  is  interesting  to  observe  in  Fig.  8  that  in  the  spherical  case  as  one  passes 
from  the  a/X  «  1  situation  (relevant  to  some  of  the  soft  radar  astronomy  targets)  over  into  the 
a/X  »  1  case  (to  which  our  discussion  of  hard  targets  has  been  directed),  the  cross  section  ap¬ 
proaches  x&2,  the  projected  area  of  the  sphere,  as  would  be  expected  from  the  discussion  in 
Sec.  II-B.J 

t  In  reflection  from  random  medio  where  there  ore  no  multiple  reflections  and  the  individual  scatteron  have  a 
strong  anisotropy,  the  depolarization  may  not  be  complete.  This  is  suggested  by  a  study  of  the  reflection  of  RF 
signals  from  clouds  of  randomly  oriented  dipoles'^  for  which  it  was  found  that  m  =  1/2. 

t  The  Stokes  vector- Mueller  matrix  formalism  might  hove  some  application  to  hard  radar  astronomical  targets. 
(Parke'' gives  a  very  complete  survey  of  this  formalism.)  This  is  suggested  by  the  following  interesting  properties. 

In  the  case  of  multiple  reflections,  the  M  describing  the  over-all  process  is  the  product  of  the  M's  describing 
the  individual  reflections,  and  when  a  number  of  wavelets  odd  to  give  a  composite  signal,  although  the  over-all 
M  Is  not  the  sum  of  the  individual  M's,  the  summation  can  be  handled  in  terms  of  the  scattering  matrices  of  the 
reflections.  Conceptually  at  least,  this  sort  of  formalism  then  provides  the  possibility  of  theoretically  studying 
the  polarization  properties  (Sr)  of  an  echo  made  up  of  a  sum  of  multiple  reflections. 

§  It  is  not  necessary  that  the  body  be  spherical  for  the  fourth-po  >r  low  to  hold. 

1  Notice  here  that  the  peculiar  method  of  defining  all  radar  cross  sections  as  a  power  ratio  timet  an  arbitrary  4e 
(referred  to  in  connection  with  Eq.  (5)]  does  hove  the  advantage  that  the  cross  section  of  a  perfectly  conducting 
sphere  equals  its  pro) acted  area.  It  it  difficult  to  imagine  a  single  other  advantage  of  ‘hit  curious  convention. 
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Fig  -  8.  Wav* length  dependence  of  the  radar 
cross  wet  ion  of  a  perfectly  conducting  sphere 
of  radius  a. 


As  we  have  just  indicated,  soft  targets  of  radar  astronomy  interest  are  limited  for  the  most 
part  to  ionized  media,  since  clouds  of  material  particles  (such  as  the  zodiacal  particles)  suffer 
the  (a/x)4  dependence  and  are  at  prohibitively  large  ranges;  and  those  effects  in  our  own  atmos¬ 
phere  not  associated  with  its  ionized  condition  are  minor. 

In  discussing  the  reflection  of  electromagnetic  signals  from  ionized  media  it  will  be  neces¬ 
sary  here  only  to  review  one  or  two  points  and  to  emphasize  the  dissimilarities  of  soft  targets 
and  hard  targets. 

A.  Reflected  Power  and  Angular  Power  Spectrum 

The  effective  radar  cross  section  of  a  given  volume  of  plasma  can  vary  over  the  enormous 
range  from  simply  the  sum  of  the  cross  sections  of  individual  free  electrons  (free  electron 
scatter)  to  the  unity  reflectivity  obtained  in  the  case  of  reflection  below  the  critical  frequency 
(in  the  absence  of  collision  losses).  Exactly  where  in  this  range  the  effective  cross  section  of 
a  unit  volume  of  plasma  will  lie  thus  depends  on  many  factors,  principally  the  operating  fre¬ 
quency,  electron  and  ion  densities,  collision  frequency,  magnitude  and  direction  of  the  magnetic 
field  and  the  temperature. 

When  the  target  presented  to  the  probing  radar  signal  is  a  large  and  inhomogeneous  body  of 
plasma,  such  as  the  sun's  corona,  then  matters  become  quite  complicated.  Unlike  the  hard  tar¬ 
gets  discussed  earlier,  the  target  material  is  in  constant  relative  motion  so  that  it  is  impossible 
to  use  most  of  the  neat  properties  one  can  attribute  to  rotating  hard  targets  (such  as  the  location 
of  loci  at  constant  range  or  doppler  offset).* 

However,  under  certain  conditions  it  may  still  be  useful  to  retain  some  of  these  ideas,  par¬ 
ticularly  the  notion  of  radar  angular  spectrum.  For  example,  although  our  own  ionosphere  is  a 
soft  target,  there  are  some  reasonable  constraints  on  shape  (particularly  the  common  assumption 
that  it  has  a  general  horizontal  stratification,  on  which  are  superimposed  statistical  irregularities), 
and  the  notion  of  angular  spectrum  of  reflection  from  the  irregularities  has  proved  quite  meaningful 


t  And  even  If  the  ionized  plasma  target  were  rigid,  its  geometrical  properties  would  be  more  difficult  to  evaluate, 
since  we  are  not  always  dealing  her*  with  reflection  from  objects  but  mors  often  with  penetration  of  the  wav*  to 
point*  where  certain  condition*  on  refractive  index  are  fulfilled  (these  conditions  being  dependent  on  the  direc¬ 
tion  of  arrival).  Thus  it  is  even  more  difficult  to  associate  the  angular  power  spectrum  with  detailed  structure  of 
the  medium  than  It  is  with  a  hard  target.  (This  association  for  the  hard  target  case  is  discussed  in  Ref.3.) 
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(see  the  last  paragraph  of  Sec.  II -D).  A  first -order  approximation  to  the  true,  highly  irregular 
state  of  affairs  in  the  sun's  corona  has  been  to  assume  that  it  is  spherically  symmetric  with  a 
general  radial  decrease  of  electron  density,  on  which,  as  with  the  ionosphere,  local  spatial  irreg¬ 
ularities  are  superimposed,  giving  rise  to  a  uniform  roughness  in  the  reflected  wavefront;  i.e., 
a  definite  radar  angular  spectrum  crQ{tp),  uniform  throughout  the  corona. 

B.  Time  Delay  (Range) 

For  soft  targets,  contours  of  equal  time  delay  will  no  longer  necessarily  be  contours  of 
equal  range  because  the  group  propagation  velocity  in  the  plasma  may  be  less  than  the  free  space 
value.  In  particular,  the  total  group  time  delay  T  will  be  twice  the  reciprocal  of  the  integrated 
product  of  group  velocity  and  distance;  i.e., 

T^wislds  . 


where  w  is  the  group  velocity  and  R  is  the  range  of  the  point  of  reflection.  In  general,  the 
medium  will  be  dispersive;  that  is,  w  will  depend  on  frequency.  In  the  important  case  for  which 
both  collision  frequency  and  gyro  frequency  are  negligible  compared  to  the  operating  frequency  f, 

w(s)  =  nc  *  c  Ji  -  [fp(s)/f]2  , 

where  n  is  the  index  of  refraction  and  f  is  the  local  plasma  frequency  in  cps; 
fp(s)  =  9.0  VN(sj  , 

N  being  in  electrons  per  cubic  meter.  If  f  «  f  also,  then  we  have  the  useful  formula 


2 
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which  expresses  the  excess  of  time  of  flight  over  that  is  vacuo  in  terms  of  integrated  electron 
density  along  the  propagation  path  (number  of  electrons  in  a  one-square-meter  column  travers¬ 
ing  the  radar-to-target  path).  • 


C.  Frequency 

We  have  just  referred  to  the  fact  that  for  a  soft  target,  both  the  reflected  power  and  the  prop¬ 
agation  velocity  can  be  frequency-dependent,  and  as  is  well  known,  the  behavior  of  these  quanti¬ 
ties  (and  polarization  too)  can  under  certain  conditions  change  violently  with  frequency.  This 
makes  it  possible,  in  principle,  to  deduce  certain  parameters  in  the  corona,  or  in  planetary 
ionospheres,  by  means  of  measurements  at  different  frequencies.  This  sort  of  technique  has 
been  richly  cultivated  for  years  in  the  study  of  our  own  ionosphere  (see,  for  example.  Ref.  19) 
whose  various  layers  are  studied  by  noting  the  frequency  below  which  a  test  signal  no  longer 
passes  through  but  is  refracted  back. 

It  is  clear  that  gross  shifts  in  the  frequency  of  the  radar  echo  are  to  be  associated  with  radial 
components  of  motion  in  the  soft  target.  For  example,  parts  of  the  solar  corona  having  the  same 
range  but  different  velocities  relative  to  the  radar  line-of-slght  may  thereby  be  distinguished  one 
from  the  other  much  as  with  various  portions  of  hard  targets  (Secs.  II-C  and  II-D).  However,  the 
lack  of  target  rigidity  destroys  much  of  the  usefulness  of  range-doppler  echo  resolution. 
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D.  Polarization 


The  way  in  which  the  polarization  characteristics  of  the  incident  radar  signal  are  altered  by 
passage  through  a  plasma  in  the  presence  of  a  magnetic  field  Is  quite  complicated  in  the  general 
case.  If  there  is  no  magnetic  field,  or  if  it  is  the  incoherent  scatter  mode  that  is  being  employed, 
then  the  polarization  is  unaffected;  it  is  as  though  the  reflection  were  from  a  smooth  hard  target. 

When  a  magnetic  field  is  present  in  the  ionized  medium,  and  for  the  common  case  of  operat¬ 
ing  frequency  considerably  larger  than  both  gyro  and  plasma  frequencies,  then  propagation  takes 
place  via  two  independent  modes  having  opposite  senses  of  circular  polarization  and  slightly  dif¬ 
ferent  phase  velocities  (neglecting  collisions).  When  present  simultaneously,  these  two  waves, 
the  "ordinary"  and  "extraordinary,"  form  an  elliptlcally  polarized  wave  whose  ratio  of  axes  and 
sense  of  rotation  do  not  change  along  the  ray  path,  but  whose  plane  of  polarization  does,  because 
of  the  different  phase  velocities.  The  rotation  is  cumulative  in  the  sense  that  ii  ;ontinues  on  the 
return  trip  from  the  target  rather  than  undoing  itself.  The  total  number  of  radians  Ax  of  this 
Faraday  rotation  is20 


N(s)  B8(s)  ds 


(20) 


where  Bg(a)  is  the  longitudinal  component  of  magnetic  field  (in  gauss).  Note  that  there  is  a  dis¬ 
persion  (strong  frequency  dependence)  of  the  total  amount  of  rotation.  In  the  special  case  in 
which  the  incident  radiation  is  linearly  (plane-)  polarized,  it  can  be  considered  the  sum  of  equal- 
strength  right-hand  and  left-hand  circularly  polarized  components.  One  propagates  by  the  or¬ 
dinary  mode  only,  and  the  other  by  the  extraordinary  mode  only,  so  that  the  composite  signal 
arriving  is  still  linearlv  polarized  -  if  not  depolarized  on  reflection  (see  Sec.  II-E)  -  but  with  the 
polarization  axis  twisted  by  Ax  radians. 

An  Ideally  smooth  target  converts  linear  polarization  into  linear,  and  right  circular  into  left 
circular,  and  vice  versa.*  In  preference  to  linear  polarization,  one  often  uses  circular  trans¬ 
mitted  polarization  in  conjunction  with  a  receiving  antenna  adjusted  for  the  opposite  direction  of 
circular  polarization  so  that  the  signal  power  received  from  a  smooth  hard  target  can  be  maxi¬ 
mized  independently  of  the  Faraday  rotation  angle  Ax- 


TV.  MEASUREMENT  ME1 HODOLOGY 

In  the  next  few  sections  of  this  report,  the  way  in  which  radar  systems  perform  measure¬ 
ments  on  radar  targets  will  be  systematized.  In  Sec.  V  we  shall  collect  all  the  preceding  notions 
and  list  the  way  in  which  quantities  of  astronomical  interest  can  be  deduced  from  these  radar 
measurements. 

The  targets  dealt  with  in  radar  astronomy  differ  from  those  discussed  in  much  of  the  stand¬ 
ard  radar  literature  by  their  sheer  physical  size  and  large  differential  velocities.  Very  rarely 
will  these  targets  be  extended  targets,  meaning  that  they  more  than  fill  the  antenna  beam  be¬ 
cause  they  are  usually  at  extreme  ranges;  but,  they  will  be  what  we  shall  call  spread  targets 


tThe  (normalized)  Mueller  matrix  of  Eq.  (17)  1$ 
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in  both  delay  and  frequency.  That  la,  since  they  are  so  deep  and  since  different  portions  are 
apt  to  be  moving  differently  relative  to  the  observer,  the  returns  are  spread  significantly  in 
both  range  and  doppler. 

If  neither  the  range  spread  nor  the  doppler  spread  is  discernible  to  the  radar  system  in  use, 
we  call  the  target  a  point  nonfluctuating  target.  If  the  bandwidth  of  the  signal  is  insufficient  to 
resolve  the  target  in  range,  and  yet  there  is  a  discernible  doppler  spread  of  the  received  signal, 
we  shall  call  the  target  a  point  fluctuating  target.  If  the  signal  resolves  the  target  in  range  but 
with  no  observable  doppler  spreading  we  can  treat  the  target  as  a  deep  nonfluctuating  target.  In 
the  general  case  in  which  observable  range  spreading  and  doppler  spreading  are  present  simul¬ 
taneously  we  speak  of  a  deep  fluctuating  target. 

A.  The  Target  Scattering  Function  <r(r,  f) 

It  has  proved  convenient  to  describe  the  average  reflecting  properties  of  a  spread  radar  tar¬ 
get  in  terms  of  its  scattering  function  <j(t,  f),  the  average  reflectivity  aB  a  function  of  the  two 
variables  delay  and  doppler  shift.  The  scattering  function  can  thus  be  depicted  as  a  surface  in 
three  dimensions,  as  in  Fig.  9.  The  two  Independent  variables  are  delay  r(seconds)  and  doppler 
shift  f  (cycles  per  second).  The  dependent  variable  o(r,  f)  drdf  is  defined  as  that  part  of  the  total 
target  cross  section  o(in  square  meters)  that  belongs  to  target  regions  that  both  reflect  a  signal 
with  delays  of  r  to  t  -t  dr  and  at  the  same  time  subject  the  signal  to  doppler  shifts  from  f  to  f  +  df. 
A  more  concise  discussion  of  the  scattering  function  follows  in  Sec.  IV-D. 

The  concept  of  a  scattering  function  has  proved  useful  in  unifying  various  techniques  of  meas¬ 
urement  (this  report)  and  detection.1  Many  of  the  measurement  techniques  we  shall  describe 
shortly  are  nothing  more  than  measurements  (i.e.,  estimates)  of  the  scattering  function  or  some¬ 
thing  derivable  from  it.  In  Ref.  1  it  is  found  that,  under  certain  reasonable  assumptions,  all  one 
needs  to  know  about  the  target  in  order  to  discuss  its  detectability  completely  is  its  scattering 
function. 

Several  cases  are  shown  in  Fig.  9.  For  example,  a  stationary  point  target  at  delay  Td  would 
have  the  scattering  function  depicted  in  Fig.  9(a),  namely  an  impulse  of  value  a  at  r  =  rd  and 
f  =  0.  If  the  point  target  were  moving  toward  the  observer  at  a  velocity  v  =  cfd/2fo  (fQ  being  the 
carrier  frequency)  the  scattering  function  would  appear  as  shown  in  Fig.  9(b).  If  the  target  (at 
delay  rd  and  moving  at  velocity  cfd/2fQ)  could  be  regarded  as  having  negligible  depth  in  delay, 
but  as  producing  rapid  random  fluctuations  in  the  returned  echo,  the  scattering  function  o(r,  f) 
would  have  the  character  shown  in  Fig.  9(c)  -  there  would  be  a  spread  only  along  the  f-axis.  The 
width  B  along  the  f-coordinate  is  proportional  to  the  fluctuation  rate.  Similarly,  if  the  fluctua¬ 
tion  rate  were  negligible,  but  the  target  were  still  quite  deep  and  thus  had  a  wide  range  spread, 
the  scattering  function  would  appear  as  in  Fig.  9(d).  The  duration  of  delays  or  "multipath  spread" 
is  labeled  L. 

Figure  9(e)  shows  the  scattering  function  of  a  hypothetical  target  having  a  significant  spread 
along  both  axes.  Figure  9(f)  shows  a  typical  such  scattering  function  for  the  important  special 
case  of  a  uniform  rotating  rough  sphere.  As  mentioned  earlier,  a  uniform  sphere  is  one  having 
the  same  radar  angular  spectrum  vQ(<p)  everywhere  on  its  surface.  Appendix  B  presents  a  deri¬ 
vation  of  the  following  equation  for  c(r,  f)  of  the  uniform  sphere  of  radius  a,  rotation  speed  0 
and  angle  a  between  polar  axis  and  perpendicularity  to  the  line  of  sight: 
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{•)  Deep  fluctuating  target. 


(f)  Uniform  rotating  iphara. 


Fig.  9.  Scattering  function  of  various  forms  of  targets  (the  total  volume 
under  the  function  Is  equal  to  e,  the  over-all  radar  cross  section). 
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where  f  is  the  carrier  frequency  and  c  is  propagation  velocity.  The  T-orlgin  is  at  the  target 
center. 

It  will  be  recalled  from  Sec.  II-D  that,  because  of  the  mapping  property  of  rotating  hard  tar¬ 
gets,  points  in  the  r,  f  plane  correspond  to  locations  on  the  surface  of  the  target.  Therefore  a 
measurement  of  the  scattering  function  a(r,  f)  represents  a  map  of  reflected  power  from  various 
regions  of  the  target  with  the  correspondence  being  that  given  in  Pig.  4.  [Note  the  occurrence 
in  Fig.  9(f)  of  the  same  semicircular  or  semielliptical  region  in  the  r,  f  plane  as  is  shown  in 
Fig.  4(b).] 

The  scattering  functions  shown  in  Fig.  9  have  been  labeled  to  indicate  roughly  the  two  param¬ 
eters  B  and  L,  the  doppler  spread  and  the  delay  spread,  respectively  [for  9(c)  we  have  L  =  0  and 
for  9(d),  B  =  0].  An  important  attribute  of  spread  radar  targets  is  the  spread  factor,  the  product 
BL.  If  BL  <  1  the  target  is  termed  underspread,  and  if  BL  >  1,  overspread.  As  pointed  out  in 
Ref.  1,  the  radar  detectability  of  a  target  of  given  total  cross  section 
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tr( r,  f)  dr  df 


(the  volume  under  the  surface  depicted  in  Fig.  9)  decreases  steadily  with  BL  if  BL  becomes 
greater  than  roughly  unity.  In  the  neighborhood  of  BL  =  1,  even  when  the  over-all  system  sen¬ 
sitivity  is  enough  to  overcome  this  loss  in  detectability,  we  shall  see  in  Sec.  VI-D  of  this  report 
that  it  may  became  difficult  to  perform  an  estimation  of  the  scattering  function  or  derived  param¬ 
eters,  a  necessity  for  studying  target  properties  in  detail. 

In  order  to  discuss  processes  of  estimating  the  scattering  function,  or  derived  properties, 
we  must  make  a  short  digression  to  Introduce  the  notions  of  a  matched  filter  and  the  associated 
ambiguity  function.  (For  a  more  detailed  discussion  see  the  paper  by  Turin.21) 

B.  Detection  of  a  Point  Nonfluctuating  Target  -  the  Matched  Filter 

The  simplest  form  of  target  is  one  whose  scattering  function  is  an  impulse,  i.e.,  a  point 
nonfluctuating  target.  As  is  well  known,  the  form  of  receiver  that  will  develop  maximum  signal- 
to-noise  ratio  from  such  a  target  is  a  matched  filter.  (Most  conventional  radar  receivers  employ 
matched  filtering  in  the  IF.)  That  is,  denoting  the  transmitted  waveform  as  x(t),  if  the  noise 
spectrum  is  uniform,  the  optimum  receiver  consists  of  a  filter  having  an  impulse  response 
m(t)  =  x(to  -  t)  followed  by  a  sampling  device  that  observes  the  filter  output  at  time  tQ  (where  tQ 
is  an  irrelevant  time  delay  introduced  to  make  the  filter  physically  realizable,  i.e.,  no  response 
before  excitation).  The  matched  filter  has  an  impulse  response  that  is  a  time-reversed  version 
of  the  signal.  Since  its  complex  frequency  function  is  thus  the  complex  conjugate  of  the  frequency 
function  of  the  signal,  a  matched  filter  is  often  referred  to  as  a  conjugate  filter. 

It  will  be  relatively  simple  and  quite  instructive  to  prove  that  this  matched  condition  maxi¬ 
mizes  the  signal -to-noise  ratio  as  follows.  Let  X(f)  be  the  Fourier  transform  of  the  signal 
x(t)  (i.e.,  its  complex  spectrum),  and  M(f)  be  the  complex  frequency  function  of  the  filter.  We 
want  to  adjust  M(f)  so  as  to  maximize  the  output  slgnal-to-noise  ratio,  defined  as  the  ratio  of 
Sout'  the  s9uare  of  the  output  voltage  at  time  tQ  due  to  signal  alone,  divided  by  NQUt,  the  output 
noise  power.  The  latter  is 
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for  input  noise  of  uniform  spectral  density  Ng/2  watts  per  cycle  per  second  (Nq  watts  per  cps 
for  a  single-sided  spectrum).  The  signal  output  is 
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X(f)  M(f)  exp(j2irfto  df 


|X(u)  exp[j2irutoJ 


du 


j”6  |M(v)  j  ‘ 


dv 


(23) 


from  the  Schwarz  inequality  for  two  complex  quantities  A  and  B: 


A(x)  B(x)  dxj  ^  y  |  A  (y)  | 2  dy  j* 


|B(z)|2dz 


Clearly,  in  order  for  the  inequality  (23)  to  be  an  equality,  we  must  have 


M(f)  ={X(f)  exp  [j2xfto)}  *  =  X*(f)  exp  [-j2irft0] 


or  in  terms  of  the  impulse  response  m(t), 
m(t)  =  x(tQ  -  t)  , 


(24) 


which  was  to  be  proved.  When  this  matched  condition  exists,  Eqs.  (22)  and  (23)  combine  to  give 
for  the  signal-to-noise  ratio 


(25) 


where  Efl  is  the  total  signal  energy  |X(f)|2  df.  This  ratio  2Eg/NQ  is  an  important  factor  in 
the  detectability  of  any  signal. 

The  output  voltage  at  time  tj  produced  from  any  filter  of  impulse  response  m(t)  due  to  an 
input  w(t)  is  given  by  the  convolution 
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w(t)  mft^  -  t)  dt 


which  for  the  matched  condition  is,  by  Eq.  (24), 

VV  =  J  w(t)  x(r  -  tj  +  tQ)  dr  .  (26) 

From  this  last  equation,  we  see  that  the  operation  of  observing  the  matched  filter  output  at  a 
given  instant  is  mathematically  equivalent  to  correlation  detection,  i.e.,  a  multiplication  of  the 
incoming  signal-plus-noise  mixture  w(t)  with  a  reference  copy  of  the  transmitted  signal  x(t)  fol¬ 
lowed  by  integration  over  the  duration  of  the  product  signal.  We  have  already  met  correlation 
operations  in  connection  with  Eq.  (7),  which  deals  with  a  complex  correlation  function.  Wr  will 
return  to  the  question  of  correlation  operations  shortly. 

2 

C.  The  Ambiguity  Function  $ 

This  matched  filtering  or  correlation  detection  operation  is  optimum  if  the  scattering  func¬ 
tion  has  the  form  of  an  impulse  located  at  the  origin.  If  the  impulse-like  scattering  function  is 
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not  at  the  origin  but  at  the  point  (r,  f),  the  output  is  in  general  diminished,  and  tha  filter  must  be 
readjusted  by  the  appropriate  frequency  offset  f  and  its  output  observed  r  seconds  later  to  re¬ 
store  the  condition  of  optimality. 

Suppose  we  plot  the  signal  output  power  from  a  matched  filter  set  for  such  a  point  non¬ 
fluctuating  target  at  r  =  0  and  f  =  0  as  a  function  of  delay  and  doppler  misalignments  r  and  f. 
The  result  is  the  quantity  <^( r ,  f)  called  the  ambiguity  function22,23  (or  the  uncertainty  function, 
resolution  function,  or  two-dimensional  autocorrelation  function) .  More  explicitly,  the  ambiguity 
function  |t2(r,  f)  is  defined  as  simply  |/0(tQ)  | 2,  the  square  of  the  envelope  of  the  output  fQ(t)  of  a 
matched  filter  at  t  =  tQ  to  an  input  waveform  x(t)  that  has  been  misaligned  by  a  lag  of  r  seconds 
and  an  upward  shift  in  frequency  f  cycles  per  second: 


=  l/.e*c 
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X*(t)X  (t  +  t)  exp[j2*ft)  dt 


(27) 


or,  alternatively. 


*2<r,f)  = 


G((p)  G*(p  +  f)  exp[-j2»ipT]  dip 
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where  Q(f)  is  the  Fourier  transform  of  X(t),  (i.e.,  G  (f)  <=  9C(t)).  As  before,  t  is  the  time 
at  which  one  optimally  samples  the  matched  filter  output  if  r  =  0  and  f  =  0.  The  ambiguity  func¬ 
tion  expresses  the  ability  of  a  given  combination  of  waveform  x(t)  and  filter  matched  to  it  to  dis¬ 
criminate  targets  in  range  and  doppler.  Ambiguity  functions  for  some  common  sorts  of  wave¬ 
forms  are  shown  in  Fig.  10. 

The  ambiguity  function  has  a  number  of  interesting  properties.  The  width  in  r  at  f  =  0  is 
roughly  the  reciprocal  of  W,  the  bandwidth  of  the  signal;  and  the  width  in  f  at  r  =  0  is  roughly 
the  reciprocal  of  the  signal  duration  T.  If  the  energy  in  the  signal  and  the  gain  of  the  matched 
filter  are  normalized  so  that 


j*  |X(f)|2=J  | M(f) | 2  df  =  1  , 

then  not  only  do  we  find  readily  from  Eq.  (27)  that  the  height  of  the  ambiguity  surface  at  the 
origin  is  unity 

«2(0, 0)  =  1  ,  (29) 

but  also  that  the  total  volume  under  the  surface  is  invariant  to  the  form  of  signal, 

I"  I"  df  dT  “  1  '  (30) 

The  "invariance"  relation  [Eq.  (30)]  is  highly  significant.  The  diagrams  in  Fig.  10  show  the 
ambiguity  functions  of  several  es  of  waveforms.  In  all  cases  the  height  at  the  origin  is  unity 
[Eq.  (29)],  and  the  volume  is  also  unity  [Eq.  (30)};  the  difference  lies  in  how  this  volume  is  dis¬ 
tributed.  The  first  three  examples  shown  concentrate  substantially  all  of  the  volume  in  one  cen¬ 
tral  peak.  Thus,  lengthening  a  simple  pulse  produces  a  broadening  in  delay  which  would  be  ir¬ 
relevant  if  only  the  doppler  were  being  measured  [Fig.  10(a)];  or  narrowing  a  simple  pulse  in 
time  [Fig.  10(b)]  produces  a  broadening  in  frequency  which  would  be  irrelevant  if  one  were  meas¬ 
uring  range  only. 
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(d)  Lorgo-TW  pul**  train. 


(•)  toudonoita  (In  TW  cantinuout 
waveform). 


Fig.  10.  Ambiguity  function!  of  typical  gouwion-ihapad  pwlw  wovcfomn. 


The  lest  two  examples  in  the  figure  repreeent  deliberate  efiorta  to  make  both  dimension*  of 
the  central  peak  small.  This  sllows  range  and  doppler  resolution  to  be  obtained  simultaneously. 

In  the  case  of  the  "large-TW  pulse  train"  [Fig.  10(d) ),  the  rest  of  the  volume  goes  into  subsidiary 
peaks.  For  the  continuous  large-TW  waveform  or  "pseudonoise"  waveform  [Fig.  10(e)].  the 
volume  is  mostly  in  the  "skirts"  away  from  the  central  peak.  A  great  deal  of  study  has  been 
given  to  ways  of  making  these  skirts  even,  so  that  there  are  no  large  subsidiary  peaks  (e.g., 

Ref.  24). 

Signals  of  the  type  shown  in  Figs.  10(c)  through  (e)  are  useful  for  another  reason  besides  the 
svallability  of  simultaneous  doppler  and  range  resolution.  For  a  given  peak  power  limitation  on 
the  transmitter,  these  waveforms  allow  an  increase  in  signal  energy  Eg  by  leaving  the  trans¬ 
mitter  on  longer,  without  a  sacrifice  in  range  resolution  as  would  be  the  case  with  the  waveform 
in  Figs.  9(a)  and  (b). 

In  estimating  the  performance  of  the  various  waveforms  of  Fig.  10  it  will  be  convenient  later 
on  to  use  the  crude  approximations  given  in  the  right-hand  diagrams  of  Fig.  10.  There  the  volume 
in  a  given  peak  is  assumed  to  be  concentrated  uniformly  within  the  half -power  contours,  and  in 
the  case  of  the  pseudonoise  waveform  the  skirt  level  is  assumed  uniform. 

D.  Interaction  of  the  Signal  and  the  Target  in  Terms  of  a  and  p 

We  have  now  introduced  two  functions  of  delay  r  and  doppler  f,  (me  being  the  target  scatter¬ 
ing  functior  <r(r,  f),  a  property  of  the  target  only,  and  the  ottu.r  being  the  ambiguity  function 
*Z(T,f)  which  describes  the  transmitted  signal  independently  of  target  properties.  These  have 
bee-'.  ...us+raie  3  in  ."vs.  sa-i  iti,  respectively.  We  biiall  now  obtain  in  expression  involving 
these  *.vro  quantities  that  shows  how  a  simple  matched  filter  radar  sys  em  responds  to  a  spread 
target.  This  expression25  is  a  convolution  of  o  and  p2.  For  the  more  complicated  radar  sys¬ 
tems  that  provide  optimum  detectability  on  spread  targets1  other  expressions  are  obtained,  but 
they,  too.  are  found  to  take  into  account  the  target  properties  only  through  the  ambiguity  func¬ 
tion,  and  to  involve  the  same  convolution  of  r  and  p2,  as  we  shall  see. 

Imagine  that  we  have  an  ensemble  of  statistically  identical  targets .  With  a  transmitted  sig¬ 
nal  of  complex  envelope  'X  (t)  emerging  from  the  transmitter,  a  signal  s(  r\ V)  X  (t)  represents 
the  random  complex  echo  signal  returning  from  the  elementary  region  having  a  span  of  delays 
lying  between  r'  end  r’  +  dr'  seconds,  and  at  the  same  time  having  a  r  *nge  of  doppler  frequency 
shifts  lying  between  V  asvd  f*  +  df  cycles  per  second.  The  target  scattering  function  ofr'.f')  is 
defined  as 


<r(r',f‘)  -  ;*<t',!’)|2  . 

where  the  bar  represents  toe  ensemble  average,  and  parameters  such  as  range,  antenna  gain, 
etc.,  are  normalized  out  so  mat  o,(r',  V)  is  in  terms  of  the  number  of  effective  square  meters  of 
cross  section  per  second  of  range  increment  per  cps  of  doppler  increment.  That  is,  on*  watt 
per  square  meter  of  incident  power  flex  will  scatter  toward  the  receiv.tr  [ff(T',  f)  dr'df  ^4*  watts 
per  steradian  of  solid  angle  subtent !i  d  at  the  target.* 

t  in  determining  the  Kettering  (unction  -  ’  soft  targets  such  as  clouds  of  scattarers,  :  t  i*  useful  to  consider 
•(V.P)  to  be  the  average  cross  section  t-  a  typical  particle  times  the  probability  H  at  such  e  particle  will  fulfill 
Ihe  condition  that  the  range  Is  In  the  ini  vrval  V  to  V  i  dV  and  the  velocity  lias  In  P  to  P  ♦  dP. 
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The  eaeiest  way  to  deduce  the  effect  of  the  apread  target  on  the  radar  aignal  is  firat  to  divide 
the  target  up  into  zonea  having  different  delays.  Suppose  we  have  isolated  such  a  region  of  the 
target  containing  all  points  lying  within  a  range  of  path  delays  r'  to  t1  +  dr'.  We  next  Imagine 
that  the  rest  of  the  target  is  absent  and  we  compute  the  contribution  of  this  zone  to  the  output  of 
the  matched  filter.  If  those  regions  within  the  zone  that  introduce  a  frequency  offset*  lying  in 
the  range  f  to  P  +  df  return  an  echo  s(t’,  f ')  9C(t),  then  the  envelope  of  the  output  from  a  filter 
that  is  matched  to  the  transmitted  waveform  (except  that  the  filter  is  misaligned  upward  by  f  in 
frequency  and  the  output  is  observed  at  t  =  tQ  +  t  instead  of  at  t  =  t0)  will  be 

i/0<yl  =  |s(t*,  f*)|!  J  X*(t>*  (t  +  r'  -  t)  exp  [j2»(f  -  f’)tj  dt|  »  P(t'  -  t,  f  -  f )  |s(t*,  f>)  |  , 

since  the  signal  arrives  with  a  total  time  offset  (t1  -  t)  and  total  frequency  offset  (f  -  f)  relative 
to  the  values  giving  maximum  filter  output.  Use  has  been  made  of  the  relation  J AB  |  =  [ A {  |B| 
for  two  complex  quantities  A  and  B.  The  average  power  [P(t,  f)JT,  appearing  in  the  output  of 
the  filter  offset  by  r  and  f  due  to  the  entire  zone  having  delay  near  r'  will  then  be  gotten  by 
integrating  the  average  of  the  square  of  this  quantity  over  all  frequency  (contributions  at  differ¬ 
ent  frequencies  always  adding  as  powe  i,  not  as  voltages): 

(P(r,  f))Tt  ■  l/0(to>|2=  J"  02(t'  -  T,  f-f’)  |s(T',f')|2  df 
*  j*  *&2(t'-T,  f-f’)  (7(T',f)  df’  . 

As  a  final  step,  we  add  up  the  contributions  due  to  all  range  zones  of  the  target.  If  we  assume 
that  the  echoes  from  different  range  elements  are  uncorrelated;*  i.e.,  |s(t^,IJ)  s(t^,  fy|  =  0 
for  all  r'j  ^  then  we  can  add  together  the  powers  returned  in  the  different  range  intervals  to 
*et 

P(r,f)  -  J  (P(T.fHT,dT' -  f-f'Mr'.ndT' df*  .  (31) 

Thus  the  power  out  of  the  matched  filter  as  a  function  of  the  time  and  frequency  offsets  is 
proportional  to  the  (two-dimensional)  convolution  of  the  scattering  function  with  the  ambiguity 
function.*  Equation  (31)  says  that  the  process  of  matched  filter  reception  isolates  certain  areas 


fThe  frequency  offset  might  occur  become  at  the  given  range  there  are  region*  In  motion  with  the  appropriate 
velocity.  On  the  other  hand,  a  frequency  apreod  of  the  energy  returning  from  a  given  range  could  just  a*  easily 
occur  as  a  result  of  other  physical  mechanisms.  For  example,  at  a  given  range  them  might  be  o  single  fixed 
scatterer  whose  reflectivity  is  fluctuating  rapidly  enough  to  cause  the  assumed  doppier  broadening. 

f  As  noted  later  at  the  close  of  Soc.V-E,  it  is  conceivable  that  the  signal  returning  from  the  target  in  one  cell  of 
the  t*,?  plane  might  be  partially  coherent  with  that  appearing  in  a  nearby  cell  at  a  later  time.  In  such  cases 
this  step  cannot  strictly  be  taken.  The  whole  question  of  such  cel  I -to -cel I  coherences  remains  to  bo  studied  In 
the  radar  astronomy  context. 

$  In  order  to  avoid  a  weighting  proportional  to  the  antenna  gain  pattern,  we  have  Implicitly  assumed  that  the  an¬ 
tenna  beam  width  Is  much  greater  than  the  angle  subtended  by  the  target  (or  equivalently,  for  targets  broader 
than  the  beam  width,  extended  targets,  that  the  antenna  radiates  uniformly  throughout  a  certain  solMangle,  and 
zero  elsewhere).  Also,  it  should  be  added  that  whenever  the  filter  and  transmission  ora  not  matchod,  f2  It  simply 
replaced  by  the  slightly  more  general  "crow  ambiguity  function.”22 
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of  the  scattering  function.  The  process  can  be  visualized  as  a  process  of  overlaying  the  three- 
dimensional  surface  +2(t,  f),  the  ambiguity  function  of  Pig.  10,  onto  a  similar  representation  of 
the  scattering  function  o(r,  f)  of  Fig.  9  (after  providing  the  appropriate  r  and  f  offsets)  and  then 
reading  off  the  total  volume  under  the  product  of  these  two  functions. 

Let  us  now  investigate  the  various  forms  of  receivers  that  are  actually  uaed  in  radar  astron¬ 
omy  work  where  adequate  signal-to-nolse  ratio  is  always  a  problem.  As  detailed  by  R.  Price  in 
the  companion  report,1  in  order  to  build  up  the  signal-to-nolse  ratio,  one  usually  builds  a  re¬ 
ceiver  structure  that  goes  beyond  simple  matched  filtering  or  simple  correlation  detection  as 
just  discussed  in  Sec.  IV-B. 

The  matched  filter  with  a  sampling  of  the  output  voltage  at  is  appropriate  only  for  a 
nonfluctuating  point  target  for  which  one  knows  target  delay  r  to  within  a  traction  of  an  RF  cycle. 
Since  this  means  a  fraction  of  a  wavelength,  in  practice,  carrier  phase  is  almost  always  unknown, 
and  so  one  samples  instead  the  envelope  at  t^.  Figure  11(a)  shows  two  alternate  versions  of  the 
matched  filter  for  known  RF  phase,  the  matched  filter  operation  itself  (as  we  have  described  it 
in  Sec.  IV-B),  and  the  correlation  detection  equivalent  (£q.  24),  The  fitter  and  ecnfcsa- 

lent  correlation  detector  of  Fig.  14(b)  are  the  receivers  of  Fig.  11(a)  suitably  modified  for  opti¬ 
mality  when  RF  phase  is  unknown.  In  the  case  of  die  matched  filter  shown  on  the  left,  the  en¬ 
velope  of  the  output  is  sampled,  and  in  the  correlation  detector  shown  at  the  right,  the  envelope 
of  the  (complex)  cross  correlation  of  arriving  signal  w(t)  and  transmission  x(t)  is  observed, 
rather  than  only  the  real  part  as  in  Eq.  (24).  The  distinction  will  be  discussed  further  in  the  next 
section. 


For  uniformity  with  the  other  more  complicated  detectors  which  we  are  about  to  discuss. 

Fig.  11  (b)  shows  the  square  of  the  envelope  being  sampled,  not  the  envelope  itself.  For  a  point 
nonfluctuating  target,  the  output  power  signal-to-nolae  ratio  is  (E,/N0)2  as  compared  with 
2(Ea/N0)2  available  when  carrier  phase  is  known  exactly.  In  terms  of  required  signal  energy 
E,.  the  lack  of  phase  knowledge  coats  3  db.  For  a  more  general  form  of  target,  the  output  power 
signal-to-noiae  ia 


p  i0r(r.  f)  »  Const  x 


r,  f-f)  e(r'.  A)  dr'  df'} 


(32) 


In  radar  astronomy  experiments  the  duration  T  of  the  transmission  is  usually  very  long 
compared  to  the  received  signal  fluctuation  period  1/B,  ao  that  one  must  regard  the  target  ar, 
being  a  point  fluctuating  target  (spread  in  frequency  f),  so  that  its  scattering  function  ia 

«r(r,f)  «  e(f)  «(t- t')  , 

where  6  ia  the  unit  impulse.  (The  condition  T  »  1/B  ia  often  unavoidable  irrespective  of  the 
target  properties  simply  because  of  limits  on  the  achievable  frequency  stability  of  the  transmit¬ 
ting  and  receiving  conversion  stages.)  The  usual  form  of  receiver  ia  then  the  one  shown  in 
Fig.  11(c),  where  the  received  signal  w(t)  is  premultiplied  or  "weighted"  by  a  replica  of  the  trans¬ 
mission  x(t)  to  undo  any  phase  modulation  imposed  at  the  transmitter  and  to  accentuate  those 
portions  of  the  received  signal  w(t)  during  which  one  expects  greater  amplitude;  the  weighting 
ia  followed  by  the  familiar  predetection  filtering,  rectification,  and  postdetection  integration 
operations.  In  the  general  theory  of  signal-to-noiae  ratio  optimization  developed  by  Price1  and 
summarized  here,  such  a  receiver  ia  called  a  weighted  radiometer  because  the  three  operations 
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(a)  Matched  filter  (left)  and  correlation  detector  (right)  for  ca»e  of  known  RF  phase. 


W  Matched  filter  (left)  and  eorrelaHen  detector  (right)  for  RF  phase  unknown. 


(c)  Weighted  radiometer,  for  point  fluctuating  targets. 


(d)  Rake  radiometer,  for  deep  fluctuating  targets. 

Fig.  II.  Various  signal-to- noise  ratio  maximizing  receivers. 
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following  the  weighting  are  precisely  thosn  used  in  a  radiometer  such  as  is  employed  in  radio 
astronomy  to  develop  maximum  signal-to-noise  ratio  from  an  incoherent  narrow-band  source 
such  as  the  hydrogen  line.  For  maximum  signal-to-noise  ratio, t  it  turns  out  that  the  predetec¬ 
tion  filter  characteristic  H(f)  must  be  chosen  to  match  the  power  profile  of  the  line,  or  in  our 
context, 

| H(f)  | 2  =  Const  x  a(f)  .  (33) 

For  the  general  case  of  a  target  that  may  or  may  not  be  spread  in  both  r  and  f,  and  for  an  ar¬ 
bitrary  H(f),  the  following  expression  gives  the  output  signal-to-noise  ratio  for  the  weighted 
radiometer 


PWp(T,  f)  =  Const  x 


| H(f •  -f)|2  df  #>2(t'  -  t,  V  -f")  «r(T',f")  dr  df" 


•  (34) 


(As  before,  r  and  f  are  the  offsets  inserted  in  the  receiver  processing,  usually  in  the  timing 
and  frequency  shift  of  the  weighting  operation.)  That  is,  the  convolution  appearing  in  Eqs.  (31) 
and  (32)  is  itself  smoothed  in  turn  frequency-wise  by  the  squared  magnitude  of  the  predetection 
filter  frequency  characteristic  H. 

The  theory  developed  in  Ref.  1  goes  on  to  say  that  if  the  target  is  significantly  spread  not 
only  in  frequency  f  but  also  in  delay  r  (deep  fluctuating  target),  then  the  optimum  receiver  is 
the  rake  radiometer  of  Fig.  11(d),  which  is  obtained,  plausibly  enough,  by  iterating  the  weighted 
radiometer  of  Fig.  11  (c)  at  a  number  of  delays  successively  spaced  by  a  small  At,  the  total  range 
of  delays  being  large  enough  to  at  least  span  L,  the  duration  in  r  of  the  target  scattering  func¬ 
tion.  For  signal-to-noise  ratio  maximization,  it  will  turn  out  that  the  frequency  functions  of  the 
predetection  filters  H^f)  =  H(r,  f)  (where  t.  =  some  tq  +  iir)  must  obey  the  condition 


|H(T,f)|2  =  Const  xqr(r,f)  .  (35) 

That  is,  the  filter  squared  magnitudes,  when  read  off  along  the  succession  of  delay-line  taps 
(and  according  to  frequency  at  a  given  tap)  must  reproduce  the  target  scattering  function.  For 
the  more  general  case  of  arbitrary  H(t,  f),  the  following  expression  gives  the  output  signal-to- 
noise  ratio  for  the  rake  radiometer: 


PRR(t,  f)  =  Const : 


-j*  |H(t'-t,  f '  — f) |  2 dr* df '  f'-f")  <r(T",f")  dt"df" 


(36) 


That  is,  the  original  convolution  of  Eq.  (31)  undergoes  a  two-dimensional  convolution  with  the 
predetection  filter  power-vs-frequency  functions  read  off  in  t. 

Thus  we  see  that  for  the  entire  class  of  receivers  considered  to  be  of  value  from  the  signal- 
to-noise  ratio  point  of  view,  the  convolution  Eq.  (31)  is  the  central  relationship.  To  see  how  the 
two-dimensional  convolution  operation  works  out  in  practice,  consider  convolving  the  scattering 
functions  of  Fig.  9  with  the  ambiguity  functions  of  Fig.  10.  In  probing  a  target  spread  in  both 


tA»  it  pointed  out  in  Ref.l,  the  retulti  quoted  Here  on  the  form  of  receiverthot  it  optimum  ore  strictly  true  only 
for  the  reotonable  condition  thot  the  receiver  Input  signal -to-noite  ratio  It  smaller  then  unity  In  a  certain  sense. 
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range  and  doppler  (deep  fluctuating  target),  such  as  Fig.  9(e)  or  9(f),  one  can  resolve  the  detailed 
structure  of  the  cr( r,  f)  surface  only  by  using  an  ambiguity  function  whose  central  peak  (with  lateral 
dimensions  at  l/W  by  l/T)  is  sufficiently  narrow.  As  we  have  seen  earlier,  such  studies  of  hard 
targets  are  very  interesting  because  of  the  simple  correspondence  between  places  on  the  target 
and  places  in  the  r,  f  plane  for  the  received  signal. 

However,  we  are  often  not  interested  in  measuring  the  whole  scattering  function  itself  but 
rather  certain  simpler  quantities  derivable  from  it.  Two  such  quantities  are  the  power  impulse 
response  <t(t)  and  the  echo  power  spectrum  cr(f).  These  quantities  were  briefly  alluded  to  in 
Secs.  II-C  and  n-D,  respectively.  (Two  further  quantities,  the  echo  correlation  function  and  the 
spaced -frequency  correlation  function,  are  described  in  the  next  section.)  If  we  are  observing  a 
spherical  body  with  uniform  surface  properties,  any  one  of  these  four  quantities  is  sufficient  to 
specify  the  angular  spectrum  <*Q(<p)  of  the  surface  (given  the  diameter  and  rotation  vector),  so  a 
measurement  of  the  entire  scattering  function  is  not  very  interesting  for  such  uniform  targets. 

For  targets  with  nonuniform  surfaces,  measurement  of  the  scattering  function  is  important. 
Moreover,  it  is  clear  that  the  design  of  the  optimum  detectors  of  Ref.  1  demands  as  complete  a 
knowledge  of  the  scattering  function  as  can  be  had. 

The  quantity  <r(r),  the  power  impulse  response  of  the  target,  is  a  plot  of  power  as  a  function 
of  range  only,  and  could  be  obtained  in  practice  by  using  simple  TW  =  1  pulses  [Fig.  9(b)]  of  such 
a  wide  bandwidth  as  to  resolve  all  the  interesting  variations  with  t  of  the  function  cr(r,  f)  and  also 
exceed  greatly  its  frequency  width  B.  Thus  <t(t)  is  defined  as  follows: 

or( t)  =  j*  <r( r,  f)  df  ;  (37) 

and  the  average  power  as  a  function  of  range  for  a  short  pulse*  will  be 

P(T)=j*  <PZ[t  -  r',  0)  o-(t’)  d-r'  ,  (38) 

the  T-convolution  at  ot(t)  with  the  f  =  0  profile  of  the  pulse  ambiguity  function.  According  to 
Eq.  (27),  this  profile,  in  turn,  is  the  magnitude  of  the  time  autocorrelation  of  the  pulse  complex 
envelope.  For  a  uniform  sphere,  the  power  impulse  response  is  given  by  Eq.  (13),  discussed 
earlier. 

Similarly,  the  echo  power  spectrum  o(f),  defined  as 

<r(f)  =  j"0  <r(r,f)  df  ,  (39) 

can  be  obtained  by  using  a  simple  TW  =  1  waveform  of  such  a  long  duration  T  [Fig.  9(a)]  as  to 
resolve  all  interesting  variations  with  f  of  the  function  a( r,  f).  The  average  power  as  a  function 
of  frequency  offset  for  a  long  transmitted  pulse  is  thus 

P(f)  =  f  *2(0,  f  -  f')  <r(P)  df  ,  (40) 

J-00 


fOne  occasionally  encounters  the  term  "modulation  loss"  In  the  literature  on  lunar  communications  techniques 
(e.g.,  Ref. 28).  Tills  term  refers  simply  to  the  average  received  power  obtained  with  the  nonzero  pulsewidth  In 
use,  relative  to  that  which  would  result  If  an  indefinitely  long  pulse  of  the  same  energy  were  used,  and  observed 
on  reception  at  rhe  most  favorable  delay  value. 
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the  f-convolution  of  ff(f)  with  the  r  =  0  profile  of  the  pulse  ambiguity  function.  Equation  (28)  says 
that  this  profile  is  the  frequency  autocorrelation  of  the  spectrum  corresponding  to  the  pulse  com¬ 
plex  envelope.  The  echo  power  spectrum  of  a  uniform  rough  rotating  sphere  was  presented 
earlier  as  Eq.  (14). 


E.  Spaced-Time  and  Spaced- frequency  Measurements 

There  is  another  important  way  of  measuring  the  distribution  of  echo  power  with  range  and 
doppler  (scattering  function)  besides  that  of  sending  a  signal  whose  ambiguity  function  probes 
the  scattering  function  itself.  This  method  involves  sending  two  sinusoids  spaced  by  Af  in  fre¬ 
quency  and  then  cross-correlating  (as  a  function  of  delay  At)  the  corresponding  received  sinus¬ 
oids.  That  is,  one  of  the  two  received  randomly  varying  sinusoids  is  delayed  by  At  seconds, 
readjusted  in  frequency  by  the  original  Af,  and  then  the  two  are  multiplied  together  and  inte¬ 
grated.  (Variations  at  the  carrier  rate  are  ignored,  as  was  done  with  the  ambiguity  function.) 

As  an  engineering  matter  the  measurement  would  probably  be  done  most  conveniently  by  passing 
the  two  received  functions  (with  one  delayed  by  Ar)  into  a  multiplying  device  and  then  measuring 
the  amplitude  and  phase  of  the  difference-frequency  sinusoid  produced  at  frequency  Af,  as  will 
be  discussed  in  the  next  section. 

The  complex  quantity  derived  in  this  way,  A(Af,  Ar),  the  two-frequency  correlation  function, 
can  be  shown  to  be  the  inverse  Fourier  transform  of  o(t,  f),  the  target  scattering  function 


A(Af,  At) 


<t(t,  f)  exp[j2ir(Afr  +  ATf)]  dT  df 


(41) 


or.  using  the  notation  introduced  in  Sec.  II-A,  A(Af,  At)  =>  ar(r,  f). 

Setting  Af  *  0  we  find  as  the  echo  correlation  function  A(0,  At),  which  we  abbreviate  A(At), 

«(At)  ==>  <r(f)  ,  (42) 

the  inverse  Fourier  transform  of  the  received  frequency  spectrum  a(f)  obtained  when  a  sinusoid 
is  transmitted  [Eq.  (39)],  just  as  one  would  expect.  Similarly,  setting  At  =  0  we  find  that  the 
spaced-frequency  correlation  function  R(Af,  0)  =  A(Af)  is 

«t(Af)  =*ff(T)  ,  (43) 

where  o(r)  is  given  by  Eq.  (37).  So,  referring  back  to  Fig.  9,  we  expect  the  maximum  At  value 
for  which  correlation  of  the  fluctuations  is  still  preserved  to  be  roughly  l/B,  the  reciprocal  of 
the  doppler  smear.  Also,  the  maximum  frequency  separation  Af  for  which  the  fluctuations  are 
still  correlated  is  roughly  l/L,  the  reciprocal  of  the  multipath  delay  smear. 

The  three  quantities  <t(t),  <r(f),  and  <t(t,  f)  are  all  real  and  positive  (since  they  are  powers), 
but  are  not,  in  general,  even  functions.  Therefore,  their  inverse  transforms,  A(Af),  A(At) 
and  A(Af,  Ar),  respectively,  will  be  even  functions,  but  will  not,  in  general,  be  real.  However, 
if  one  of  the  a's  happens  to  be  even  [as  would  be  the  case  for  e(f)  of  a  rotating  uniform  sphere] 
then  the  corresponding  transform  A  will  be  real. 

Figure  1 2  presents  a  summary  of  the  Interrelations  between  the  scattering  function  v(r,  f), 
the  two-frequency  correlation  function  A(Af,  At),  and  the  four  derived  quantities  <r(r),  cr(f),  A(Af) 
and  A(At).  The  method  of  measuring  each  of  these  will  be  summarized  presently  (in  Table  I). 

The  reader  may  wonder  at  this  point  why  the  scattering  function  has  been  singled  out  as  the 
basic  quantity  and  all  the  others  (including  the  two-frequency  correlation  function)  spoken  at  as 
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TWO- FREQUENCY  CORRELATION  FUNCTION 

Fig.  12.  Summary  of  the  interrelationship  of  the  largo!  scattering  function  and  rotated  quantifies.  (The 
quantities  A  can  all  bo  complex  In  general .  For  clarity  they  are  drawn  as  positive  real  quantities.) 
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Function 

flower  Impulio  response  o(t) 


TABLE  I 


Echo  power  ^ectrum  a(f) 


Scattering  function  o(r,f) 


Two-frequency  correlation 
function  A(Af,AT) 

(a  complex  quantity) 


Spaced  frequency  correlation 
function  A  (Af)  (complex) 


tie  ll.  ||  .1  *  * - - A. 

ISfmoQ  pf  MWfyfyftfnr 

Average  the  power  vt  range  (r)  of  a  large  number 
of  short  poises. 

Spectral  analysts  of  a  long  record  of  the  signal 
received  when  a  sinusoid  Is  transmitted. 

Tran  an  it  repeatedly  a  waveform  x(t)  with  duration 
T  equal  to  the  reciprocal  of  the  desired  frequency 
resolution,  and  with  a  bandwidth  W  equal  to  the 
reciprocal  of  the  desired  range  resolution.  At  the 
receiver,  either  (1 )  build  a  filter  matched  le  x(t) 
and  observe  average  envelope  squared  of  the  out¬ 
put  at  desired  values  of  delay  and  frequency  offsets 
f  and  t,  or  (2)  use  the  correlation  detection  equiv¬ 
alent;  i«.,  multiply  the  Incoming  signal  by  a  rep¬ 
lica  of  the  Irammlttian  x(t),  except  that  the  replica 
Is  offset  by  t  In  delay  and  f  +  fj  in  frequency.  The 
product  Is  integrated  (in  a  bandpass  filter  tuned  to 
fj)  over  the  arriving  signal  duration  (“T)»  and  the 
mean  square  of  the  envelope  of  all  such  integrations 
Is  determined.  Rapaat  far  other  t  and  f  settings 
over  the  range  af  t  and  f  values  covered  by  the 
target. 


Transmit  two  tine  waves  of  known  spacing  Af .  Arrange 
to  receive  the  two  sinusoids  separately;  insert  a  de¬ 
lay  At  In  one,  and  measure  average  amplitude  and 
phase  of  the  difference  frequency  tone  (at  frequency 
Af)  generated  by  passing  the  two  into  a  mtxer(mul- 
tiplier).  lepeat  for  other  values  of  Af  and  At. 

(Those  is  measured  relative  to  that  which  would  result 
from  a  point  target  observed  with  the  Af  in  use.)  The 
averaging  af  phase  and  amplitude  of  the  difference- 
frequency  tone  is  done  by  patting  the  mixer  output 
into  a  narrow-band  filter  tuned  to  Af. 

Same  at  above  except  use  zero  relative  delay  At. 


Transmit  a  single  sinusoid.  Past  received  signal  into 
one  Input  of  mixer  and  alto  into  the  other  mixer 
input,  but  with  a  delay  of  At  and  a  frequency  off¬ 
set  of  tome  ffj-  Examine  average  amplitude  and 
phase  of  mixer  output  at  fj. 


Echo  correlation  function  A(At) 
(complex) 


being  derivative  from  it.  Measurement  of  the  entire  target  scattering  function  o(t,  f)  itself  [or 
equivalently,  A(Af,  At)]  may  not  be  very  interesting  in  some  circumstances  (notably  where  the 
target  is  uniform).  However,  by  starting  our  presentation  with  this  basic  quantity  and  then  dis¬ 
cussing  the  five  derived  functions,  a  more  unified  presentation  has  been  possible.  And  more 
important,  as  we  have  seen  from  our  earlier  discussion  of  hard  targets,  features  in  the  scattering 
function  a  may  be  associated  directly  with  physical  features  on  the  object;  this  direct  connection 
is  not  available  with  the  correlation  function  A. 

However,  in  certain  circumstances,  the  two-frequency  method  leading  to  A  might  be  pre¬ 
ferable;  the  scattering  function  can  then  be  obtained  by  Fourier  transformation.  The  choice  of 
method  of  observation  depends  on  the  characteristics  of  the  transmitter  (e.g.,  how  it  is  limited 
in  power  and  amplitude  form  factor  of  the  exciting  signal)  and  receiver  signal  processing  equip¬ 
ment  (e.g.,  how  many  variables  of  what  fluctuation  rate  can  be  operated  upon  simultaneously). 

F.  Some  Equipment  Problems 

It  is  now  time  to  deal  in  more  explicit  detail  with  methods  of  implementing  some  of  the  op¬ 
erations  we  have  discussed.  In  radar  astronomy  practice,  the  choice  between  matched  filtering 
and  the  equivalent  correlation  detection  is  usually  decided  in  favor  of  the  latter.  For  certain 
simple  waveforms  liVe  a  single  square  pulse  it  often  pays  to  synthesize  a  filter  having  approxi¬ 
mately  the  desired  response.30  However,  for  waveforms  with  large  time -bandwidth  products, 
only  rarely31'32  does  it  pay  to  synthesize  a  filter  with  such  a  response  in  preference  to  simply 
generating  the  waveform  and  using  it  in  a  correlation  device. 

So  at  this  point  we  will  focus  our  attention  on  the  practical  problems  of  doing  the  cross- 
correlation  operation,  with  particular  attention  to  the  implications  of  the  statement  that  the  de¬ 
sired  cross -correlation  function  is  sometimes  a  complex  quantity.  We  have  met  the  complex 
correlation  function  of  signals  at  spaced  antennas  ( Z )  and  at  spaced  frequencies  (A);  we  shall 
presently  discuss  a  similar  correlation  function  in  connection  with  polarimetric  observations. 

First  we  must  explain  the  operation  depicted  in  Fig.  11(b).  It  is  well  known  that  when  the 
target  is  a  point  nonfluctuating  target  and  carrier  phase  is  unknown,  detectability  is  maximized 
by  observing  the  envelope  of  the  (complex)  cross  correlation  of  w(t),  the  incoming  Bignal,  and 
the  transmission  x(t).  A  straight  multiplication  of  w(t)  and  x(t)  followed  by  integration  will  pro¬ 
duce  a  number  f0(tj)  of  Eq.  (26)  which  is  the  real  part  of  the  correlation,  and  so  will  the  matched 
filtering-plus-sampling  operation.  But  the  desired  envelope  |-/0(tj)|  =  |  f0(t4 )  |  is  the  square 
root  of  the  sum  of  squares  of  this  quantity  fQ(tj )  =  Re  {/Q(t^ )}  and  its  Hilbert  transform, 

Im  {/0(tj)},  and  we  have  no  way  of  observing  anything  but  the  former  for  the  schemes  described 
in  Fig.  11(a).  However,  as  the  reader  may  recall  from  a  statement  made  in  the  introduction  to 
this  report,  if  we  can  somehow  contrive  to  make  the  correlator  output  f0(tj )  a  bandpass  function, 
then  Im  {/Q(tj)}  i»  observable  as  the  quadrature  component  of  this  output.  And  this  is  easily  done 
by  frequency-offsetting  either  w(t)  or  x(t)  by  an  amount  f^,  doing  the  integration  in  a  narrow  band¬ 
pass  filter  tuned  to  f^,  instead  of  the  low-pass  filter  implied  by  the  integral  sign,  and  finally  ex¬ 
amining  the  in-phase  and  quadrature  components  of  the  f^  tone  appearing  at  filter  output  at  time 
tj.  The  in-phase  and  quadrature  components  are  Re  {^(tj)}  and  Im  {/0(t^)},  respectively,  and 

the  desired  |/  (t. )  |  is  of  course  the  envelope  of  the  tone  observed  at  t. . 

O  *  33  1 

Such  a  device  is  called  a  bandpass  correlator,  and  is  not  only  capable  of  providing  the 

correlation  function,  its  Hilbert  transform  and  its  envelope,  but  is  also  much  simpler  to  build 
than  the  straight  correlator  Indicated  at  the  right  in  Fig.  11(a),  for  which  the  integrator  is  a 
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low-pass  rather  than  bandpass  filter.  The  true  correlator  requires  an  accurate  multiplier  of 
wide  dynamic  range  and  low  zero-drift,  and  the  integrator  (low-pass  filter)  must  be  built  with 
small  DC  zero  drift.  With  the  bandpass  correlator  the  multiplier  can  be  such  a  simple  element 
as  a  mixer  tube  biased  to  minimize  undesired  intermodulation  products  at  frequency  fj,  and  the 
integrator  can  be  any  form  of  narrow-band  element,  for  example,  a  quartz  crystal  filter. 

As  for  the  matched  filter  equivalent  shown  at  the  left  in  Fig.  11(b),  the  observation  of  the 
envelope  at  tj  is  done  by  building  the  filter  at  bandpass  (and  without  a  frequency  offset),  and 
observing  the  envelope  of  the  output  time  function  in  the  neighborhood  of  tj. 

The  envelope  detection  that  follows  either  the  bandpass  correlator  or  matched  filter  can  be 
performed  by  using  the  usual  combination  of  diode  and  RC  circuit  to  follow  the  crests  at  the  si¬ 
nusoidal  tone.  At  high  output  signal-to-noise  ratios  this  registers  the  envelope  exactly. 

The  correlation  function  expressions  presented  in  this  report  have  involved  infinite  time 
averages,  an  obviously  academic  restriction.  In  actual  practice,  not  only  is  the  duration  of  the 
experiment  finite,  but  the  phenomenon  being  observed  may  not  persist  for  very  long.  Moreover, 
lack  of  complete  equipment  phase  stability  over  a  long  period  may  make  it  difficult  to  separate 
the  real  and  imaginary  parts  of  the  correlator  output  by  in-phase  and  quadrature -phase  detection 
over  more  than  just  a  local  portion  of  the  experiment.  ThAt  is,  the  integrating  filter  can  be  no 
narrower  than  the  reciprocal  of  the  stability  time  of  the  conversion  equipment. 

The  effect  of  noninfinite  integration  time  is  to  introduce  noise  of  two  kinds  into  the  correlator 
output.  The  additive  noise  that  is  picked  up  by  the  receiving  antenna  and  also  generated  in  the 
receiver  front  end  will  not  be  integrated  out  completely.  Also,  since  the  processes  being  cross- 
correlated  are  noise-like  fluctuations  (usually  gaussian),  the  output  resulting  from  a  noninfinite 
integration  time  will  not  converge  to  the  true  value.  The  fluctuation  about  the  true  value  due  to 
finite  integration  time  is  classed  as  self-noise,  a  term  describing  any  fluctuation  components  in 
the  output  that  do  not  go  to  zero  as  the  additive  noise  is  removed.  The  effect  of  finite-time  in¬ 
tegration  in  bandpass  correlators  using  filters  of  arbitrary  shape  has  been  discussed  in  Ref.  33. 

Often  the  integration  time  provided  by  the  methodology  just  outlined  is  not  enough  to  produce 
an  adequately  high  output  signal-to-noise  ratio.  One  must  then  resort  to  postdetection  integra¬ 
tion;  i.  e.,  integration  of  the  output  of  the  bandpass  correlator  after  envelope  detection  (which  of 
course  destroys  the  resolution  into  in-phase  and  quadrature  components).  Figure  11(c)  shows 
how  in  the  weighted  radiometer  the  detectability  of  a  point  fluctuating  target  is  maximized  by  in 
effect  following  a  bandpass  correlator  (with  integration  time  equal  to  the  fluctuation  period  l/B 
of  the  target)  with  a  postdetection  integration  of  the  squared  envelope. 

Often  it  is  desirable  to  use  a  linear  envelope  detector  rather  than  a  square-law  envelope 
detector  in  systems  like  the  weighted  radiometer.  The  only  effect  of  this  is  a  small  reduction 
of  about  0.2db  in  the  postdetection  output  signal-to-noise  ratio.34 

Another  artifice  that  is  sometimes  resorted  to  in  practice,  and  one  that  has  more  serious 
consequences,  is  the  cross  correlation  of  the  envelopes  of  two  signals  instead  of  the  signals 
themselves.  This  has  been  done,  for  example,  with  interferometry  because  of  the  difficulty  of 
transmitting  the  two  signals  to  the  multiplier  point  with  the  phase  relations  still  preserved.  If 
the  two  signals  are  gaussian  and  S  (r)  is  the  complex  cross-correlation  function  (normalized  to 
the  product  of  the  two  powers),  then  the  (real)  cross-correlation  function  of  the  two  envelopes 
^ENV  i8  a  funotion  °*  envelope  of  S  only35 
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*>ENV<r)  =  Const  x  [2E(|-?|)-  (i  -|S|2)  K(|5|)]  =  Const  x  [|5|2  +  +  l^6  +. .  .J  , 

where  E  end  K  ere  elliptic  integrals.  Since  1 <1,  <?KNV  elweys  differs  from  |<?(r)|2  by  less 
then  6  percent.  However,  only  the  megnitude  of  the  complex  correlation  function  is  preserved, 
end  elso  there  is  cleerly  a  serious  loss  in  signal-to-noise  ratio  whenever  such  a  nonlinear  oper¬ 
ation  as  envelope  detection  is  carried  out  on  a  weak  signal  mixed  with  noise.  Specifically,  if  the 
signal-to-noise  ratio  p  of  one  member  of  the  pair  is  less  than  unity,  then  after  the  envelope  de¬ 
tection  the  signal-to-noise  ratio  of  that  member  will  be  p2.  This  is  the  well-known  weak-signal- 
suppression  effect  (see,  e.g.,  Davenport  and  Root3*1). 

The  methods  of  instrumenting  measurements  of  <r(r,  f)  and  related  quantities  are  summarized 
in  Table  I,  which  supplements  Fig.  1 2  by  giving  a  short  recipe  for  the  measurement  of  the  target 
scattering  function  c(t,  f)  and  the  five  derived  quantities  ®(Af,  At),  ®(Af),  <R(Ar),  cr(f)  and  <r(r). 

G.  Interferometry 

In  Sec.  H-D  it  was  pointed  out  that  by  computing  Z  (cf,  At),  the  cross  correlation  of  spaced 
antenna  outputs,  it  was  possible  to  determine  the  strip-integrated  brightness  distribution,  or 
what  is  equivalent  to  it  for  a  uniform  sphere,  the  echo  power  spectrum.  Moreover,  it  was  shown 
that  measurements  on  ground  pattern  velocity  could  be  made  by  determining  the  value  of  Ar  giving 
maximum  cross  correlation  for  a  given  spacing  cf.  The  ground  plane  was  assumed  to  be  perpen¬ 
dicular  to  the  line  of  sight. 

We  now  touch  briefly  on  some  of  the  methodology  involved  in  using  radio  interferometers  to 

study  the  received  signal.  The  value  of  radar  interferometric  measurements  for  radar  astronomy 
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was  first  pointed  out  by  Manasse,  who  suggested  interferometric  methods  at  studying  surface 
roughness  and  rotation  vector.  Here  we  will  mention  a  few  of  the  characteristics  of  interferom¬ 
eters  and  defer  until  Sec.  V-D  and  V-F  a  discussion  of  the  specific  experiments  proposed  by 
Manasse. 

In  actually  carrying  out  interferometer  measurements  the  stations  will  not,  in  general,  lie 
on  the  plane  perpendicular  to  the  line  of  sight;  therefore,  an  obvious  projective  correction  is 
necessary.  As  long  as  the  station  separation  as  projected  along  the  line  of  sight  is  not  an  appre¬ 
ciable  fraction  of  the  wavelength  in  space  of  the  modulation  bandwidth,  the  results  of  Sec.  Q-D 
still  apply.  (That  is,  the  attempted  range  resolution  of  the  modulation  must  not  exceed  the  range 
difference  at  the  two  stations.)  This  condition  will  be  fulfilled  until  radar  range  resolution  of, 
say,  several  tens  of  miles  or  less  is  attempted.  Let  us  then  refer  to  cf  as  the  actual  physical 
baseline  as  projected  on  this  perpendicular  plane. 

The  fact  that  the  baseline  orientation  changes  as  the  earth  rotates  can  be  used  to  advantage. 

In  Sec.  U-D  it  was  pointed  out  that  the  behavior  of  the  correlation  Z  is  interestingly  different 
when  the  baseline  <T  is  chosen  to  lie  in  different  directions.  In  practice,  provided  opportunities 
to  observe  are  sufficiently  frequent,  it  may  not  be  necessary  to  build  more  than  a  two-station 
interferometer  since  different  directions  and  magnitudes  of  <f  become  available  as  the  earth 
presents  different  aspects. 

The  directivity  pattern  of  any  antenna  system  is,  of  course,  the  Fourier  transform  cf  its 
aperture  distribution  (as  we  have  seen  in  regarding  the  distant  target  as  a  transmitting  antenna). 
By  reciprocity,  the  same  thing  holds  true  for  the  interferometer,  as  for  any  receiving  antenna 
the  outputs  of  whose  elements  are  combined  coherently.  Thus,  it  is  easy  to  visualize  why  the 
interferometer  output  |  Z  |  drops  as  the  spacing  in  wavelengths  Increases  by  more  then  R/D, 
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where  R  is  the  target  range  and  D  is  the  diameter  of  the  central  patch  at  the  brightness  distri¬ 
bution  in  which  most  of  the  reflected  power  is  concentrated.  The  power  will  drop  as  the  width 
of  the  antenna  lobe  becomes  narrower  than  D  because  of  phase  cancellation  (contributions  from 
alternate  lobes  combine  in  antiphase). 

Another  fact  that  can  be  easily  visualized  using  this  Fourier  transform  way  of  looking  at  the 
directivity  pattern  is  that  for  planetary  radar  astronomy  experiments  we  can  neglect  the  finite 
size  of  the  individual  antennas  and  regard  them  as  point  receivers,  a  condition  we  have  been  as¬ 
suming  throughout  Secs.  II-A  and  II-D.  The  planets  subtend  angles  of  about  one  minute  of  arc  at 
most,  and  until  the  individual  elements  making  up  the  interferometer  attain  this  order  of  beam- 
width  the  nonzero  size  can  be  neglected.  The  sun  and  moon,  with  their  one-half  degree  angular 
widths,  do  require  antenna  widths  to  be  taken  into  account. 

A  unified  discussion,  using  the  Fourier  transform  picture,  of  the  effects  of  nonzero  element 
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size  is  given  in  papers  by  Bracewell,  '  whose  treatment  also  allows  for  such  variations  in 
realization  as  multielement  interferometers  and  "crosses,"  phase -switching  techniques,  and  the 
fact  that  one  may  not  be  measuring  Z  (both  in-phase  and  quadrature  components)  but  only  ft t(Z). 

In  practice,  the  complex  quantity  Z  would  appear  as  a  difference-frequency  sinusoid  (of 
frequency  f^)  generated  by  offsetting  the  local  oscillator  injections  of  one  of  the  two  receivers  fay 
f .;  the  in-phase  and  quadrature  components  of  this  sine  wave  then  correspond  to  Re  (Z)  and 

U  in 

im  (Z).  Gehrels  and  Parsons  have  described  the  construction  and  use  of  an  interferometer 
which  is  associated  with  the  Millstone  Hill  radar  and  which  haB  a  sufficiently  long  baseline  for 
planetary  studies.  The  necessary  phase  identity  of  the  two  local  oscillator  injections  is  obtained 
by  using  a  special  line-of-sight  radio  link  between  stations. 

H.  Polarimetry 

The  object  of  a  radar  experiment  using  polarimetry  will  usually  be  to  measure  total  Faraday 
rotation  x  “  X0  +  Ax  of  a  soft  ionized  target  [xQ  is  the  polarization  orientation  in  the  absence  of 
the  ionization  and  Ax  is  given  by  Eq.  (20)]  or  the  degree  of  polarization  m  of  the  echo  from  a  hard 
target,  or  perhaps  both  for  a  hard  target  echo  passing  through  an  ionized  region.  As  mentioned 
in  Sec.  H-E,  there  are  several  schemes  in  use  for  specifying  the  polarization  of  an  electromag¬ 
netic  signal,  for  example  the  Stokes  vector,  Eq.  (16).  For  radar  astronomy  the  problem  of  meas¬ 
uring  any  or  all  of  the  four  parameters  of  the  Stokes  vector  does  not  differ  appreciably  from  the 
radio  astronomy  case,  and  there  is  a  fairly  large  literature  available  on  such  radio  polarimetry 
instrumentation.  The  reader  is  therefore  referred  to  the  survey  paper  by  Cohen,10  which  dis¬ 
cusses  both  theory  and  technique  in  some  detail. 

For  example,  if  the  receiving  antenna  structure  consists  of  two  elements  responding  to  left- 
circular  and  right-circular  polarization,  the  four  Stokes  parameters  are  given  by 

I.<|VlI2>+<|Vr|>2  =  pl  +  pr  . 

Q  ■  2  Re  <VlVr>  *  2A  coa? 

U  >  2  Im  O'lV'iP  =  2A  siny 

V"  <|Vl|2>-<|Vr|2>-Pl-Pr  .  (44) 
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The  two  quantities  appearing  in  the  equations  for  I  and  V  are  the  average  powers  of  the  two 
antenna  outputs.  The  quantity  <  V^VR )  is  the  cross  correlation  of  the  two  outputs,  and  can  be 
obtained  in  practice  by  passing  the  two  aqtenna  outputs  into  the  two  inputs  of  a  mixer  after  first 
frequency-offsetting  one  of  the  two  by  f^;  the  output  tone  at  frequency  f^  is  then  averaged.  This 
tone  will  have  an  average  amplitude  and  phase  denoted  above  by  A  and  y,  respectively.  In  prac¬ 
tice  the  averaging  interval  would  not  be  infinite  as  implied  by  the  abbreviation  <  >,  partly  be¬ 
cause  it  may  be  desirable  to  study  the  slow  time  variations  in  the  Stokes  vector,  in  particular 
the  quantities  Q  and  U  which  place  in  evidence  the  angle  y.  If  the  degree  of  polarization  m  is 
the  only  interesting  quantity  and  one  is  not  trying  in  addition  to  follow  slow  variations  of  x  in 
Fig.  7,  then  the  averaging  interval  can  be  extended. 

Note  from  Eq.  (18)  that  the  degree  of  polarization  m  is  not  given  simply  by  (PL  -  PR)/ 

(PL  +  PR)  unless  one  can  be  sure  that  the  right-  and  left-hand  received  signals  are  uncorrelated 
so  that  Q  =  U  =  0.  In  practical  terms  this  means  that  in  order  to  determine  the  depolarization 
introduced  by  a  target,  one  must  use  a  true  polarimeter  which  preserves  the  relative  phases  of 
the  two  received  signals,  rather  than  working  only  with  the  powers  of  the  two  received  compo¬ 
nent  signals.  The  expression  (PL  -  PR)/(PL  +  PR)  gives  a  lower  bound  on  m. 

V.  MEASURING  QUANTITIES  OF  ASTRONOMICAL  INTEREST 

It  is  time  now  to  take  the  material  on  the  effects  of  target  properties  on  the  electromagnetic 
signal  (Secs.  II  and  in)  and  turn  the  question  around  to  inquire  what  astronomically  interesting 
target  properties  can  be  deduced  from  the  received  echo. 

A.  Ranging 

The  measurement  of  interplanetary  distances  and  their  various  order  time  derivatives  is, 
of  course,  a  central  astronomical  problem.  The  basic  yardstick  of  the  solar  system  is  the 
Astronomical  Unit  (a.  u.)  —  the  mean  radius  of  the  earth’s  orbit  -  whose  accuracy  has  been  sub¬ 
ject  to  steady  refinement  by  optical  means  over  a  period  of  several  centuries.  Knowledge  of  the 
position  of  bodies  in  the  solar  system  to  one  part  in  at  least  106  is  presently  possible,  so  long 
as  the  results  are  expressed  in  terms  of  the  a.u.  (At  present  10~6  represents  the  accuracy  with 
which  free-space  light  velocity  is  known.)  However,  for  some  years  various  attempted  measure¬ 
ments  of  this  distance  spread  over  a  range  of  about  a  part  in  to3.  The  extent  of  this  disagree¬ 
ment  is  evident  from  Fig.  13,  which  lists  various  suggested  values  of  the  a.  u.  The  difficulty  lies 
partly  in  the  fact  that  all  previous  nonradar  methods  of  fixing  the  size  of  the  a.  u.  depended  on 
measurements  of  angle,  and  thus,  by  virtue  of  the  short  baseline  obtainable  on  the  earth  and  be¬ 
cause  of  atmospheric  refraction,  were  likely  to  be  permanently  limited  to  about  this  order  of 
precision.  Recently,  measurement  of  the  orbit  of  a  deep  space  probe  with  a  radar  repeater 
aboard  has  been  employed**  (see  the  value  labeled  "Pioneer  V"  in  Fig.  13).  Figure  13  also 
shows  the  results  of  several  Venus  radar  experiments.  The  experiments  labeled  with  an  aster¬ 
isk  used  the  time  of  flight  to  deduce  the  a.  u. ;  those  labeled  with  a  dagger  employed  doppler  shift 
(next  section).  Recent  radar  range  measurements*2, *3'**'*5  of  the  a.u.  agree  within  several 
parts  in  10*. 

As  radar  sensitivity  improves,  more  subtle  things  are  being  sought  in  interplanetsry  distance 
measurements.  In  addition  to  providing  refinement  of  the  single  quantity  embodied  in  the  a.u., 
it  is  possible  to  use  radar  data  to  refine  the  orbital  elements  of  various  bodies*^’*7'*8  so  that 
previously  undiscernible  effects  may  be  measured,  for  example  the  effect  of  higher  order 
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perturbation*  due  to  other  planet*.  Among  the  interesting  thing*  to  look  tor  in  the  future  are  the 
relativistic  perihelion  motion  of  an  eccentric-orbit  planet  such  a*  Mars,  the  depth  below  the  cloud 
cover  of  Venus'  reflecting  surface,  and  other  quantities.  If  a  given  planet  has  an  appreciably 
dense  ionosphere,  then  the  retardation  effects  discussed  in  Sec.  III-B  will  impair  the  accurate 
measurement  of  range,  unless  one  has  some  knowledge  of  the  target  ionosphere  density. 

From  the  discussion  in  Sec.  IV-D,  it  is  clear  that  to  measure  the  range,  one  uses  a  waveform 
with  as  great  a  bandwidth  W  as  possible.  The  accuracy  of  the  measurement  depends  on  W  and 
the  signal-to-noise  ratio  as  will  be  discussed  quantitatively  in  Sec.  Vl-B. 

B.  Velocity 

Precise  measurements  of  the  velocity  cn  a  planet  can  also  serve  to  provide  a  determination 
of  the  a.  u.  and  this  technique  was  recently  employed  on  Venus.45  The  accuracy  goes  up  with  the 
ratio  of  total  relative  velocity  between  earth  and  target  and  the  frequency  resolution  of  the  sys¬ 
tem.  A  slight  disadvantage  of  this  method  ever  direct  ranging  is  that  the  required  velocity  ac¬ 
curacy  increases  whan  the  body  is  closest  to  the  earth  (since  the  total  relative  velocity  is  a  mini¬ 
mum  at  this  time)  and  yet  this  may  be  the  only  time  when  the  object  is  detectable. 

Clearly,  to  measure  the  velocity  most  accurately  one  usee  a  waveform  with  as  great  a  dura¬ 
tion  T  as  possible.  The  accuracy  of  tha  measurement  depends  on  T  and  slgnal-to-nolse  ratio, 
aa  will  be  discussed  in  Sec.  Vl-B. 

C.  Planetary  tadtua 

Radar  measurement  of  the  radius  of  a  target  planet  or  satellite  is  not  of  pressing  interest, 
since  usually  this  quantity  is  already  known  optically.  To  measure  target  radius  from  the  total 
extent  at  e(r)  or  eff)  is  unpromising,  since  the  energy  returned  from  the  limb  will  be  very  small 
because  at  e  grazing  incidence  of  the  line  of  sight  from  the  radar,  and  will  thus  usually  be 
buried  in  noise. 

It  is  worthwhile  to  recall,  however,  that  a  measurement  of  tha  target  scattering  function  of 
a  hard  target  gives  information  on  shape  and  aims.  As  Pig.  4(b)  shows,  with  proper  normalisation 
of  the  frequency  axis,  nonzero  values  of  the  scattering  function  in  the  r,  f  plane  occupy  e  region 
that  reproduces  tha  shape  of  the  target  body,  as  viewed  from  the  side.  Therefore  the  radius  of 
a  spherical  planet  can  be  determined  from  the  high  signal-to-noiee  ratio  portions  at  the  observed 
scattering  function  (independently  of  knowledge  of  rotation  vector  and  thus  doppier  width)  by  fitting 
an  ellipse  to  the  perimeter  of  observed  nonzero  values  of  received  power  In  the  r,  t  plane  and  read¬ 
ing  off  tha  length  at  the  r  semiaxis. 

In  what  follows,  we  shall  assume  that  the  planetary  radius  a  is  one  of  the  known  quantities. 

D.  Radar  Angular  Spectrum  of  a  Uniform  Hard  Target 

The  link  between  the  detailed  character  of  a  portion  of  the  rough  surface  of  the  target  end  the 
received  signal  It  tha  radar  angular  power  spectrum  aQ  at  that  portion.  As  w*  saw  in  Sac.  H-A, 
the  angular  power  spectrum  expresses  the  statistical  character  of  tha  corrugations  of  phase,  am¬ 
plitude  and  vector  direction  of  tha  electromagnetic  signal  returning  from  a  plane  region  near  to 
the  surface  of  the  target  on  which  a  plana  wave  is  incident.  One  of  the  most  useful  types  of  ob¬ 
servation  of  hard  targets  to  be  made  by  radar  is  tha  determination  of  the  angular  power  spectrum 
for  different  surface  regions. 
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The  astronomical  significance  of  radar  Information  on  the  target  surface  is  obvious.  Indeed, 
it  is  interesting  to  note  how  significant  the  radar  observations  on  the  moon  have  turned  out  io  be 
considering  the  moon's  accessibility  to  optical  observation. 

If  the  spherical  surface  can  be  assumed  uniform  (i.e..  »0(<p)  the  same  everywhere),  than  de¬ 
termination  of  »o  is  straightforward.  A  narrow-pulse  experiment  to  determine  the  power  impulse 
response  e(r)  [Eq.  (37))  gives  *0(p)  directly  by  the  inverse  of  Eq.  (13). 

•0M  =  Tic  c0**’  •  («> 

In  practice,  as  a  result  of  nonzero  pulse  width  one  will  actually  observe  not  <r(r)  but  P(r)  of 
Eq.  (38). 

This  short-pulse  method  of  determining  e(r)  is  probably  the  simplest  method  at  all.  A 
mathematically  almost  equivalent  method  is  to  measure  the  correlation  between  sinusoids  spaced 
by  Af  as  a  function  of  Af.  This  gives  ®(Af),  which  then  can  be  Fourier  transformed  to  give  o(r), 
provided  one  has  been  able  to  measure  the  complex  A  and  not  simply  the  magnitude. 

Obviously  <roi<p)  can  also  be  obtained  from  a  measurement  of  the  echo  power  spectrum  <r(f) 

Eq.  (39).  since  the  two  are  uniquely  related  fay  the  integral  equation  (14)  (where  the  projected 
rotation  speed  Q  cos  a  is  known).  Determination  of  e(f)  is  accomplished  by  determining  P(f) 

(Eq.  (40))  from  the  received  version  of  a  long  pulse  of  sine  wave  either  fay  direct  spectral  analy¬ 
sis  or  by  determination  of  <*(At),  the  echo  correlation  function,  followed  by  Fourier  transforma¬ 
tion.  Talcing  this  data  on  off)  is  simple  enough,  but  the  inversion  of  (14)  is  more  tedious  than 
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solving  Eq.  (45)  for  the  short-pulse  method.  The  inversion  to  Eq.  (14)  is 


It  in  worth  pointing  out  that  if  the  target  is  a  uniform  sphere,  all  four  at  these  measure¬ 
ments  (*(r),  a  if),  A( Af)  and  A  (At))  are  mutually  redundant.  That  is,  any  oaa  of  them  suffices 
(if  a  and  0  cos  a  are  known)  to  uniquely  specify  tbs  scattering  function  r(r,  f)  or,  equivalently, 
the  correlation  function  t(Af,  At).  This  is  clear  from  the  development  of  Appendix  B  in  which 
it  is  shown  that  the  angular  power  spectrum,  target  radius  and  projected  rotation  speed  com¬ 
pletely  define  the  target  scattering  function  if  the  sphere  is  uniform. 

A  fifth  method  of  measuring  »#(f)  is  the  interferometer  method37  of  correlating  the  output 
of  spaced  receivers  separated  fay  distance  7  and  with  relative  delay  At  ■  0  to  give  Z(  J,  0)  as  in 
See.  n-D  and  Fig.  6(b).  If  the  target  is  uniform,  then  the  brightness  distribution  P(x,  y)  of 
Fig.  3(b).  will  be  circularly  symmetric  and  Z  becomes  independent  at  the  orientation  of  if.  That 
is,  Z(«f,  0)  ■  Zf d,  0).  Repeating  the  measurement  at  a  number  of  spacings  to  get  Z  (d,  0)  will  give 
the  brightness  distribution  1 1 1  ■  P(x,  y)  of  Eq.  (7)  and  thence  e0(y)  from  Eq.  (•).  However,  one 
does  not  need  to  go  to  all  the  trouble  of  using  variable  spacing,  because  the  variable  Ar  is  avail¬ 
able.  The  echo  power  spectrum  «r(f)  at  a  single  station  is  the  inverse  Fourier  transform  of  the 
echo  correlation  function  *(At),  which  is  recognised  to  be  none  other  than  Z( 0,  Ar).  Thus, 
Manasse  suggests  that  measurement  of  the  power  spectrum  of  the  signal  received  at  a  single 
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station  gives  the  desired  radar  angular  spectrum  aQ{+)  through  Eq.  (46)  unless  perchance  one 
does  not  know  the  projected  rotation  rate  0  cos  a,  which  is  a  scale  factor  on  the  bandwidth  of 
ait),  in  which  case  determination  of  Z  <d,  At)  for  a  single  value  of  d  resolves  this  scale  factor. 
That  Is,  both  ?0(p)  and  0  cos  a  can  be  determined  by  a  cross  correlation  between  outputs  of  two 
stations  plus  an  autocorrelation  at  one  of  them.  The  experiment  gives,  then,  information  about 
surface  roughness  and  rotation  vector.  We  shall  say  more  about  measurement  of  the  rotation 
vector  3  in  Sec.  V-F. 

The  surface  roughness  of  a  hard  target  can  also  be  inferred  grossly  from  measurement  of 
the  degree  of  polarisation  m  of  the  returning  echo  signal,  as  was  mentioned  in  Sec.n-E.  This 
type  of  measurement  has  been  made  on  Venus.50 

In  case  the  body  can  be  regarded  aa  composed  partly  of  an  Ideally  smooth  and  specularly 
reflecting  material  and  partly  of  rough  material,  the  amplitude  probability  density  distribution 
affords  a  crude  method  of  deducing  the  ratio  of  powers  due  to  the  two  components.5*  For  a  com¬ 
pletely  specular  smooth  surface,  the  received  signal  would  net  fade  at  eB  but  would  maintain  a 
constant  value.  A  rough  surface  will  produce  a  gauss lan  distribution  of  voltages  about  a  aero 
mean  (i.e.,  a  Rayleigh  distribution  of  the  envelope).  In-between  cases  give  an  intermediate  den¬ 
sity  distribution.52  This  first-order  statistic  Is  relatively  useless  compared  to  the  second -order 
quantities  we  have  been  discussing,  since  it  is  relatively  insensitive  to  the  ratio  of  specular  and 
scattered  power  and  thin  ratio  would  thus  be  quite  difficult  to  measure  accurately  with  noise 
present. 

E.  Isolating  Individual  Surface  Features 

A  more  realistic  model  of  the  surfaces  of  planets  and  the  moon  la  achieved  fay  assuming  that 
they  are  not  uniform  in  surface  properties.  Ideally  one  would  like  to  resolve  individual  surface 
regions  and,  by  studying  them  at  different  aspect  angles  p,  learn  something  about  their  specific 
nature  through  determination  of  their  angular  spectra  oQ(p);  thence,  using  the  theories  discussed 
in  Ref.  3,  one  could  decide  on  reasonable  physical  models  for  the  terrain  fantures.  One  might 
also  study  the  polarization  properties  and  the  amplitude  probability  distribution  of  returns  from 
isolated  regions. 

The  moat  obvious  approach  is  to  point  a  sufficiently  narrow  antenna  beam  at  the  region  in 
question,  but  this  requires  such  enormous  antenna  gains  as  to  appear  unpromising  at  the  moment, 
except  in  the  case  of  the  moon  or  sun,  where  it  is  only  now  becoming  possible  to  resolve  to  sev¬ 
eral  tenths  of  a  target  diameter. 

Interferometry  next  suggests  itself,  but  if  the  interferometer  consists  of  only  a  pair  of  an¬ 
tennas  with  fixed  spacing  to  give  the  desired  resolution,  the  echoes  picked  up  in  the  undesirod 
lobes  render  the  received  data  too  noisy  to  isolate  a  given  target  feature.  Two  measures  are 
then  possible:  (1)  the  original  element  spacing  d  is  filled  with  n  - 1  other  equlspaced  antennas, 
whereupon  there  is  only  one  major  lobe  every  nk/2d  radians  instead  of  every  k/2d,  as  would  be 
the  case  with  only  two  elements;  or  (2)  many  interferometer  records  are  made,  using  many  dif¬ 
ferent  spacinge,  so  that  the  results  may  be  combined  to  give  a  Fourier  synthesis  of  tbs  distribu¬ 
tion  of  power  over  the  target  ("aperture  synthesis”).  The  two  schemes  are  roughly  equivalent. 

However,  it  la  really  not  necessary  to  go  to  all  this  trouhle,  at  least  for  hard  targets,  be¬ 
cause  the  variables  time  delay  and  doppler  offset  provide  isolation  of  various  surface  regions, 
as  was  discussed  in  Seca.U-C  and  Q-D.  This  mapping  property  gives  a  simple  geometrical 
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correspondence  between  points  on  the  target  and  regions  In  the  r,  f  plane  of  the  target  scattering 
function  <r(r,  f).  That  is,  energy  returning  with  a  certain  doppler  shift  and  time  delay  imposed 
on  it  could  only  have  come  from  certain  places  on  the  target  surface.  To  isolate  the  pair  of 
points  heavily  shaded  on  the  surface  of  Fig.  4(a),  one  usee  a  transmitted  waveform  x(t)  of  dura¬ 
tion  T  and  bandwidth  W  such  that  the  ambiguity  function  42  has  the  desired  frequency  width  1/T 
and  range  width  l/W;  observation  of  the  power  received  at  the  appropriate  value  of  delay  r  and 
with  the  proper  doppler  offset  f  then  ensures  the  desired  isolation.^  In  order  to  be  able  to  aver¬ 
age  a  number  of  readings  of  power  in  the  r-f  cell  l/T  X  i/w  in  sine,  the  transmission  should  be 
repeated  many  times  (see  Sec.  VT-D). 

It  remains  to  resolve  the  twofold  nr  ibigutty  between  the  pair  of  points  in  Fig.  4(a).  This  can 
be  done  by  repeating  the  measurement  at  a  time  when  die  apparent  rotation  axis  orients  differ¬ 
ently  relative  to  the  observer.  Or,  spaced  receivers  will  allow  an  interferometer  null  to  be  laid 
onto  one  of  the  two  regions  unless,  fortuitously,  the  baseline  lies  perpendicular  to  the  ground 
projection  of  the  polar  axis.* 

The  range -doppler  mapping  technique  has  been  used  by  Pettengill  in  investigating  the  char¬ 
acter  of  die  lunar  surface,8'^  and  a  preliminary  such  observation  on  Venua  has  been  reported 
by  Smith.88 

So  far  we  have  discussed  isolation  of  a  fixed  region  on  die  body  by  that  part  of  the 

return  arriving  in  the  neighborhood  of  fixed  values  of  range  r  and  doppler  f .  Actually,  as  the 
body  rotates,  the  return  associated  with  a  given  target  feature  lies  at  different  places  in  the 
range-doppler  plane  aa  time  progresses.  That  la,  there  will  be  e  characteristic  trajectory  or 
range-doppler  history  through  the  r-f  plane  which  will  be  followed  fay  die  echoes  from  each  point 
on  the  target  aa  time  progresses.  In  the  spherical  case,  these  trajectories  are  unambiguous 
with  one  another,  unless  the  polar  axis  happens  to  lie  perpendicular  to  the  radar  line  of  sight, 
whereupon  there  is  still  the  twofold  ambiguity  just  referred  to.  In  Fig.  14  are  sketched  the  tra¬ 
jectories  of  several  points  on  a  typical  spherical  surface. 

Knowing  the  range-doppler  history  of  a  given  point  on  the  target,  one  can  process  returns 
taken  at  different  times  during  a  rotation  period  and  build  up  at  least  a  partial  picture  of  <r0(p), 
the  radar  reflectivity  a a  a  function  of  Incidence  angle.  It  should  be  possible  also  to  use  this 
correspondence  to  build  up  the  slgnal-to-noiae  ratio  and  also  the  degree  of  range-doppler  reso¬ 
lution  (l.e.,  the  fineness  of  the  "resolution  cells"  of  the  range-doppler  map).  Roughly  speaking, 
both  these  quantities  increase  with  increased  observing  time  on  a  given  return.  An  ingenious 
ground-mapping  radar  using  just  these  range-doppler  histories  has  been  built  fay  Cutrona,  at  tl.53 
Because  there  is  one  trajectory  for  each  of  many  thousands  of  observed  target  points,  techniques 
of  parallel  processing  of  thousands  of  hypotheses  were  necessary  in  this  equipment.  Two- 
dimensional  optical  methods  were  chosen  and  developed  to  s  high  degree.  This  sort  of  parallel 
processing  warrants  much  further  investigation  for  radar  astronomical  mapping  and  also  for 
soma  of  the  other  less  elaborate  signal  processing  procedures  ws  have  discussed. 

The  possibility  that  the  echo  in  one  range-doppler  cell  will  be  correlated  with  that  in  another 
cell  after  a  suitable  time  delay  (and  possibly  a  frequency  shift)  means  that  certain  previous 


t  Provided  the  target  It  net  ovanpraod,  at  will  be  dltcuned  In  Sec.VI-D. 

t  Exprimd  In  termt  of  voltage,  the  nullt  are  narrower  than  the  peoki  became  the  pattern  hot  the  character  of  the 
magnitude  of  a  tlnmotd.  Complete  tupprenlan  would  retulv  then  only  for  very  imall  resolution  celli.  But  thk  re¬ 
quire*  large  bandwidtht,  which  Itself  comet  a  fl!llng-ln  of  the  null,  unlett  the  center  frequency  It  wfficlently  high 
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W  <awg>  dapfikr  plana. 


Fig.  14.  Tra^actoria*  through  Hw  r-f  plan*  (ranpa-AppUr  hbtoriaa)  of  «<w  achoa*  from 
four  point*,  P ,  through  P^,  on  a  ^horioal  twrfaca. 
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assumptions  may  ultimately  have  to  be  re-examined.  As  pointed  out  in  the  footnote  preceding 
Eq.  (31),  the  scattering  function  er{r,  f)  then  no  longer  suffices  to  completely  specify  the  echo  be¬ 
havior  and  one  must  extend  the  definition  of  this  function  to  take  into  account  these  "cell-to- 
cell  coherences."  Measurement  techniques  of  astronomical  interest  that  take  these  coherences 
into  account  (as  does  the  ground -mapping  radar  just  referred  to)  await  a  systematic  investigation. 
Likewise,  not  much  is  known  about  the  effect  on  signal  detectability  of  cell-to-cell  coherences, 
although  Kailath54  has  given  the  problem  some  attention.  It  is  clear  that  known  systematic  co¬ 
herences  in  the  echo  signal  should  make  possible  an  increase  in  detectability. 

F.  station  Vector 

The  quantity  fl,  indicating  the  rotation  velocity  in  radianB  per  second,  has  appeared  in 
several  of  the  expressions  developed  in  Sec.  IV  of  this  report.  We  shall  now  inquire  how  one 
might  measure  '.he  speed  R  and  the  axial  direction.  In  the  case  of  several  planets,  notably  Venus, 
fl  is  at  present  largely  unknown. 

The  width  of  thedoppler  spectrum  constitutes  a  measure  of  the  projected  rotation  rate  0  cos  a. 
where  a  is  the  tilt  of  the  axis  from  perpendicularity  with  the  line  of  sight.  The  total  over¬ 
all  width  of  the  echo  power  spectrum  cr(f)  gives  fl  cos  a  immediately,  since  this  width  is  equal  to 
(4a/cf  )  fl  cosa.  (a  is  the  target  radius,  fQ  the  carrier  frequency.  and  c  is  propagation  veloc¬ 
ity.)  However,  the  outer  limits  of  the  echo  spectrum  are  usually  buried  in  noise  and  it  is  nec¬ 
essary  to  infer  the  spectrum  width  from  the  width  of  the  central  high-energy  portions.  But  if 
this  is  so,  one  clearly  requires  more  information  than  data  on  the  spectrum  <r(f)  alone,  since  a 
given  width  of  the  center  of  the  spectral  peak  could  have  been  caused  by  either  a  smooth  target 

rotating  rapidly  or  a  rougher  target  rotating  more  slowly. 

Information  on  surface  roughness  is  clearly  necessary  and  there  are  several  ways  of  getting 
it.  Manasse37  suggested  that  a  single  reading  of  Z(cf,  0)  (the  correlation  of  spaced  receiver  sig¬ 
nals  as  a  function  of  spacing  at  zero  delay)  would  suffice,  since  for  a  uniform  spherical  target, 
the  shape  of  Z  is  determined  once  cr(f)  Is  known,  except  for  a  scale  factor  on  width  (sue  Sec.  II-D 
and  Fig.  6).  A  simpler  method55  that  does  not  involve  a  separate  receiving  facility  is  to  deter¬ 
mine  the  roughness,  i.e.,  the  angular  power  spectrum  oQ(ip),  from  the  power  impulse  response 
ct(t)  by  Eq.  (45);  Eq.  (14)  then  gives  the  value  of  the  projected  rotation  rate  fl  cosa  that  corre¬ 
sponds  to  the  observed  and  the  observed  doppler  spectrum  o(f). 

If  the  target  is  nonuniform,  the  magnitude  fl  of  the  rotation  velocity  can  also  be  deduced 
from  periodicities  observed  in  long-term  records  of  echo  power. 

Cohen56  has  pointed  out  that  one  should  observe  a  systematic  change  of  the  width  of  the  echo 
spectrum  tr(f)  as  observed  over  many  months  because  the  rotation  vector  if  remains  fixed  in 
inertial  space  as  the  aspect  from  which  the  planet  is  viewed  changes.  Cohen  shows  how,  by  ob¬ 
serving  the  "phase  angle"  and  "amplitude"  of  this  quasi-periodic  oscillation  of  spectrum  width, 
it  is  possible  to  infer  not  only  the  rotation  speed  but  also  to  establish  the  tilt  of  the  axis.  The 
ambiguity  in  axis  direction  that  remains  is  that  the  direction  of  rotation  might  bo  clockwise  or 
counterclockwise  and  the  axis  might  be  tipped  at  some  angle  q  or  -q  to  the  plane  of  the  orbit. 

In  his  discussion  of  radar  interferometry,  Manasse37  suggested  a  way  of  resolving  the  first 
of  these  two  ambiguities,  namely  whether  it  is  the  left  or  the  right  side  of  the  apparent  disk  of 
the  target  that  is  approaching  the  observer.  One  simply  narrows  the  input  filters  of  the  two 
spaced  receivers  so  that  they  both  pass  only  the  energy  arriving  from  a  strip  of  the  apparent  disx 
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of  the  target,  detunes  both  input  filters  together  so  that  this  strip  lies  either  in  the  approaching 
side  or  the  receding  side,  and  then  steers  the  interferometer  lobe  to  maximize  received  power. 
The  direction  of  movement  of  the  best  position  of  the  lobe  can  be  compared  with  the  direction  of 
frequency  offset  to  get  the  sense  of  rotation. 

The  sense  of  rotation  can  also  be  obtained  from  interferometer  experiments  in  which  time 
delays  Ar  are  allowed,  using  the  effects  described  in  Sec.  II-D.  As  Fig.  6(b)  indicates,  for  a 
given  receiver  spacing  d,  there  is  a  value  of  Ar  that  maximizes  the  correlation  Z,  namely 
At  of  Eq.  (15).  The  true  velocity  vector  v  of  the  ground  pattern  expresses  the  project*'’ 
target  rotation  velocity  U  cos  <v,  and  can  be  found  by  making  two  observations  of  rmgx  using 
crossed  baselines,  and  solving  two  equations  like  Eq.  (15) 

As  a  practical  example  of  the  application  of  this  technique,  suppose  we  have  a  radar  trans- 
mitter  at  point  A,  and  receivers  at  points  A  and  B  (A  and  B  might  be  Millstone  Hill  and  Lex¬ 
ington,  respectively,  20km  apart).  Suppose  the  transmission  x(t)  is  a  train  of  phase-coherent 
pulses  2  msec  in  duration  occurring  every  30  msec,  and  suppose  the  planetary  target,  Venus,  is 
overhead.  We  desire  to  determine  the  ground  diffraction  pattern  velocity  as  projected  along  the 
baseline  AS  =  d.  First,  the  outputs  of  the  two  receivers  would  be  gated  into  the  two  inputs  of  a 
correlator  for  2  msec  every  30  msec  (with  the  range  of  the  planet  known)  to  give  Z(£,  0).  Then, 
to  try  different  values  of  Ar  to  find  Armax,  the  one  for  which  Z  («T,  At)  is  greatest,  one  would  de¬ 
lay  one  of  the  correlator  inputs  by  various  integral  multiples  of  30msec.  If  Venus'  rotation 
period  were  225  days  (n  ~  4  x  10"7  radians/sec),  then  the  ground  pattern  speed  would  be  some 
15  km/sec,  and  one  would  find  that  for  the  Lexington-Millstone  baseline  lying  along  the  direction 
of  pattern  motion,  maximum  correlation  would  be  obtained  for  values  of  Ar  around  1.3  seconds. 

In  summary,  it  appears  that  by  the  combined  use  of  several  types  of  radar  experiments,  it 
is  possible  to  determine  rotation  speed  (I  and  resolve  all  the  ambiguities  in  the  direction  of  the 
polar  axis  except  whether  it  is  tipped  toward  the  observer  or  away  from  him  (a  positive  or 
negative). 

G.  Ionospheres  and  Interplanetary  Ionization 

It  is  reasonable  to  suppose  that  above  certain  levels  in  planetary  atmospheres,  recombina¬ 
tion  rates  are  low  enough  and  the  incident  solar  flux  strong  enough  that  a  significant  density  of 
free  electrons  exists.  As  discussed  in  Sec.  Ill,  if  losses  are  not  excessive,  such  regions  inter¬ 
posed  between  the  earth  and  u  i ,  fleeting  hard  target  can  be  studied  in  terms  of  the  time  delay 
suffered  by  the  echo  components,  or  in  terms  of  the  critical  frequencies  at  which  the  attenuation 
changes  abruptly,  or  (if  a  magnetic  field  is  present)  in  terms  of  the  accumulated  Faraday  rotation. 

Equation  (19)  gives  an  expression  for  the  round-trip  time  of  flight  of  a  signal  reflected  from 
an  object  at  range  R.  The  second  term,  the  excess  of  time  of  flight  over  that  in  vacuo,  is  seen  to 
be  a  measure  of  ele'.tron  density  integrated  along  the  path,  and  this  quantity  appears  weighted  by 
the  inverse  square  of  the  frequency.  This  strong  frequency  dependence  suggests  several  methods 
of  measuring  integrated  electron  density.  The  most  sensitive  would  be  a  pair  of  separate  radar 
experiments  spaced  as  widely  in  frequency  as  possible,  and  with  the  lowest  frequency  as  low  as 
possible.7  The  two  measurements  of  time  of  flight  at  two  such  widely  spaced  frequencies  will  give. 


t  But  not  so  low  that  the  signal  does  not  escape  the  terrestial  Ionosphere.  Critical  frequencies  of  the  F2  layer  vary  > 

from  about  2  to  15  Mcps.  If  the  beam  goes  through  the  ionosphere  at  on  angle  x  the  zenith,  the  current  value 
of  critical  frequency  must  be  multiplied  by  sec  y  to  give  the  minimum  useful  frequency. 
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by  simultaneous  solution  of  two  equations  like  Eq.  (19),  both  the  true  range  R  and  the  integrated 

electron  density.  (It  is  assumed  that  light  velocity  c  and  freauency  f  are  known  accurately.) 

57 

A  more  convenient  procedure  is  to  send  a  modulated  signal  and  then  examine  the  relative 
delays  of  the  various  sidebands  in  the  received  signal,  but  the  limited  frequency  spread  renders 
this  method  less  sensitive  than  two  separate  experiments  at  widely  separated  frequencies.  Also, 
as  we  have  seen  in  Secs.IV-D  and  IV-E,  rapid  target  fluctuations  (doppler  spread  B)  or  long 
multipath  duration  L,  or  both,  have  their  effect  on  the  fading  of  the  sideband  components. 
Reference  57  presents  a  discussion  of  a  number  of  single -station  experiments  of  ‘his  typ**  em¬ 
ploying  various  forms  of  modulation. 

In  either  class  of  experiments,  the  fraction  of  the  integrated  electron  density  due  to  the 
interplanetary  and  target  regions  is  not  known  until  the  terrestrial  component  is  subtracted  out. 
The  interplanetary  component  is  likely  to  be  much  smaller  than  the  terrestrial  component.  Usu¬ 
ally  the  terrestrial  component  is  estimated  from  HF  ionograms,  but  this  gives  only  an  approxi¬ 
mate  figure  because  of  the  variations  of  electron  density  with  time  and  place.  Briscoe58  has  made 
the  ingenious  suggestion  that  during  a  planetary  radar  experiment,  or  immediately  before  it  or 
after  it,  one  can  examine  the  Faraday  rotation  of  incoherent  electron  scatter  returns  from  a  re¬ 
gion  well  above  the  F2  maximum  of  our  ionosphere  exactly  along  the  planetary  radar  line  of  sight. 
Knowledge  of  the  earth's  magnetic  field  then  gives  an  accurate  figure  for  the  terrestrial  compo¬ 
nent  by  using  Eq.  (20). 

If  radar  equipment  becomes  available  that  will  do  accurate  ranging  over  a  wide  tunable  fre¬ 
quency  range,  planetary  ionospheric  densities  can  be  studied  by  noting  the  critical  frequency  be¬ 
low  which  the  reflection  is  from  the  ionosphere  and  above  which  the  reflection  is  from  underlying 
regions.  The  transition  will  evidence  itself  most  easily  by  a  discontinuity  in  range  or  in  other 
parameters,  such  as  effective  surface  roughness. 

H.  Magnetic  Field 

Faraday  rotation  has  been  used  to  deduce  the  total  electron  content  of  the  earth-moon  region 
1 8 

using  Eq.  (20).  Similarly,  using  a  radar  polarimeter  for  reception  of  a  planetary  echo,  and  sub¬ 
tracting  this  same  near-the-earth  component,  it  is  possible  to  find  the  integrated  field-electron 
density  product  Rg(a)  N(s).  Then,  if  the  integrated  N(s)  has  been  obtained  by  some  alternative 
means  (see  the  preceding  section),  one  can  say  something  about  the  magnetic  field  between  the 
vicinity  of  the  earth  and  the  target  Unfortunately,  it  is  the  product  that  appears  in  Eq.  (20),  so 
it  is  not  possible  to  localize  the  field  contributions  even  if  the  location  of  the  electi  on  density 
contributions  is  known,  unless  collateral  information  of  some  sort  is  available.  The  depolariza¬ 
tion  (Sec.  II-E)  due  to  target  surface  roughness  may  impede  the  measurement  of  total  Faraday 
rotation  <p. 

Experiments  with  tunable  radars  to  discover  critical  frequencies  can  provide  very  interest¬ 
ing  information  on  the  magnetic  field.  Specifically,  if  the  ordinary  and  extraordinary  modes 
can  be  separated  either  according  to  direction  of  circular  polarization  or  propagation  delay 
imes.t  then  the  critical  frequencies  for  the  two  modes  can  be  isolated.  One  can  then  set  crude 
bounds  on  the  value  of  the  local  gyro  frequency,  which  is  proportional  to  the  magnetic  field  in¬ 
tensity. 


t  The  indices  of  refraction  are  slightly  different,  a;  discussed  in  Ref.  19. 
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A.  ftaUstteal  Bsttofti  b 

The  final  sections  of  t  •' \  report  will  bo  devoted  to  various  matter*  relating  to  errora  in 
making  the  meaaaremeut  d  curbed  earlier.  In  planning  and  carrying  out  physical  Invest lgetiona, 
it  is  important  not  only  tc  j 1  <I«  to  conceive  the  sort  of  •xptrtmint  that  will  place  in  evidence 

the  variables  to  be  mea*  ec  but  also  to  knew  whet  the  sources  of  error  are,  the  expected  size 
of  the  errors,  and  whetk*.  the  xperhnant  conceived  is  apt  to  minimize  them,  relative  to  other 
possible  alternative  ways  of  conducting  the  experiment,  hi  tee  radar  astr  jnomy  experiments  we 
have  discussed,  the  noise  added  into  the  echo  signal  is  the  principal  factor  influencing  measure¬ 
ment  error.  O  the  target  itself  produces  fluctuations  in  the  echo,  this  too  is  a  source  of  error. 

The  abstract  methodology  of  suiting  up  expe  riments  to  minimise  error  n  the  presence  of 
random  influences,  the  enleulatfen  of  what  sorts  of  error  behavior  saw  »«hir  rl  1  and  eo  forth  - 
these  are  tea  domain  of  mathematical  *t*tistic«  59  Some  at  the  particular  theories  developed  by 
ttstisticane  hash  proved  te  be  extremely  useful  to  the  communications  and  **d»r  engineer.  An 
octets  ntllng  example  is  the  theory  cf  bwethesis  testing,  which  treats  the  rater  detection  ope  ra¬ 
ttan,  far  example,  as  a  teat  between  these  two  mutually  exclusive  hypotheses;  "the  received  slg- 
aal  is  noise  plus  a  target  echo,"  and  “the  respired  signal  is  antes  alone."  Ik  the  catppgpg0n  de¬ 
port.1  this  approach  la  applied  to  the  dess  at  targets  of  interest  to  the  radar  aetrosumrf  (those 
that  are  spread  in  r  or  f  or  both)  ta  such  a  way  as  to  dertan  tee  item  tf  ttataflipltted  signal  and 
rtiettvsr  praoaasing  eqtfpaaeat  that  will  mtntoniaa  the  prctehtlftfru  at  fatijMf  to  detect  the  pres¬ 
ence  at  the  ohfect  and  of  falsely  concluding  that  it  is  present.  The  the r><p,»too.ttf>|  how  good  the 
system  will  be,  teat  is,  to  what  love*,  three  probabilities  will  be  rednoed'lit  thegpimi  system  in 
use. 

There  is  an  analogous  atatteUoal  theory,  estimation  teeaey  (see,  for  example.  H.ilstrotn,60), 
which’  te  sgpUoahto  to  super  tarn  ntti  that  ara  intended  to  measure  something,  for  exarqg’e,  a  radar 
target  parameter.  This  theory,  too,  has  proved  useful  to  the  radar  experfawater  since,  in  those 
eases  wteeb  are  umthenmtieaUy  tractable,  it  provides  him  with  the  form  tf  transmitted  signal 
and  receiver  processing  equipment  that  will  aatehnlse  the  error  with  which  fas  desired  quantity 
is  registered.  The  theory  goes  on  to  say  mv  good  the  meaenramank  is,  fate  ta,  ptat  the  average 
■pried  (variance)  in  observation  values  -.•/  :>aUI  have  been  if  tea  eqpertmata  oauld  have  has it  re- 
pasted  e  Urge  number  of  times  aaiw  su-i.rticelly  identical  oondtttans  (enaaaaUe  average).  In 
this  section  we  shall  be  interested  ia  preoeidiag  axpreeatane  far  error  vertaneeJ 

The  gild  tf  tee  statistical  estimation  approach  is  as  follows.  The  trap— Ht*f  waveform  x(i) 
m  reflected  from  the  target,  one  of  whose  properties  a  we  desire  to  measure.  The  property  <* 
might  1m  e  single  parameter  like  range  - ,  it  might  be  a  parameter  set  like  fas  pair  r  and  f.  or 
it  might  be  e  .tarnation  li'te  o(r)  or  air,  f).  Stgpcoe  tee  true  value  tf  this  gpROdhp  fa  ten  ane  de¬ 
pot -!  by  oQ.  The  reflected  signal  y(t»  arrives  at  th  1  receiver  mix-*?  with  notes  jiff);  tee  received 
0ig>*.  t  rfti  1  y(  '  4-  n(t).  C*.  A'*'  md  x)  ia  the  probability  that  the  actual  o  nerved  reoeiver 
ibp  *  a-avaTorm  wN)  did  .  ri.se  S.-om  Ihe  transpniasion  :c(t)  and  e  certain  s  lying  b.  the  Marvel 
s  t»  a  ♦  da.  (Often  ilu»  r>*  afas.mitty  can  actually  be  known  to  the  receiver  U  the  ntstteHoa  of  tee 

fitietfer  to  mdsnte  dcssitfeie  V  STe.  **  wl..  holly  be  dhtrlbeatf,  am  nett  (new  ’hetrpwfahflby 
dhkfad tan,  but  far  purposes  of  Ihe  present  evr-*y  treolm*..  w-  vH!  be  eonfant  tedtaMSoaly  die  mttaawef 
Aiew. 
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noise  n(t)  are  known.)  The  function  of  the  receiver  Is  to  produce  an  output  or  an  estimate  a* 
which  we  would  like  to  have  correspond  to  the  actual  a^,  but  which  will  not,  of  course,  because 
of  the  noise.  One  logical  procedure  for  the  receiver  to  carry  out  is  to  produce  as  the  output 
a  *  that  value  of  a  for  which  p(w/cr  and  ac)  is  a  maximum  .t  Such  a  receiver  la  called  a  Bmimimu 
nfcairti/wS  in  practice,  the  way  It  is  designed  is  first  to  determine  analytically  what 

p(w/a  and  x)  is.  Then  it  often  turns  out  that  the  expreseion  for  this  probability  can  be  interpreted 
in  terms  of  simple  signal  processing  operations  on  the  received  signal  w(t)  that  produce  an  a* 
far  which  this  probability  is  maximized,  when  the  noise  statistics  and  the  form  of  the  transmis¬ 
sion  X<t)  are  known.  After  this  is  done  it  is  usually  possible  to  go  back  and  perform  alterations 
on  the  tsniwmission  x(t)  to  still  further  improve  the  performance. 

The  maximum-likelihood  procedure  is  not  just  heuristically  plausible;  when  the  a  to  be 
estimated  is  a  parameter  or  a  parameter  set,  it  has  been  shown  that  the  a  *  based  on  n  independ¬ 
ent  such  observations  approaches  the  true  a  for  all  practical  purposes  as  n  -*  «•  (•consistent" 
and  "asymptotically  unbiased”  estimator),  and  the  varianew  of  m*  appsowehes  aero  as  that  as 
any  other  possible  form  of  estimator  ("asymptotically  efficient"  estimator).  And,  what  is  more 
important  practically,  for  a  finite  n  (e.g.,  unity),  if  the  estimator  exists  giving  the  theoretical 
■if.i-n.-i  achievable  error  variance  (an  "efficient"  estimator)  as  calculated  from  the  so-called 
Cramer-Rao  lower  bound  (on  error  variance),  then  that  method  is  a  maximum-likelihood  estima¬ 
tor.  However,  not  all  maximum -likelihood  procedure#  achieve  this  lower  limit.  Indeed,  in 
mw«  onsag,  the  Cramer-Rao  lower  bound  is  lower  than  that  achievable  by  any  realizable  estima¬ 
tor;  for  example,  sometimes  it  is  zero. 

Suffice  it  to  say  that  a  maximum-likelihood  strategy  is  the  most  appropriate  approach  under 
moat  conditions,  with  an  exception mentioned  in  the  next  paragraph.  It  is  often  mathemati¬ 
cally  tractable  and  gives  estimators  that  have  been  found  by  experience  tegs  useful  and  are  often, 
in  fact,  minimum -variance  estimators.  In  measuring  the  evens  section  of  a  point  eeaOaetuating 
target,  tor  the  minimum-variance  estimator  exists  and  is  a  iW— Imjji i -likelihood  esti- 

mt  .vr.  For  eatimat.  ng  the  scattering  function  <r(r,  f),  on  the  other  hand,  methods  dP  getting  the 
:m.  -j  mum-likelihood  mtimatee  have  recently  been  found,  but  at  present  it  is  not  known  whether 
f,--  variance  of  the  ertlmate  is  as  small  as  is  theoretically  achievable.  For  measurements  of 
,-ivre  involved  quantities,  estimation  theory  has  been  little  need  up  to  now.  For  example,  no 

las  worked  out  a  ^■»<m,«m-jik«tiiwwi  estimator  of  3,  the  vector  target  rotation  rate. 

The  above  diacuuiion  is  appropriate  to  tbs  often-encountered  situation  in  which  the  receiver 
^rations  take  place  in  Ignorance  of  any  prior  knowledge  of  the  relative  probability  that  various 
T.\i.i*e  of  a  will  occur.  Sometimes  one  does  have  such  information  in  the  form  af  an  a  priori 
^.v'tMkbUity  distribution  p(a).  Than,  instead  of  having  the  receiver  produce  as  a  *  that  value  of 
,>  lor  which  p(w/a  and  x)  is  maximum,  one  uses  instead  p(a)  p(w/a  and  x).  Also,  ono  may 
vrtnt  to  assign  "costa"  to  estimation  errors;  for  example,  a  10 -mile  error  in  target  range  might 
U  considered  twice  as  undesirable  in  the  vicinity  of  2000  mile*  range  as  at  2500  miles  range.  A 
second  modification  car  be  made  to  take  care  of  this.  The  form  of  estimation  procedure  that 
t£kes  into  account  cost*  and/or  a  priori  probabilities  is  called  Bayes  estimation.  If  the  cost 
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assignment  and  the  a  priori  probability  p(a)  are  both  eprecd  broadly  over  valuta  of  a,  whereat 
the  function  p(w/a  and  x)  it  relatively  aharply  peaked,  then  the  Bayee  and  maximum-likelihood 
procedures  ere  equivalent.  The  reader  la  referred  to  Ref.  60  for  a  discussion  of  Bayes  estima¬ 
tion;  we  shall  assume  here  that  the  a  priori  probability  and  coat  functions  are  too  broad  to  be  a 
fkctor  in  our  measurement  procedures.  For  example,  if  the  reader  will  refer  to  Fig.  13,  he  will 
probably  agree  that  previous  measurements  of  the  astronomical  unit  have  been  so  widely  spread 
as  to  discourage  assigning  a  peak  of  a  priori  probability  in  a  region  of  a.  u.  value#  any  narrower 
than  several  times  the  probable  error  of  a  new  measurement. 


B.  Point  Nonfluctuating  Targets 

When  the  target  has  a  scattering  function  that  is  much  narrower  along  both  the  r-  and  f-axes 
than  is  the  ambiguity  function  of  any  signal  we  might  be  interested  in  using  to  observe  it,  then 
we  can  say  that  it  is  a  point  nonfluctuating  target. 

The  effect  of  a  point  nonfluctuating  target  on  the  transmission  x(t)  is  simply  to  change  its 
amplitude  in  proportion  to  '/tr  where  <r  is  the  target's  actual  radar  cross  section,  its  delay  by 
the  actual  round-trip  delay  r,  and  its  frequency  by  f,  the  actual  doppler  offset.  Maximum- 
likelihood  estimators  for  all  these  quantities  have  been  found  and  described  in  the  literature. 
Helstram,*0  among  others,  gives  a  reasonably  complete  treatment  of  this  problem,  and  we  shall 
reproduce  his  results  here,  stating  first  the  signal  processing  operations  implied  by  the  expres¬ 
sions  for  the  maximum  -likelihood  estimate,  and  then  giving  the  variance  of  the  estimate.  White 
gaussian  noise  will  be  assumed. 

When  t  and  f  are  known,  and  one  wants  to  obtain  an  estimate  of  received  signal  amplitude 
(which  will  be  proportional  to  Vo),  one  forms  the  quantity 


(V?)*  =  J"  w(t)  x(t)  dt/J  " 


x2(t)  dt 


(47) 


which,  according  to  Sec.  IV- B,  can  be  interpreted  as  the  suitably  normalized  output  at  time  tQ 
at  a  filter  matched  to  the  transmission  [l.e.,  m(t)  -  x(tQ  -  t)  as  in  Eq.  (26)] .  As  one  might  expect, 
the  variance  is  the  noise-to-signal  power  ratio  times  the  amplitude  squared: 


Var  [(«/?)*  ]  s  U-JZ)*  -  <Vw>0|2  * 


(48) 


where  the  subscript  zero  indicates  the  true  value. 

It  can  be  shown  that  Eq.  (48)  is  the  minimum  achievable  variance  of  estimation  of  amplitude 
(by  the  Cramer-Rao  inequality),  so  no  other  estimator  of  V<r  can  do  any  better  than  the  matched 
filter. 

Assuming  that  the  echo  is  of  unknown  amplitude  and  buried  in  white  gaussian  noise,  the 
maximum-likelihood  estimate  of  the  parameter  pair  r  and  f  is  the  value  at  r  and  f  that 
maximizes 


Q(r.f) 


1C 


X*(t)V(t  +  r) 


2 

exp[2*jft]  dt 


(49) 


which  is  the  square  of  the  envelope  at  time  t  +  tQ  of  a  filter  matched  to  x(t),  except  that  the  filter 
is  offset  by  f  cycles  per  second.  As  in  Eq.  (27),  t0  is  the  time  at  which  the  filter  would  have 
maximum  response  if  target  delay  r  were  zero.  In  practice,  one  implements  Eq.  (49)  by  building 
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a  aet  of  matched  filter*,  each  followed  by  an  envelope  squarer,  and  all  having  different  frequency 
offsets  with  the  total  frequency  span  encompassing  all  reasonable  values  of  f.  One  then  examines 
all  the  outputs  over  an  interval  of  time  large  enough  to  span  all  likely  values  of  r,  and  selects 
the  delay  and  frequency  of  the  target  peak  as  representing  the  estimated  values  r  *  and  f* .  Since 
r  and  f  are  continuous  variables,  strictly  speaking  there  should  be  a  continuum  of  stagger-tuned 
matched-filter-envelope-squarer  combinations,  and  also  the  continuum  of  r-values  should  be  ex¬ 
amined.  In  practice,  of  course,  the  f-spacing  and  r-spacing  would  be  set  for  some  fraction  of 
the  expected  rms  error  (square  root  of  the  variance)  in  f  and  t,  respectively. 

Formulas  for  the  range  and  doppler  error  variances  as  a  function  of  duration  TQ  and  band¬ 
width  W_  have  been  found  to  be 
o 
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and 

Var  (f*)  =  . 
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where 
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That  is,  Tq  and  WQ  are  the  "rms"  duration  and  bandwidth,  respectively.  The  variance  expres¬ 
sions  (50)  and  (51)  are  valid  if  three  conditions  are  satisfied:  (1)  the  ambiguity  function  of  x(t), 
the  transmission,  can  be  approximated  near  the  origin  by  parabolas  along  the  r-  and  f-asae; 

(2)  the  axes  of  symmetry  of  the  ambiguity  function  lie  along  the  r-  and  f-axis;  (3)  die  signal-to- 
noise  ratio  is  sufficiently  large  that  the  error  arises  from  random  shifts  of  the  peak  of  the  Amo¬ 
tion  Q  near  the  true  value  of  the  parameter,  and  not  from  mere  remote  peaks  becoming  suffi¬ 
ciently  large.  The  first  condition  is  usually  met,  except  for  such  things  as  ideally  square  pulses - 
The  second  condition  is  met  fay  all  the  function  of  the  sort  typified  in  Fig.  10  except  for  the 
•chirp"  (swept  frequency)  waveform  of  Fig.  10(c).  For  waveforms  that  obey  condition  (2),  the 
range  error  is  not  a  function  of  signal  duration  T  and  the  doppler  error  ie  not  a  function  of  band¬ 
width  WQ.  Also,  the  range  error  la  not  increased  by  Ignorance  of  the  velocity,  and  vice  versa. 

Expressions  (50)  and  (Si)  are  seen  to  involve  only  the  signal-to- noise  ratio  2E#/No  and  the 
widths  of  the  ambiguity  surface.  They  express  the  way  in  which  both  the  notes  end  the  widths  of 
the  ambiguity  function  conspire  to  render  the  estimate  of  the  true  values  of  f  and  r  imprecise 
by  causing  the  maximum  of  Q  in  Eq.  (49)  to  shift  to  the  wrong  piece.  Were  it  not  for  the  noise, 
the  smearing  effect  of  finite  signal  bandwidth  or  duration  would  be  reversible  and  the  position  of 
the  peak  could  be  located  to  arbitrarily  high  precision. 

One  may  wonder  whether  the  processor  of  Eq.  (49)  gives  error  variances  (Eqs.  (SO)  end  (51)] 
that  are  the  minimum  achievable  variances.  Manases6*  has  shown  that  in  a  point  target-gaussian 
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nois*  case  like  the  present  one,  at  large  signal -to-noise  ratios  the  estimates  are  jointly  minimum- 
variance  estimates  (l.e.,  "jointly  efficient")  so  long  as  the  transmission  x(t)  is  a  well-behaved 
function  of  each  parameter  a.  An  interesting  titration  in  which  the  Cramer-Rao  lower  hound  on 
error  variance  is  Qg  physically  achievable  is  the  ease  where  the  transmission  is  a  single  recta  n¬ 
gular  pulse  of  sinusoid  having  a  duration  T.  Again,  one  would  expect  the  range  and  dappler  er¬ 
rors  to  depend  on  pulse  width  and  signal -to-noise  ratios  and,  if  the  noise  is  small  enough,  for 
range  error  to  be  considerably  smaller  than  pulse  width,  as  before.  But  when  one  computes  the 
theoretical  lower  bound  for  range  error  it  turns  out  to  be  aero  for  any  finite  noise  level,  and  this 
does  not  make  physical  sense.  The  derivation  leading  to  Eqa.  (50)  and  (51)  is  not  usable  because 
Condition  (1)  is  violated.  However,  Manases62  has  developed  a  different,  more  heuristic  approach 
to  find  the  error  variances  from  a  matched  filter  processor  end  has  found  that  the  error  variances 
in  this  square-pulse  esse  are 

Var  (r*)  «  (Na/2Eg)2«/2Ty  .  (SOit) 

Vsr  <T* )  =  N{J/2E>  3/»ZT2  .  (51a) 

When  a  number  n  of  similar  but  statistically  independent  measurements  are  edded  together, 
the  variance  is  decreased  by  a  factor  n.  An  example  of  this  would  be  n  pulses,  each  measuring 
range  or  doppler  independently.  However,  In  many  practical  circumstances  the  n  measurements 
may  not  be  independent.  For  example,  thejnilae  train  illustrated  in  Fig.  10(d)  has  the  phase  pre¬ 
served  from  pulse  to  pulse,  so  the  floppier  resolution  available  is  that  appropriate  to  IhedbFfUon 
at  the  train  (as  expressed  by  the  frequency  width  of  the  ambiguity  function),  rather  than  only  l/n 
times  the  variance  to  be  obtained  with  one  pulse.  The  computation  of  error  variances  of  range 
and  dappler  becomes  somewhat  complicated  when  the  target  has  an  appreciable  acceleration  com¬ 
ponent,  when  we  desire  to  measure  this  acceleration,  qr.urhm  skewed  wrfturved  ambiguity  uur- 
faces  S£f  sjlsMNd  (•  teewed  ease  &  iSoemln  Fig.  10(c)).  These  bases  are  particularly  hanl  to 
handle  when  large-TW  signals  such  as  pulse  trains  era  considered  rather  then  single  pulses. 

The  best  available  treatment  of  these  more  eompltoated  esses  has  been  given  by  Kelly.63 

C.  Targets  Spread  only  in  Daisy  or  Doppler 

After  the  point  nonfluctuating  target,  the  next  class  of  radar  targets  in  order  of  complexity 
is  that  for  which  there  is  a  spread  along  either  the  r-axls  or  the  f-axls  (relative  to  the  corre¬ 
sponding  spread  of  tits  signal  ambiguity  function),  but  hot  both. 

For  example,  with  a  point  fluctuating  target,  there  is  a  spread  only  along  the  f-axia,  nnd  the 
scattering  function  is  o(r,  f)  *  <r(f)  g(r  -  r')  where  t'  Is  the  range.  With  such  a  target  then  are 
at  least  four  quantities  one  might  be  Iterated  in  estimating:  (i)  tbs  range  r',  (2)  the  curve  <r(f) 
to  whatever  frequency  resolution  is  interesting,  (5)  the  balk  dappler  shift  of  the  target  as  a  whole, 
given  Hie  shape  aif)  apart  from  the  unknown  translation,  and  (4)  the  width  of  the  dappler  spectrum, 
given  both  the  center  frequency  and  tbs  spectral  shape  o(f)  apart  from  tbs  unknown  width  factor. 

With  the  classical  nonfluctuating  point  target,  the  real  source  of  error  (other  than  such  things 
as  ealttratlon  errors)  was  the  receiver  notes,  represented  by  the  density  NQ.  For  the  fluctuating 
target,  the  echo  fluctuations  art  an  additional  causa  of  error.  As  we  shall  sea  shortly,  those 
fluctuations  place  a  limit  on  the  measurement  accuracy  even  when  additive  received  note)  is  aero. 
Such  contributions  to  the  random  variation  of  the  estimates  will  be  designated  as  "self -noise"  aa 
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before,  according  to  the  convention  employed  here  that  any  random  variations  in  the  output  that 
remain  after  Nq  goes  to  zero  are  to  be  so  designated. 

To  estimate  tie  range  t\  it  appears  plausible  to  proceed  as  with  the  point  nonfluctuating 
target,  namely,  to  build  a  parallel  set  of  signal-to-noise  ratio-maximizing  processors  that  are 
.spaced  in  delay  and  to  pick  as  the  estimate  r'*  the  delay  of  the  processor  exhibiting  the  largest 
output.  Sines  the  target  is  fluctuating,  the  processors  will  now  be  weighted  radiemeters  rather 
that;  the  matched-f  lter-envelope-detector  devices  implied  by  Eq.  (49).  The  r-spscing  between 
adjacent  processor  3  will  be  a  convenient  fraction  of  the  rms  range  error.  For  small  input  S/N, 
and  large  output  S/.'i,  R.  Price*’-*  has  found  the  error  variance  to  be 
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“mere  ip  is:  the  ambiguity  function  of  the  transmission,  as  before,  and  the  output  signal-to-noise 
.  “tio  (S/N)0  is  given  by  Eq.  (34) 

The  e<  tlmation  of  the  function  cr({)  is,  of  course,  just  the  problem  of  spectral  analysis;  i.e., 
G’ven  a  sample  segment  of  a  random  process  of  unknown  spectral  density,  one  is  to  estimate  the 
sp  ,?ct:  »il  d«  osity,  which  in  our  context  is  S(f)  =  c(f)  +  Nq,  where  Nq  is  the  noise  density.  Although 
there  h»  .a  sizable  literature-  on  '.he  subject  of  spectral  analysis  (e.g.,  Refs.  65  and  66),  as  matters 
new  stand  there  is  no  veil -formulated  procedure  leading  to  minimum -variance  estimates  of  the 
spectrum  3(f).  However,  various  recipes  for  estimation  have  been  analyzed  in  detail  and  it  is 
f.’ftrn  possible  to  give  ipprcucimate  expressions  for  the  error  variance. 

Suppose  we  have  a  sample  of  received  signal  T  seconds  long  on  which  we  desire  to  perform 
0  spectral  analysis.  Such  a  sample  might  be  obtained  from  a  radar  experiment  in  wliich  the  trans¬ 
mission  ii!  a  r  inusoid  «l  duration  T.  If  one  either  performs  a  Fourier  series  analysis  of  the  sam¬ 
ple  or  Fouriei'-transfoims  the  autocorrelation  of  the  sample  (these  are  the  two  obvious  ways  of 
proceeding),  the  resulting  estimated  spectrum  will  have  violent  variations  every  1/Tcps,  no  mat¬ 
ter  svhut  ’he  true  spectrum  of  the  random  process  from  which  the  sample  came.  Clearly,  if  one 
insists  on  a  frequency  resolution  ff  tiat  is  as  small  as  l/T,  then  this  jagged  picture  is  the  only 
estimate  possible,  but  is  rneanit  css.  In  practice,  the  desired  frequency  resolution  df  must  be 
made  larger  than  l/T  so  that  a  number  of  different  samples  (or  "degrees  of  freedom"  as  we  shall 

call  them)  will  contribute  to  each  po  nt  on  the  estimated  spectrum. 
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olaokman  and  Tukey  suggest  he  following  procedure  which  tends  to  give  good  estimates 
of  smoothed  values  of  the  true  spectral  densities:  (1)  the  autocorrelation  of  the  sample  is  com¬ 
puted,  ( ;)  this  even  function  <p(r)  is  then  multiplied  by  some  suitably  chosen  even  weighting  func¬ 
tion  or  window  function  D(t)  which  has  a  maximum  at  zero  and  is  zero  outside  the  interval 
(—  t/’flf,  i/fif)  where  fif  represents  tne  desired  frequency  resolution  of  the  estimate,  and  finally 
(5)  the  product  <p{r)  D(t)  in  Fourier  -transformed  to  get  the  spectral  estimate  S*(f).  This  is 
equivalent  to  convolving  the  transform  of  D(t)  with  the  transform  of  the  specimen  correlation 


l  Thta  expression  is  valid  only  for  targets  having  a  nonzero  doppler  spread;  It  does  not  reduce  to  Eq.  (50)  as  doppler 
width  900s  to  zero. 
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function  <?(t).  The  proper  choice  of  D(t)  is  a  compromise  between  effectively  throwing  away 
some  of  the  data  by  suppressing  the  extremities  of  <p[ r)  and  retaining  the  violent  variations  of 
estimated  spectrum  with  frequency.  For  gaussian  processes,  Blackman  and  Tukey  give  the  fol¬ 
lowing  expression  for  error  variance  of  the  spectrum  at  a  given  frequency  when  the  true  S(f) 
varies  slowly  with  frequency  {relative  to  i/T): 

Var  [S*{f)J  -  —■ jff*  (f)-l  ■  (53) 


The  average  of  the  estimate  S*(f)  is  essentially  the  true  SQ(f).  Much  the  same  will  be  ob¬ 

tained  if  the  spectral  estimate  is  performed  by  making  a  Fourier  analysis  of  the  sample  and  then 
adding  together  the  powers  in  the  T  flf  frequency  increments  contained  in  Of  cps. 

Equation  (53)  places  in  evidence  the  well-known  rule  of  thumb  that  the  ratio  of  variance  to 
squared  mean  is  inversely  proportional  to  the  number  of  independent  samples  or  degrees  of  free¬ 
dom  (roughly  2T ..Of  for  a  spectral  interval  Of  wide  observed  for  T  seconds).^  Thus,  if  one  is 
measuring  S(f)  of  a  radar  astronomy  target  he  must  try  to  make  the  total  observation  time  T 
several  times  longer  than  the  reciprocal  of  Of  the  desired  frequency  resolution,  as  as  to  increase 
the  number  of  degrees  of  freedom  contributing. 

For  our  radar  astronomy  purposes  the  usur!  recipe  for  spectral  analysis  (of  which  we  have 
just  presented  one  that  has  been  elaborately  studied  for  practical  applicability)  suffers  from  the 
defect  that  it  provides  an  estimate  of  ST)  =  <r(t)  +  Nq  and  not  a(f)  Itself.  Without  careful  calibra¬ 
tion  it  may  be  difficult  to  know  just  what  baseline  level  to  subtract  out  of  the  resulting  plot  S*(f) 
in  order  to  isolate  the  echo  component  <r(f). 

The  problem  of  estimating  the  frequency  shift  f^  of  a  spectrum  of  known  shape  S(f)  from 
sample  of  durr.tion  T  of  a  random  gaussian  function  of  true  spectral  density  S(f  -  frf)  has  been 
treated  by  Swerling^7  and  Levin.^®  Swerling  has  suggested  a  procedure  in  which  r  <-  .  es¬ 

timate  S*(f)  obtained  by  the  Blackman-Tukey  method  just  described  is  operated  on  by  a  frequency 
weighting  factor  or(f)  so  that  the  estimate  f^  is  given  by  f  S*  (f)  a(f)  df.  The  functior  ■[!)  is  chosen 
from  prior  knowledge  of  the  S(f)  shape  to  minimize  the  error  variance  of  f^  in  terms  of  error 
variances  of  S*(f).  Since  the  latter  quantity  S*  (f)  is  obtained  to  begin  with  by  a  more  or  less 
ad  hoc  procedure,  as  we  have  seen,  the  use  of  S*  to  get  f^  is  not  guaranteed  to  produce  minimum- 

variance  estimates  of  f .. 
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Levin  sought  to  find  a  maximum -likelihood  estimator  of  f^  directly  and  has  examined  the 
Cramer-Rao  lower  bound  on  error  variance.  He  finds  that 
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where  <r'(f)  is  the  derivative  of  o  with  resp'’0*  *c  frequency.  Note  that  rapid  variations  of  o(f) 
with  f  tend  to  make  estimates  of  frequency  shift  more  reliable,  as  one  would  expect.  The 
maximum-likelihood  estimator  of  is  that  which  picks  that  value  of  i'd  for  which 
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t  The  factor  2  occurs  because  at  each  frequency  the  signal  is  characterized  by  two  quantities,  e  g.,  an  amplitude 
and  a  phase. 


is  a  minimum,  where  a_  and  b_  are  the  Fourier  line  and  cosine  coefficients  of  a  Fourier  series 
expansion  at  the  sample  function.  In  other  words,  one  processes  the  received  sample  to  gen¬ 
erate  a  Jagged  function  representing  the  power  in  the  successive  Fourier  components  and  than 
convolves  this  with  the  reciprocal  at  the  known  shape  S(f)  *  NQ  +  at t)  to  find  the  shift  fd  for  which 
a  minimum  output  is  observed.  As  the  product  at  observation  time  T  and  B  at  S(f)  becomes 
large,  this  maximum-likelihood  estimate  is  a  minimum-variance  estimate  and  Eq.  (54)  converges 
to  tile  equality  condition. 

The  frequency  width  at  a  spectrum  can  be  estimated  by  a  related  procedure  also  worked  out 
fay  Levin.**  (Recall  that  one  method  at  determining  the  rotation  vector  at  a  planet  is  to  measure 
tile  width  at  the  echo  power  spectrum  at  various  times  during  the  synodic  period.)  Let  the  true 
spectral  shape  be  Sfh(f  -  fj)  J  -  efh(f  -  fd)  J  +  where  fd  is  the  center  frequency  and  h  is  a  acale 
factor  on  width.  If  only  h  is  unknown  and  is  to  be  estimated.  Levin  finds  that  the  Cramer-Rao 
lower  bound  is 
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is  very  similar  to  Eq.  (54).  Note  again  that  rapid  variations  at  target  behavior  a  wtth  fre- 
fhetHtate  measurement  of  the  width  factor  h;  also  that  the  quadratic  factor  means  that 
at  frequencies  distant  from  ths  band  esnter  are  mere  helpful  than  i 
The  maximum  -likelihood  esilsisliw  of  h  picks  thet  value  at  h  for  which 
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is  s  K  being  a  certain  constant,  hi  other  wards,  cm  processes  ths  received  sample 

to  generate  the  power  in  the  eucceeaive  Fourier  components  (ter  example  fay  using  s  narrow-band 
filter),  and  then  computes  the  convolution  for  fixed  shift  at  this  Jagged  function  with  the  reciprocal 
at  the  profile  sfh(f  -  fd)  J.  choosing  various  values  at  h  until  the  consolation  is  minimised.  Again, 
the  bound  [Eq.  (54)]  applies  with  the  equality  sign  when  MT  id  large. 

The  other  form  of  target  spread  in  one  domain  only  is  the  deep  aeaftaetmting  target  whose 
scattering  function  is  <r(r,  f)  «  s(r)  ttt  —  f).  All  ths  notions  Just  presented  abost  errors  at  esti¬ 
mating  t',  a(f).  fd  and  h  apply  through  a  certain  duality  relationship2*  to  the  estimation  of  A, 
a(r),  rd,  etc.  for  the  deep  nonfluctuating  target.  Instead  of  dealing  with  s  sample  function  T 
seconds  long,  obtained  perhaps  by  illuminating  the  target  with  a  sins  wavs,  the  target  is  affec¬ 
tively  Illuminated  with  a  single  impulse  end  phase  end  amplitude  at  the  echo  components  at  vari¬ 
ous  frequencies  are  measured  over  s  bandwidth  W.  Instead  at  having  a  narrow-band  received 
signal  exhibiting  fluctuations  of  amplitude  and  phase  every  1/B  seconds  in  time  throughout  T,  we 
have  fluctuations  of  amplitude  and  phase  every  l/Lcps.  throughout  W..  In  situations  where 
previously  we  dealt  with  a  Fourier  series  analysis  of  a  time  function  T  seconds  long  we  deal  with 
time  samples  of  the  echo  waveform  every  1/2  W  seconds.  The  detailed  application  of  this  duality 
to  gat  answers  for  the  deep  nonfluctuating  target  will  not  be  pursued  any  further  here.  This  duality 
has  been  developed  in  detail  by  Bello.** 
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D.  Dmp  TlMkating  Targets 

As  the  remder  will  eee  from  c'enparing  the  materiel  in  the  previous  two  sections  the  appli¬ 
cation  of  statistical  estimation  theory  to  the  classical  nonfluctuating  point  target  is  rather  com¬ 
plete,  but  our  present  understanding  of  the  mere  complicated  singly  spread  target,  although 
filling  out  rapidly,  is  still  somewhat  fragmentary.  For  the  doubly  spread  or  deep  fluctuating 
target  the  situation  is  even  worse.  Only  a  few  miscellaneous  results  are  available. 

In  practice  one  might  want  to  estimate  the  following  attributes  at  the  true  scattering  func¬ 
tion  o(r,  f):  (1)  the  center  frequency  and/or  center  delay,  given  the  exact  shape  of  e(r,  f);  (2)  the 
frequency  width  and/or  delay  width,  given  all  other  imormatlon  on  e(r,  f);  (3)  the  function  e(r,  f) 
itself,  given  no  a  priori  data  about  it.  Of  these,  only  the  first  and  last  have  received  much 
attention. 

The  simplest  method  of  pleasuring  e(r,  f)  is  the  rather  direct  one  that  was  discussed  in 
Table  I,  the  use  of  a  simple  matched  filter  or  correlation  detector  with  preset  r  and  f  offsets 
to  probe  the  scattering  function  at  the  point  r,  f .  This  is  the  method  employed  in  Fotteagttl’e 
moon-mapping  work.*'*  Suppose  we  dee  ire  to  do  such  an  experiment  and  achieve  s  frequency  reso¬ 
lution  at  St  and  s  delay  resolution  Sr;  that  is,  we  with  to  measure  the  average  eche  power  In  a 
cell  of  the  r  -  f  plane  having  the  dimensions  gf  and  Sr.  If  we  make  the  signal  bandwidth  t/Sr  sad 
the  duration  t/st.  than  ww  jpt  in  tile  output  squared  envelope  just  such  an  ohssrvutlen  of  (he  de¬ 
sired  cell  [Eqs .(31)  end  (32)],  but  the  number  of  degrees  of  freedom  of  the  observation  is  only  two 
end  the  variance  equals  the  mean,  ea  with  the  stapler  spectral  analysis  ones  described  If 
Eq.  (53).  H,  hesrever,  a  bandwidth  m/Sr  and/or  a  duration  a/St  la  used  (m  and  n  >  i^aai  the  am 
separate  received  squared  envelopes  ere  added  together  and  the  square  root  tabes,  then  the  num¬ 
ber  of  degrees  at  freedom  is  2mn  and  the  variance  at  estimation  of  e(r,  f)  is 

Ver  [«t*(t, f)J  «  (3«) 

as  with  Eq.  (53). 

One  defect  of  this  simple  procedure  is  its  behavior  when  the  target  approaches  the  over¬ 
spread  condition,  that  is,  when  the  product  of  total  doppler  spread  B  and  total  multipath  spread 
L  is  near  unity.  As  we  have  seen  earlier,  the  operation  of  such  a  matched  filter  or  correlation 
detection  method  can  be  visualized  as  overlaying  onto  the  surface  <r(r,  f)  the  surface  $  (r,  f)  (with 
appropriate  delay  and  frequency  offsets)  and  then  adding  up  all  the  volume  under  the  produet 
function.  As  long  ee  the  ambiguity  function  g2  has  most  of  ita  volume  in  the  central  peak,  very 
little  output  power  comes  from  any  piece  but  in  or  very  near  the  desired  cell  in  the  r,  f  plane, 
no  matter  how  large  the  area  BL  covered  by  the  function  a.  However,  as  the  invariance  relation 
Eq.  (30)  assures  us,  such  a  concentration  of  the  volume  in  g2  will  occur  only  when  the  base  of  the 
central  peak  has  an  area  of  approximately  unity,  as  with  the  waveforms  of  Figs.  10(a),  (b)  and  (c). 
But  suppose  it  la  necessary  to  us*  one  of  the  large-TW  waveforms  of  Figs.  10(d)  and  (e).  This 
would  be  the  procedure  if  it  is  desired  to  resolve  s  cell  in  the  r,  f  plane  of  size  Sr  Sf  ■  i/TW  «  1 . 
Because  at  the  secondary  peaks  at  Fig.  10(d)  or  the  low-lying  skirts  of  Fig.  10(e),  it  the  area  BL 
covered  by  the  scattering  function  is  great  enough,  then  the  processor  output  power  due  to  the 
desired  region  of  the  r.  f  plane  (the  region  at  the  central  peak  of  g2)  will  have  added  to  it  gglf- 
nolee  power  from  undesired  regions.  (Once  again  we  use  the  term  self -noise  to  designate  output 
components  that  do  not  go  to  sero  even  if  the  additive  input  noise  does.)  Clearly,  whether  such 
components  will  appear  in  the  output  depends  on  the  scattering  function  area  BL  and  on  the  type 
of  ambiguity  function  chosen  to  probe  it. 
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Let  ua  illustrate  these  ideas  by  picking  some  numbers  from  a  specific  radar  astronomy  ex¬ 
periment,  Suppose  we  wish  to  make  a  radar  map  of  Mars  using  a  pulse  train  [Fig.  10(d)],  The 
multipath  depth  L  from  nose  to  limb  is  24  msec,  which  means  that  observations  of  the  received 
signal  made  more  frequently  than  once  every  24  msec  will  be  noisy  due  to  time  overlap  of  the 
echoes.  The  echo  power  near  the  limb  is  so  small  relative  to  that  from  the  nose,  that  a  periodic 
sampling  of  the  echo  from  the  nose  would  not  be  significantly  corrupted  with  echoes  from  the  limb 
until  the  sampling  interval  became  much  leas  than  24  msec,  but  the  converse  would  not  be  true; 
samples  of  the  echo  from  near  the  limb  would  be  badly  corrupted  by  strong  self -noise  from  the 
nose  as  soon  as  the  sampling  period  was  reduced  below  24  msec.  Now  the  carrier  frequency 
above  which  the  total  doppler  spread  of  Mars  is  1/24  kcps  turns  out  to  be  26.5  Mcpe.  In  order  to 
compute  a  clean  echo  spectrum  l/24kcps  wide,  the  ability  to  get  a  clean  sample  ones  every 
24  msec  must  be  preserved,  and  we  see  that  this  condition  begins  to  be  violated  at  a  carrier  fre¬ 
quency  of  25.5  Mcps;  i.e.,  this  is  the  lowest  frequency  for  which  BL  «  1,  if  fay  B  we  mean  the 
total  limb-to-llmb  doppler  spread  and  by  L  we  mean  the  total  multipath  duration.  (The  fre¬ 
quencies  for  which  BL  (so  defined)  equals  unity  are  given  in  Table  II  for  other  radar  astronomy 
targets.) 

Now  suppose  we  are  operating  at  440  Mcps  so  that  BL  is  very  high,  actually  18,  B  being 
750  ops.  If  we  do  net  need  more  than  about  16  cells  in  our  radar  map,  we  can  send  a  simple 
TW  =  i  pulse  (Figs.  10(a)  and  (b)]  of  duration,  say,  6  msec.  and  by  averaging  the  received  power 
at  the  desired  receiver  r-  and  f-offseta  get  an  obeervation  of  echo  power  hi  a  call  of  the  r-f 
plane  that  if  dr  *  6msec  by  M  ■  l/6kcps  in  area. 

We  conclude  that  when  BL  is  much  smaller  Bam  unity,  a  large-TW  signal  oaa  ha  nhnseo  such 
that  the  self-noise  lc  negligible.  For  example,  in  Fig.  10(d),  so  long  aa  the  scattering  function 
does  not  spread  outside  the  area  enclosed  by  the  dotted  lines,  self -noise  do*  to  nonoentral  peaks 


TASLE  II 

CAMtMK  FMQUENCIES  AiOVE  WWCH  VARIOUS  TARGETS  AXE  OVEXSWAO 
(XL  Product  far  an  Arbitrary  Frequency  It  1^) 

Target 

■SSI: 

Mean  (at  maximum  librettos  tats) 

f§  ■  3130 

Mercury  (amuming  88-day  rotation  period) 

4280 

Venus  (amuming  225-day  rotation  period) 

1180 

Earth 

39.7 

Mart 

26.5 

Jupiter 

1.58 

Saturn 

0.280 

Uranut 

2.08 

Neptune 

3.77 

Pluto  (amuming  16  hour  rotation  period) 

2. 25 
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will  be  negligible.  We  further  conclude  that  when  BL  is  much  greater  than  unity,  if  one  can  be 
satisfied  with  a  resolution  cell  at  size  dr  df  »  1,  then  a  TW  =  1  pulse  [Fig.  10(a),  (b)  or  (c»  can 
be  used,  and  because  of  the  central  concentration  of  the  volume  in  the  ambiguity  function,  once 
again  the  self -noise  will  be  negligible .  It  is  the  in-between  values  of  BL  of  the  order  at  unity 
that  produce  self-noise  difficulties.^  For  illustration,  we  have  plotted  in  Fig.  15  several  curves 


li-w-mti 


Fig.  15.  Hie  manner  in  which  ielf-noiie  corrupt!  the  estimate  of  the  Mattering 
function  vdien  the  simple  matched  filter  procedure  Is  utad  and  the  target  is  moder¬ 
ately  overspread. 


of  setf-noiae-to-algnel  ratio  as  a  function  of  BL,  assuming  that  tha  scattering  function  being  ob¬ 
served  haa  a  gauasian  shape  in  both  f  and  r  directions  {i.s.,  cr(r.t)  •  exp  -  [t2/L2  +  f2/B2)}, 
and  that  the  siae  of  the  desired  resolution  cell  dr  St  is  0.1  BL  to  10~*  BL.  It  is  assumed  that 
the  waveform  in  use  is  the  pseudonolae  whose  ambiguity  function  ia  idealised  at  the  right  of 
Fig.  10(e). 

Use  two-frequency  method  of  measuring  the  scattering  function  also  encounters  difficulties 
as  the  target  becomes  overspread,  as  is  quite  easily  visualized.  By  referring  to  the  fourth 
entry  of  Table  I  we  are  reminded  that  one  first  obtains  41  (Af.  At)  by  cross-correlating  the  outputs 


t  We  may  conjecture  that  It  is  Impossible  to  realize  a  physical  waveform  whoso  ambiguity  function  has  loss  than 
a  certain  asount  of  volume  in  an  area  muah  larger  than  tha  order  af  unity  surrounding  the  central  peak.  Tha 
volume  may  occur  in  the  central  peak,  as  in  Fig.  10(a),  (b)  or  (c);  in  separate  peaks,  as  In  Fig.  10(d),  or  In  wide 
foothills  as  In  10(a).  (Such  an  area  with  a  low  none  antral  volume  is  tha  dotted  region  of  Fig.  10(d).)  This  con- 
lecture  remains  unproved  to  dote. 
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of  two  filters  tuned  to  soperetely  receive  two  sinusoid*  transmitted  Af  apart,  and  then  Fourier-  ' 
transforms  A  (Af,  At)  to  get  <r(r,  f).  Now  A  varies  with  Af  at  roughly  a  rate  l/L;  that  ie,  we 
will  want  to  measure  A  for  frequency  spacing*  Af  as  low  as  i/Lcps.  But  each  sinusoid  is 
smeared  by  B  ops,  so  a  clean  (self-noiseless)  signal  separation  will  be  possible  only  if  B  >  l/L, 
i.e.,  if  the  target  ie  underspread. 

There  is  a  class  of  methods  of  measuring  the  scattering  function  <t(t,  f)  that  appears  to  be 
devoid  of  the  problem  of  self -noise  due  to  overspreading,  and  has  other  advantages  over  the 
matched  filter  (or  correlator)  method  we  have  described,  even  though  the  latter  is  easier  to  im¬ 
plement.  These  methods  are  the  fourth  moment  method,70,71  the  S/N-maximization  method.72 
and  the  least-squares  method.  Each  requires  elaborate  receiving  equipment  and/or  computer 
processing.  These  three  approaches  are  worth  taking  a  bit  of  space  to  describe  here,  because, 
although  their  advantages  and  usability  are  not  completely  understood  at  present,  they  are  being 
actively  investigated  and  are  not  gj)  hoc  methods  (as  is  the  matched  filter  approach)  but  have  op¬ 
timum  propertied,  as  we  shall  now  explain. 

The  fourth  moment  method  proposed  by  Kallath  involves  a  receiver  which  processes  the  re¬ 
ceived  signal  wft)  so  as  to  compute  *he  quantity 

<w(t)  w(t  ♦  or)  xft -0)  x(t  +  a -d)>  , 

in  which  the  transmission  xft)  is  a  white  noise.  As  before,  <  >  represents  the  infinite  time 
average.  Karting  from  each  a  measurement  of  this  fourth  moment,  we  can  than  the 

estimated  scattering  function  s*(r,  f). 

Levin's  least -squares  method  is  a  more  general  procedure  of  which  the  fourth  moment  method 
is  the  special  case  for  very  long  integration  time  and  for  x(t)  a  white  noise,  hi  die  least  squares 
procedure  see  processes  the  finite  length  sample  wft)  to  get  the  same  fourth  moment  (moept  that 
the  integration  time  is  nqw  finite).  This  fourth  moment  ie  then  openated  on  by  an  explicit  linear 
transformation  Involving  tba  known  transmission  xft)  and  the  result  t  is  a  quantity  which,  when 
Fourier-tranef earned  once,  is  the  desired  <t*(t,  f). 

Price's  B/M  smnhnlaaMoa  method  consists  of  building  a  rate  radiometer  (Fig.  11(d))  with 
adjustable  fitter  characteristic*  |H|(Q|  at  each  tap.  These  ehsracterietioe  ere  then  «f)nstad  by 
some  trial -and-error  pro endure  until  the  trial  output  signal -to-aotee  ratio  is  maximised.  Re¬ 
calling  that  the  [■t(f)|2|s  that  nModmisn  signal-to-ndse  ratio  ere  a  repredaeHau  of  the  aeattwlng 
function  efr.f)  (Eq.  (35)),  we  have  thus  orsatad  an  eethnste  e*(r,f)  la  (he  "test*  settings  af  the 
fitters.  Price  shows  diet  at  lour  input  aljgml  to  notes  ratio*.  this  method  produoee  a  maxhusn- 
likelihood  estimate.  Levin  has  shown  that  the  least  squares  end  S/M-maxhnisatica  methods  are 
equivalent,  end  therefore  that  at  low  Input  si^ml-to- noise  ratios  they  are  both  maxim um- 
litelihood  estimators. 


t  It  turn  out  that  at  ail  dpsNf  nolle  mdas  Ibl*  proeedure  b  equivalent  le  perfermlng  e  Isait-iquorw  fit  be¬ 
tween  die  ceverlenee  setrix  of  the  maple  uwvelens  (sell  it  W)  end  e  trial  oeverlence  netrix  S .  That  h,  the 
eleeuMto  *t|  of  S  ere  adfinied  le  srinlmiss 

evV  • 

S  Involve*  both  die  knew*!  end  fixed  tens*  hd an  xft)  and  the  edjisleble  trial  target  teeter  hq  function. 
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At  present,  not  much  is  known  about  the  variance  of  estimates  obtained  by  these  methods, 
it  seems  plausible  that  some  specification  of  the  intended  degree  of  range  resolution  tr  and  dcp- 
pler  rendition  af  must  be  made,  and  that  then  the  variance  of  the  estimates  will  go  down  as 
fir  Si  'W  the  number  of  degrees  of  freedom  per  resolution  cell,  goes  up  (T  being  the  sample 
duration  and  W  being  its  bandwidth). 

Prici.^4  has  investigated  the  problem  of  estimating  r  ,  the  center  range  of  a  deep  fluctuating 
target  whose  scattering  function  is  known  apart  from  this  range  parameter.  The  result  is  a  more 
general  case  of  that  given  in  Eq.  (52),  and  amounts  to  the  following.  The  estimator  is  an  array 
of  rake  radiometers  set  for  various  delays.  The  estimate  r*  is  that  corresponding  to  the  proc¬ 
essor  exhibiting  the  largest  output.  Under  the  conditions  of  low  input  and  high  output  signal-to- 
ttoise  ratios,  and  assuming  that  the  output  in  the  absence  of  noise  is  parabolic  in  the  neighborhood 

vi  the  trw  r  ,  Price  findr  that 
o 


Var  <r*)  *  <S/N)J4 


vh.  suiotiunr  and  conclusions 

In  this  report  the  author  has  attempted  to  summarize  the  techniques  by  which  various  astro¬ 
nomically  Interesting  properties  at  remote  objects  may  be  studied  by  processing  at  received  ra¬ 
dar  signals,  hi  the  esse  of  "soft*  targets  such  as  ionized  regions,  clouds  of  small  particles, 
sod  so  forth,  little  that  is  unfamiliar  waa  said  here;  conventional  Bounding  technique#  are  avail¬ 
able  lor  estimating  electron  density  and  magnetic  field. 

Hard  targets  are  quits  interesting  from  the  signal  processing  viewpoint,  since  by  the  choice 
of  suitable  transmitted  signals  sad  receiver  operations  a  number  of  target  properties  may  be  de¬ 
termined  .  One  of  the  meet  important  at  these  is  the  statistical  character  of  the  rough  target 
eurface.  In  this  report  a  number  of  procedures  were  discussed  for  determining  the  radar  angu¬ 
lar  power  spectrum  <rQ{f)  at  either  a  portion  of  the  target  or  the  whole  surface  (the  Utter  under 
the  ajeumption  at  statistical  uniformity  at  the  surface).  The  function  rQ  expresses  directly  the 
second-order  statistics  of  the  random  spatial  irregularities  in  the  echo  wavefront,  and  it  ex- 
j,  ,-es  tea  indirectly  the  second-order  statistics  of  the  target  terrain  Itself.  The  reUtion  between 
sUtl.  lies*  properties  of  wavefront  and  terrain  was  not  dealt  with  here. 

Measurement  of  or0(?)  of  uniform  objects  can  be  made  by  short -pulse,  long-puUe,  spaced - 
frequency  or  interferometer  experiments;  also,  the  roughness  can  be  inferred  grossly  from 
poUrimetry  measurements.  This  last  type  of  procedure  requires  at  this  writing  much  further 
theoretical  work  before  one  can  say  Just  how  the  amount  of  signal  depolarization  depends  on  sur¬ 
face  roughness  (suitably  defined). 

?or  nonuniform  objects  one  is  interested  in  determining  oD(p)  for  various  regions.  Isolation 
of  various  target  regions  in  antenna  beamwidth  is  usually  out  of  the  question  because  of  the  greet 
target  ranges;  also,  it  was  found  to  be  impossible  to  perform  the  laolation  at  the  basis  of  differ¬ 
ent  nations  of  the  corresponding  parts  of  the  diffraction  pattern  observed  at  the  earth.  However, 
resolution  of  target  regions  can  be  effected  by  using  large  time-bandwidth  signals  such  as  pulse 


////  p2(r,f)  u(t'  -  t,  f  +  f)  ir(r'.f)  dr,  dfdr’df' 


d‘ 

‘5? 


(59) 


////  »<*'  -  r.  f  *  «  »(r'  -  X,  V)  dTdfdr'df' 


X»0 
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trains  or  pseudonoise  signals  to  isolate  echo  power  returning  at  various  points  in  a  hypothetical 
range-vs-doppler  plane  (that  is,  to  estimate  the  "scattering  function"  of  the  target).  Moreover, 
there  is  a  simple  distortionless  geometrical  mapping  between  points  on  the  target  and  points  in 
this  plane. 

There  are  a  number  of  unanswered  questions  about  processing  techniques  for  isolation  of 
target  regions.  For  example,  the  way  in  which  one  might  apply  in  practice  certain  measurement 
procedures  that  are  known  to  be  statistically  optimum  is  unclear.  It  is  unclear  to  what  extent 
one  can  avoid  the  added  noise  injected  into  the  receiver  output  when  examining  an  overspread 
target.  The  added  complexities  of  processing  and  the  added  possibilities  for  sensitive  mapping 
experiments  using  range -doppler  trajectories  have  not  been  adequately  studied. 

The  measurement  of  target  rotation  speed  and  axial  orientation  was  discussed  and  various 
fairly  straightforward  methods  were  outlined  for  partially  determining  these  quantities  by  proc¬ 
essing  of  single  station  or  interferometer  echo  signals. 

The  radar  measurement  of  several  other  target  attributes  (e.g.,  target  shape)  was  discussed 
briefly. 

Explicit  procedures  for  actually  carrying  out  the  observations  were  dealt  with  in  Sec.IV-F, 
and  particularly  in  Table  I. 
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APPENDIX  A 

A  DISCUSSION  OF  CORRELATION  FUNCTIONS  AND  ANGULAR  SPECTRA 
OF  INTEREST  IN  RADAR  ASTRONOMY 


by 

Martin  Balser 

Consider,  for  simplicity,  the  one-dimensional  problem  of  a  component  of  the  field  £  (x,  0) 
given  on  the  plane  z  =  0  as  a  function  of  the  coordinate  x  along  the  plane.  Here,  x  and  z  are 
distances  expressed  in  multiples  of  a  wavelength.  We  assume  £  and  £ 2  are  integrable;  i.e., 
the  field  is  essentially  limited  to  a  given  region  or  aperture  in  the  plane.  Thus  £  has  a  Fourier 
transform,  given  by 

'W'(l)  =  f  £ (x,  0)  exp [— j27rix]  dx  (A-l) 

and,  inversely, 

<?(x,°)  =  ^j'  lV(t)  exp[j2rfx]  df  (A-2) 

It  is  easy  to  see  that 

£(x,  z)  =  =rj-  J  exp[j2ir(fx  +  nz)]  ds  ,  (A-3) 

where  n  =  —  f2,  is  a  solution  to  Maxwell's  equations,  being  a  superposition  of  plane  waves, 

the  one  of  amplitude  X^f)  df  traveling  in  the  direction  ip,  where  f  =  sin  ip  and  n  =  cos  ip.  The  well- 
known  result  that  structure  in  the  aperture  field  which  is  finer  than  the  wavelength  does  not  prop¬ 
agate  is  seen  formally  in  the  observation  that  such  structure  corresponds  in  Eq.  (A-2)  to  values 
of  |l|>  1,  which  results  in  Eq.(A-3)  in  imaginary  values  for  n  and  hence  in  damped  waves. 

A 

It  can  further  be  shown  that  at  great  distances  (compared  to  the  size  of  the  aperture)  from 
the  source  plane,  the  field  in  Eq.  (A-3)  becomes  £(x,  z)  ~"W\t ),  where  f  =  z/*Jx2  +  z2;  i.e., 
I^(f)|2  is  the  angular  spectrum,  indicating  the  power  radiated  in  the  direction  ip  =  sin"1  i  . 

(The  aperture  distribution  is  assumed  to  be  in  the  vicinity  of  the  origin.)  Also,  the  correlation 
in  the  angular  pattern 

£{x,  0)  £  *  (x1,  0)  exp(— j2xpx  +  j2ir(p  +  f )  x']  dx  dx'  dp 

|  £  (x,  0)  |  2  exp  [j2jrfx]  dx  ,  (A-4) 

i.e.,  the  Fourier  transform  of  the  aperture  power  distribution  or  illumination.  Since  the  values 
of  f  can  be  immediately  related  to  distances  x  along  the  observation  plane  z  =  z  (denoted  zQ  in 
Fig.  2  of  the  text),  this  theorem  relates  a  "distance  correlation"  to  the  aperture  brightness. 

A  somewhat  different  concept  of  distance  correlation  along  the  observation  plane  is,  how¬ 
ever,  usually  considered  in  practical  situations.  First  of  all,  the  measurement  is  taken  quite 
locally,  not,  as  in  Eq.(A-4),  through  all  values  of  f  (or  x  in  the  observation  plane).  Second, 


J  TV(p)W*(p  +  t)  dp  =  JjJ 


c 


2jt 
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the  averaging  (for  most  practical  cases,  in  time)  used  to  determine  the  correlation  implies  a 
statistical  process,  whereas  the  previously  cited  results  are  precisely  correct  for  any  individual 
distribution  <f(x,  0).  Accordingly  it  may  be  asked  whether  the  same,  or  a  similar,  result  exists 
for  a  correlation  function  more  nearly  related  to  this  practical  concept.  To  see  whether  this  is 
so,t  consider  the  aperture  field  distribution  to  be  a  sample  drawn  from  an  ensemble  of  the  form 

£ (x,  0)  =  S(x)  3C(x)  ,  (A-5) 

where  JC(x)  is  a  stationary  random  process  in  x,  and  S(x)  is  a  fixed  envelope  which  determines 
the  power  distribution  in  the  aperture  plane  ( |  S(x)  | 2  is  the  brightness  distribution).  The  random 
function  is  assumed  to  have  the  ensemble  averages  (expected  values)  <3C(x»  =  0,  <3CZ(x)>  =  i, 
<JC(x)  X*  (x  +  f)>  =  p(t).  If  now  we  form  the  product  ' 

OO 

^(p)^*(p  +  i)  =  S(x,  0)  £*{x\  0)  exp [— j2rpx  +  j-*(p  +  t)  x1]  dx  dx' 

-CO 

to 

=  Jj*  S(x)  S*(x')  X(x)  X*  (x  ’)  exp  [— j2irp(x  — x1)]  exp(  j2xfx)  dx  dx'  (A-6) 

—  OO 

and  take  the  ensemble  average,  we  find 

to 

<K/(p) >/*(p  +  1)>  =  JJ  S(x)  S*(x')  p(x-  x')  exp [— j2rp(x  —  x') ]  exp{  j2jrfx]  dx  dx1  .  (A-7) 

•OO 

It  is  reasonable  to  assume  that  the  over-all  aperture  power  distribution  |  S(x)  | 2  changes  very 
little  over  a  correlation  distance  (i.e.,  the  range  of  non-negligible  values  of  p),  so  that  in  Eq.  (A-7) 
we  can  replace  S(x')  S*(x)  =  |  S(x)  |  2.  Changing  the  variable  x>  -  x  =  u,  we  find 

0/(p)X^*(p  +  l)>  =  C  p(u)  exp[j2rpu]  du  f  |  S(x)  | 2  exp  [j2xfx]  dx  (A-8) 

03  v.eo 

An  incidental  result  of  Eq.(A-8)  is  found  by  setting  I  =  0;  then 

<lMp)  1 2>  =  K  J  p(u)  exp[j2rpu)  du  ,  (A-9) 

(where  K  =  f_n  \  S(x)  |  2  dx),  which  states  that  the  angular  spectrum  radiated  from  the  aperture 
is  the  Fourier  transform  of  the  correlation  function  of  the  field  in  the  aperture.  (This  was  also 
shown  by  Bramley7*).  The  principal  result  is,  however,  the  other  integral,  i.e., 

<l/(p)?J/*(p  +  l)>  =  <|H/(p)|Z>  j*  I S(x) |  2  exp [j2xfx]  dx  (A-10) 

This  result  is  very  similar  to  the  one  in  Eq.  (A-4)  in  that  it  relates  the  angular  (or  distance)  cor¬ 
relation  to  the  Fourier  transform  of  the  aperture  illumination,  but  in  Eq.  (A-10)  the  statement  is 
made  for  each  value  of  direction  cosine  p  separately,  rather  than  integrated  over  all  values  of  p. 


t  ThU  it  etMntiolly  the  approach  uted  by  Ratcliffs  (Ref.  2,  p.  223). 
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For  different  value*  of  p,  the  angular  correlation  function  has  the  same  form  and  differs  only  by 
the  constant  '■■f  proportionality,  the  angular  spectrum. 

At  this  point  the  effect  of  the  assumption  that  S(x)  S*  (x1)  ~  |  S(x)  |  2  can  be  seen.  [It  was  this 
assumption  that  allowed  the  separation  into  a  product  of  the  double  integral  in  Eq.  (A-7).]  To  ne¬ 
gate  this  assumption,  i.e.,  to  state  that  the  correlation  distance  of  the  aperture  field  is  compa¬ 
rable  ;.o  the  total  size  of  the  aperture,  would  mean  that  the  angular  correlet’on  distance  is  com¬ 
parable  to  the  width  of  the  angular  spectrum.  Thus  [from  Eq.  (A-10)],  the  received  power  level 
in  the  observation  plane  would  change  over  distances  comparable  to  the  space  correlation  distance, 
destroying  the  normal  assumption  of  "local  stationarity "  used  in  measuring  correlation  functions. 

Note  that  all  the  previous  analysis  has  based  the  angular  coordinate  system  at  the  origin,  lo¬ 
cated  at  the  aperture.  Actually,  by  defining  a  new  set  of  angles  from  the  point  of  observation  as 
origin,  it  can  be  seen  that  Eq.  (A-10)  is  the  same  as  Eq.  (A-9),  which  states  that  at  a  given  plane 
the  angular  spectrum  and  space  correlation  function  are  a  Fourier  transform  pair.  It  has  already 
been  mentioned  that  the  ground  pattern  <f  (x,  z)  at  the  observation  plane  is  simple  related  toU^fs), 
the  angular  spectrum  based  at  the  aperture.  Thus  the  left  side  of  Eq.  (A-10)  is  essentially  the 
space  correlation  function  at  the  observation  plane.  Also,  it  is  easily  seen  that  the  (received)  an¬ 
gular  spectrum  is  proportional  to  the  illumination  distribution  in  the  aperture  (clearly,  the  radi¬ 
ation  is  limited  to  the  angle  subtended  at  the  receiver  by  the  source).  This  duality  between  angle 
at  one  terminal  and  distance  along  the  plane  at  the  other  terminal  thus  allows  the  inference  of  ei¬ 
ther  Eq.  (A-9)  or  Eq.  (A-10)  from  the  other. 

It  might  also  be  mentioned,  finally,  that  most  experimental  measurement*  consist  of  time, 
rather  than  ensemble,  averages.  This  is  essentially  an  assumption  of  ergodicity,  or  an  assump¬ 
tion  that  the  fields  being  observed  are  being  replaced  in  time  by  other  members  of  the  ensemble, 
so  that  the  one  average  is  equivalent  to  the  other.  We  can  thus  regard  the  abbreviation  (  >  used 
in  this  Appendix  as  representing  a  time  average. 
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APPENDIX  B 

SCATTERING  FUNCTION  FOR  UNIFORM  ROUGH  SPHERE8 


1 

i 


f 


« 
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The  problem  treated  here  is  the  determination  of  the  distribution  in  range  and  doppler  of 
the  radar  return  from  a  uniform  rough  rotating  sphere.  This  quantity  (defined  as  a)  will  depend 
on  the  diameter  of  the  body,  its  rotation  velocity,  axial  tilt  relative  to  the  line  of  sight,  and  the 
carrier  frequency  in  use.  It  will  also  be  a  function  of  the  radar  reflecting  characteristic  <rQ(ip) 

(to  be  defined)  of  the  material  of  which  the  surface  is  composed,  and  we  will  be  interested  in  see¬ 
ing  the  relationship  between  the  radar  return  and  this  quantity  <j  describing  the  surface  properties. 
We  define 

a  =  tilt  of  axis  from  perpendicularity  with  line  of  sight, 

SI  =  rotation  velocity  in  radians  per  second, 
f  =  carrier  frequency, 

a  (ip)  =  reflecting  cross  section  per  unit  area  of  surface  as  a  function  of  tilt 

of  the  radar  line  of  sight  from  normal  (radar  angular  power  spectrum.) 

We  seek  to  find  the  expression  showing  how  these  quantities  affect  the  scattering  function 

<t(t,  f)  =  time  average  radar  energy  returned  in  a  hypothetical  element  of  time 
delay  (range)  dr  wide  at  delay  r  and  df  wide  at  doppler  frequency  f. 

We  shall  assume  that  the  position  of  the  center  of  the  sphere  corresponds  to  the  r-origin  and  that 
the  velocity  of  the  center  of  the  sphere  (assumed  in  straight-line  motion)  corresponds  to  the  f-origin. 

Figure  B-l  depicts  the  rotating  body  cut  by  a  plane  Q  containing  both  line  of  sight  and  axis. 

An  arbitrary  point  P  on  the  surface  is  shown  having  coordinates  (0,y).  Also  shown  is  the  angle  <p 
between  line  of  sight  and  radius  vector  from  the  center  of  the  sphere  to  the  point  P. 

The  first  step  is  to  find  the  projection  on  the  line  of  sight  of  V,  the  velocity  of  P.  First, 


magnitude  of  V,  V  =  O  a  sin  0  (B-l) 

The  projection  of  V  on  the  plane  Q  is 

il  a  sin0  siny  ;  (B-2) 

the  projection  of  this  quantity  in  turn  along  the  line  of  sight  is 

J1  a  sin  0  sin  y  cos  a  ,  (B-3) 

and  the  doppler  frequency  offset  of  the  echo  from  P  is 
2f 

f  =  fi  a  sin  0  sin  y  cos  a  .  (B-4) 

But  the  perpendicular  distance  x  from  point  P  to  plane  Q  is 

x  =  a  sin  y  sin  0  ,  (B-5) 


so 


S2  x  cos  a 


(-  a  <  x  <  a) 


The  range  (delay)  of  point  P  is  clearly  given  by 
,  2a 

t  -  —  cos  <p  =  —  t  , 
c 


(B-6) 


(B— 7) 
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where  for  present  purposes  r'  reads  in  a  positive  direction  toward  the  observer. 

Now  we  have  defined  aQ(<p)  to  be  the  cross  section  per  unit  area  of  target  surface  {not  per 
unit  area  of  incident  wavefront) ,  if  is  at  shown  in  the  figure.  If  we  assume  unit  power  flux  den¬ 
sity  throughout  the  incident  wavefront,  the  total  power  returned  with  delay  r  and  doppler  f  in 
a  range  of  delay  dr  and  doppler  df  is 

<r[ t,  f)  =  (number  of  square  meters  of  wavefront  corresponding  to  the 

element  dfdr)  x  (cross  section  per  unit  area  of  wavefront).  (B-8) 

The  second  factor  is  aJi<p)/coB(p. 

Call  the  first  factor  A.  Figure  B-2  shows  a  head-on  view  of  the  target  (looking  along  the 
line  of  sight  with  the  line  corresponding  to  plane  Q  in  the  up-down  direction).  The  two  regions 
making  up  area  A  are  shown  shaded.  The  width  Ax  (in  meters)  of  the  strip  containing  all  the 
energy  lying  in  unit  bandwidth  at  frequency  f  is  [from  Eq.  (B-6)] 

Ax  =  c/(2fQfl  cos  a)  .  (B-9) 

The  width  of  the  annular  strip  containing  all  the  energy  lying  in  unit  spread  of  delays  at  delay  r' 
we  will  call  Ay.  From  Fig.  B-3, 

Ay  =  j  cos  ip  (B-iO) 

The  area  A  is  given  by 

A=2AxAy/cos4  ,  (B-ll) 

where  6  is  indicated  in  Fig.  B-2,  and 

sine .  V^-Tct1/*)2  •  (**-«> 

o 


Combining  (B-12)  with  (B-ll)  and  (B-7)  with  (B-10)  and  then  substituting  (B-9),  (B-10)  and  (B-ll) 
in  (B-8)  we  have  finally 


o(r',  f)  - 


c2o  [cos_1(CT'/2a)]  ,  , 

- zT{l  -Qs- -  [1-  (cr'/2a)  -  (fc/2foa«  cosc)  )  ' 


lo  for  imaginary  values 

From  Eqs.  (37)  and  (39),  respectively,  power  as  a  function  of  delay  only  is 
(t(t)  =  t  ac  (To(cos"1(cT'/2a)) 
and  power  as  a  function  of  doppler  offset  alone  is 


for  real  values 
(B-13) 

(B-14) 


(B-13)  df 


(B-15) 
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